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INTRODUCTION 

The  Center  for  Thin  Film  Studies  was  formed  in  1986  under  the  University 
m  Research  Initiative  Program,  funded  by  the  Air  Force  Office  of  Scientific  Research. 

This  report  contains  a  sununary  of  the  research  performed  under  this  umbrella 
grant;  separate  sections  cover  work  on  growth  and  characterization  of  thin  films  by 
different  methods,  modeling  of  thin  film  growth,  and  preparation  and  characterization 
of  substrates  for  growth. 

The  task  numbers  have  been  reassigned  from  the  previous  year. 


The  work  described  herein  was  performed  under  Contract  it  F4%20-86-C- 
0123P(X)001. 


r 


TASK  1 
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NUCLEATION  AND  GROWTH  STUDIES 
BY  CONVENTIONAL  AND  ATOMIC  LAYER  EVAPORATION 

Principal  Investigator:  U.  J.  Gibson 

Project  Goal 

The  goal  of  the  project  is  to  study  the  details  of  nucleation  and  film  growth, 
primarily  in  systems  that  lend  themselves  to  deposition  by  different  evaporation 
methods.  The  research  over  the  last  year  had  two  major  thrusts:  1)  study  of  the 
properties  of.  and  manipulation  of  growth  habits  in,  conventional  and  ion-assisted  ZnS 
deposition,  and  2)  detailed  study  of  the  formation  of  initial  layers  of  ZnS  by 
molecular  beam  epitaxy  (MBE)  and  atomic  layer  evaporation  (ALE)  on  substrates 
prepared  using  different  protocols. 

Work  Accomplished 

ZnS  Films 

ZnS  is  an  interesting  material  for  many  reasons,  including  its  high  refractive 
index,  wide  transparency  range,  and  a  large  thermal  nonlinearity  that  is  useful  for 
demonstration  experiments.  It  is  also  an  interesting  system  for  studies  of  growth, 
since  it  is  one  of  the  few  materials  that  form  compressive  thin  films  when  deposited 
at  room  temperature.  This  is  believed  to  occur  because  the  Zn  and  S  dissociate 
during  evaporation,  and  the  individual  mobilities  of  the  elements  dominate  the  growth 
kinetics  at  the  substrate.  It  thus  fails  to  fit  into  the  empirical  structure  zone  model, 
and  may  provide  more  rigorous  tests  of  first  principle  theories  of  film  growth. 

The  properties  of  ZnS  films  were  studied  using  waveguiding  experiments,  x-ray 
diffraction,  electron  microscopy,  and,  in  collaboration  with  G.  Stegeman,  Brillouin 
scattering  and  second  harmonic  generation  measurements.  Rutherford  backscattering 
spectroscopy  (RBS)  was  also  an  invaluable  resource. 

Despite  extensive  development  for  integrated  optics  applications,  waveguiding 
techniques  have  been  used  only  recently  for  analysis  of  the  properties  of  thin  films. 
Single  layers  can  be  analyzed  readily  to  d.;termine  thickness  and  refractive  index,  and 
waveguiding  techniques  provide  an  extremely  sensitive  measure  of  losses  caused  by 
both  absorptive  and  scattering  mechanisms. 
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Structural  and  waveguide  studies 

The  waveguide  loss  effects  of  water  vapor  and  oxygen  partial  pressures  during 
evaporation  were  reported  on  previously.^  Changes  in  the  orientation  and  average 
sizes  of  the  grains  within  the  film  were  correlated  with  these  partial  pressures. 
Those  experiments  indicated  that  the  dominant  loss  mechanism  in  ZnS  waveguides 
was  scattering  from  large  crystallites.  The  crystallite  size  is  strongly  dependent  on 
substrate  temperature,  and  implementation  of  substrate  cooling  has  allowed 
achievement  of  essentially  amorphous  films.  ZnS  crystallite  size  (determined  by 
XRD)  as  a  function  of  deposition  temperature  is  shown  in  Fig.  1.  Using  cooled 
substrates,  reductions  of  waveguide  losses  from  approximately  10  dB/cm  to  less  than 
0.5  dB/cm  in  the  red.  and  from  25  dB/cm  to  less  than  5  dB/cm  in  the  green,  have 
been  achieved.  These  lower  losses  may  permit  further  investigation  and  exploitation 
of  the  nonlinear  properties  of  21nS  in  integrated  optics  applications. 


Figure  1.  ZnS  crystallite  size  versus  the  substrate  temperature  during  deposition. 

The  films  made  at  reduced  temperatures  also  have  reduced  indices  of  refraction, 
as  shown  in  Fig.  2.  The  observed  reduction  for  cooled  films  may  be  caused  by 
formation  of  ZnO  during  deposition,  or  by  a  fibrous  or  columnar  structure  with  a 
large  volume  of  pores.  At  present,  there  is  evidence  for  each  of  these  models.  RBS 
analysis  reveals  a  stoichiometric  Zn  to  S  ratio,  favoring  the  porosity  model.  High- 
power  waveguiding  experiments  also  displayed  long-term  instabilities,  which  may  be 
caused  by  water  desorption;  vacuum  waveguiding  experiments  will  be  undertaken  to 
clarify  these  results.  Alternatively,  neither  electron  microscopy  nor  second  harmonic 
generation  studies  revealed  evidence  of  columnar  growth  for  the  cooled  films. 
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Temperature  of  Substrate  (K) 

Figure  2.  Index  of  refraction  versus  substrate  temperature. 


The  structure  zone  model  of  Movchan  and  Demchishin  predicts  a  changeover 
from  a  columnar  film  to  an  equiaxed  grain  structure  at  a  temperature  of 
approximately  0.45  T„.  where  T,„  is  the  melting  point  of  the  material  being 
deposited.  For  ZnS,  the  melting  temperature  is  1830*>C,  implying  a  transition 

temperature  of  673'’C.  whereas  films  dep<^ited  at  temperatures  as  low  as  30"C 

display  the  dense  structure  expected  above  the  transition.  If  for  Zn  is  used,  the 
transition  would  be  expected  at  approximately  38°C.  It  is  more  difficult  to  assign  a 
value  for  the  sulfur  transition,  since  sulfur  clusters  of  different  sizes  (2,  4,  8)  are 
readily  formed.  For  S,.  the  transition  temperature  would  be  about  -98°C. 

Obviously,  assigning  these  individual  values  is  less  than  entirely  appropriate.  Work 
on  use  of  mobilities  and  chemical  potentials  to  describe  the  competing  mechanisms  in 
growth  of  these  films  is  underway  in  Task  4. 

Index-of-refraction  data  for  ion-assisted  deposition  at  low  temperatures  are 
shown  in  Fig.  2.  In  this  case,  the  reduction  of  index  is  much  less  masked. 
Chemical  changes  induced  by  the  bombardment  compete  with  densification  effects. 
RBS  indicates  the  presence  of  significant  amounts  of  argon  in  the  films,  as  well  as  a 

sulfur  deficiency.  These  results  suggest  use  of  HjS  as  the  bombarding  species  in 

future  work. 


Stress  changes 

A  significant  change  in  the  stress  properties  of  ZnS  has  also  been  observed  on 
cooling,  with  the  stress  going  from  compressive  at  room  temperature,  to  strongly 
tensile  at  -lOO^C.  At  about  -S(y*C.  the  stress  is  minimized,  and  much  thicker  layers 
than  normally  achievable  are  readily  attained.  This  may  have  significant  impact  on 
the  deposition  of  narrow-band  filters,  where  stress  effects  often  limit  the  choice  of 
materials  and  the  maximum  spacer  layer  thicknesses.  Compressive  stress  problems 
are  particularly  severe  in  ion-beam-sputtered  coatings;  the  results  of  this  work 
suggest  the  use  of  a  cooled  substrate  for  these  depositions.  Further  studies  on  ion- 
assisted  deposition  are  aimed  at  understanding  the  potential  of  the  combination  of 
energetic  bombardment  with  a  reduced  temperature  substrate. 

Second  harmonic  studies 

Earlier  results  on  the  sensitivity  of  the  second  harmonic  generation  (SHG) 
process  to  columnar  structure^  led  to  similar  studies  in  ZnS.  a  material  in  which 
second  harmoiiic  generation  is  a  symmetry-allowed  process.  However,  ambient 
temperature  depositions  often  resulted  in  preferential  crystal  orientation  of  the  ZnS, 
which  dominated  the  column-induced  anisotropy  in  the  second  harmonic  signal.  Ion 
bombardment  during  deposition  was  used  in  an  attempt  to  quench  the  preferential 
orientation;  this  merely  introduced  additional  texture  into  the  SHG  signal.  Recently, 
ZnS  was  deposited  at  low  temperature  (-SO^Q.  with  a  substrate  incidence  angle  of  30 
degrees.  Conventional  wisdom  implies  that  this  would  result  in  a  highly  columnar 
microstructure,  with  a  clear  anisotropy  in  SHG.  It  was  found,  however,  that  the 
SHG  signal,  while  still  strong,  gave  no  indication  of  the  formation  of  columns.  The 
measurement  geometry  is  shown  in  Fig.  3a.  The  results  of  SHG  measurements  for 
1)  ZnS  deposited  at  room  temperature  (without  preferential  orientation),  2)  ion-assisted 
ZnS  deposition,  with  evidence  of  additional  symmetry  breaking  caused  by  the  beam. 
3)  the  room-temperature  film  measured  with  9*0,  and  4)  a  cooled  film  measured  at 
the  same  9  are  shown  in  Fig.  3b.  For  0-0.  there  was  no  signal  from  the  cooled 
film,  indicating  symmetry  in  the  plane  of  the  film.  All  films  were  deposited  with  a 
vapor  incidence  angle  of  30  degrees  to  the  substrate  normal. 

Brillouin  Samples 

ZnS  samples  were  also  prepared  for  Brillouin  adhesion  studies,  as  described  in 
the  Task  3  report. 


UHV  Studies 


Investigations  of  the  nucleation  and  growth  behavior  of  ZnS  films,  deposited  in 
I  ultra-high  vacuum  (UHV)  from  a  molecular  source  and  from  reactive  materials,  also 

have  been  performed. 

ZnS  films  were  deposited  onto  oxidized  single-crystal  Si  under  UHV  conditions 
for  comparison  with  conventionally  evaporated  films.  However,  the  large  grain  size 
I  and  high  internal  stress  resulted  in  high  scatter  and  short  lifetimes.  Epitaxial  films 

have  been  deposited  onto  Si.  and  investigated  by  a  number  of  techniques.  One  such 
film  was  mechanically  damaged,  thus  inducing  microcracks,  and  submitted  for 
Brillouin  testing,  as  described  in  the  Task  3  report. 

I  Studies  of  the  nucleation  and  substrate  interaction  behavior  were  undertaken  to 

elucidate  the  role  of  cleaning  procedures  in  ZnS  growth  on  Si  wafers,  and  to 
understand  the  critical  first-layer  growth  in  ALE  of  ZnS.  In  particular,  the  role  of 
exposure  of  the  Si  to  HjS  during  heating,  prior  to  deposition,  was  of  interest.  It  was 
^  found  that  the  presence  of  HjS  did  not  significantly  affect  the  temperature  at  which 

RHEED  spots  were  first  observable,  but  it  did  significantly  affect  the  surface 
impurity  concentrations,  as  observed  by  XPS.  Lower  concentrations  of  both  C  and  O 
were  observed  with  HjS  treatment. 

>  Depositions  of  thin  Zn  layers  onto  Si  substrates  prepared  with  the  HjS  treatment 

gave  a  strong  indication  of  a  preference  for  a  layer-by-layer  growth  habit,  whereas 
high-aspect  ratio  islands  were  observed  without  the  HjS.  Angle-resolved  XPS  data 
were  taken  on  the  ratio  of  the  Zn  to  substrate  (Si)  signal,  and  analyzed  according  to 
^  the  Fadley  model  for  the  two  cases,  as  shown  in  Figs.  4a  and  4b.  For  the  substrates 

exposed  to  HjS  after  cleaning  and  heating,  the  Zn  formed  a  thin  layer  over 
approximately  80%  of  the  substrate  surface.  The  incomplete  coverage  is  attributed  to 
oxygen  impurities  blocking  sulfur  bonding  at  the  interface.  For  the  substrates 

^  exposed  to  Zn  following  attainment  of  a  RHEED  pattern,  but  with  no  HjS  exposure, 

the  slope  of  the  curve  of  the  Zn:Si  ratio  versus  theta  is  discontinuous,  indicating  that 
there  are  large  aspect  ratio  islands  of  Zn  forming  in  preference  to  continuous  layers. 

^  Plans 

In  the  coming  year,  additional  work  will  be  performed  to  determine  the 
structure  of  cooled  ZnS  films,  using  both  optical  and  electron  microscopy  techniques, 
I  and  to  investigate  the  simultaneous  use  of  ion  bombardment  and  substrate  cooling  for 

production  of  low-loss,  stable  ZnS  material. 
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Figure  4.  Zn:Si  XPS  data:  (a)  substrate  exposed  to  HjS  -  the  squares  are  data  points, 
the  dashed  line  is  the  Fadley  fit  for  50%  coverage:  (b)  substrate  not  exposed  to  H^S  -  the 
crosses  are  data  points,  no  Fadley  fit  possible. 

p  It  is  hoped  that  comparison  of  these  and  SHC  studies  on  a  material  such  as 

MgFj.  which  evaporates  molecularly.  can  be  made  to  determine  the  importance  of  the 
evaporative  dissociation  of  ZnS  in  determining  its  growth  patterns. 

The  nucleation  studies  will  also  be  pursued  to  illuminate  the  formation  of 
k  subsequent  layers  in  ALE  deposition.  The  treatment  chamber  will  be  used  for  these 

studies,  and  the  MBE  system  will  be  used  for  thick  layers  and  deposition  of  other 
materials. 
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Experiments  are  also  planned  to  test  some  of  the  theoretical  predictions  made  in 
the  other  Tasks  regarding  the  development  of  columnar  structure.  In  particular.  A1 
films  will  be  studied  to  test  the  dependence  of  column  size  on  deposition  rate. 

Personnel 

U.  Gibson 
J.  Watanabe 
J.  Ruffner 
B.  Williams 
F.  Suits 


Collaborations 

Macleod:  XPS  analysis  of  oxynitride  films 

Leavitt  and  McIntyre:  RBS 

Stegeman  and  Seaton:  Brillouin  and  SHG  studies 

Lichter:  Columnar  growth  model,  chemical  potential  model 

Lindley:  Preparation  of  opal,  thin  film  samples  for  EM  and  STM  studies 
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ION-ASSISTED  DEPOSITION  AND  MODELING 

Principal  Investigators:  Bertrand  G.  Bovard  and  H.  Angus  Macleod 

Project  Goals 

Investigate  the  use  of  ion-assisted  deposition  (lAD)  for  modification  of  film 
structure,  and  determine  its  utility  for  improving  the  performance  of  optical  films. 

Develop  a  model  of  film  growth  that  includes  the  range  of  effects  present  in  an 
actual  coating  chamber. 

Work  Accomplished 

Ion-Assisted  Deposition 

Work  on  the  ion-assisted  deposition  (lAD)  of  metals  has  continued.  Some  of  the 
results  quoted  here  were  reported  previously,  but  are  included  for  the  sake  of  clarity 
and  continuity.  Metal  films  are  important  in  a  wide  range  of  optical  coatings  and 
the  effects  of  lAD  on  the  optical  properties  is  therefore  of  some  importance. 

From  the  point  of  view  of  optical  properties,  the  results  of  argon  lAD  of  silver 
are  disappointing.  The  real  part  of  the  refractive  index  rises  and  the  imaginary  part 
falls.  This  measurement  required  great  precision  to  be  certain  of  the  consistency  of 
the  small  shifts  observed  in  the  constants.  Some  wide-spectral-range  ellipsometry  was 
performed  for  this  Task  by  Professor  Vedam’s  group  at  Pennsylvania  State 
University,  but  the  variation  of  optical  constants  is  small  and  so  a  more  sensitive 
technique  for  a  single  wavelength,  involving  the  measurement  of  a  surface  plasmon 
resonance,  has  also  been  employed.  The  results  show  a  clear  increase  in  the  real 
part  of  the  refractive  index  and  a  reduction  in  the  imaginary  part,  both  of  which 
correlate  well  with  the  rate  of  momentum  transferred  to  the  growing  film.  DC 
resistivity  increases  with  bombardment,  and  a  change  in  lattice  constant  occurs  that  is 
consistent  with  a  shift  from  tensile  stress  to  compressive  stress  with  increasing 
bombardment.  Rutherford  backscattering  (RBS)  shows  implanted  argon  increasing 
linearly  with  the  bombarding  momentum  transfer  rate.  TTie  shifts  in  refractive  index 
and  conductivity  are  consistent  with  the  effects  expected  from  the  implanted  argon. 
Aluminum  is  affected  in  a  similar  way.  although  the  optical  constant  variation  is 
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much  less  linear  and  so.  for  light  bombardment,  the  optical  constants  are  less 
affected.  The  work  is  summarized  in  the  attached  Applied  Optics  submission. 

Some  of  the  implications  are  clear.  Silver  is  a  material  often  used  for  high- 
power  infrared  coatings,  usually  with  some  dielectric  overcoats  to  boost  reflectance. 
Ion  assistance  to  improve  adhesion  and,  after  Sainty  et  al..^  to  increase  chemical 
resistance  would  seem  a  very  attractive  technique.  Our  findings  suggest  that  the 
effects  of  the  ion  bombardment  are  not  all  beneficial  and  that  the  reduction  in 
performance  must  be  taken  into  account.  Improved  adhesion  without  the  overall 
reduction  in  optical  performance  may  be  possible  if  ion  assistance  is  limited  to  only 
the  very  early  stages  of  film  deposition.  Work  is  in  progress  on  that  aspect,  with 
some  encouraging  results. 

lAD  of  Dielectrics.  The  problems  of  dielectric  deposition  with  ion  assistance  are 
quite  different.  Here  the  implantation  of  the  bombarding  species  appears  to  be  quite 
benign.  Experiments  involving  high-temperature  baking  of  films  indicate  that  even 
implanted  argon  is  firmly  embedded  and  shows  no  tendency  to  evolve.  When  the 
ions  are  reactive  (e.g.,  oxygen  or  nitrogen),  beneficial  chemical  effects  can  result. 

The  fluorides  represent  a  very  important  class  of  thin-film  materials,  with  a 
much  wider  transmittance  region  than  other  materials,  and  yet  their  iisefulness  in 
thin-film  coatings,  especially  in  the  infrared,  has  been  limited,  especially  by  their 
generally  high  intrinsic  tensile  stress.  Hie  studies  of  the  effects  of  ion  assistance  on 
the  properties  of  the  fluorides  of  the  lanthanide  series  have  continued.  Ion-assisted 
deposition  is  of  great  benefit  in  reducing  the  tensile  stress,  but  it  also  damages  the 
films  by  affecting  their  stoichiometry.  This  damage  causes  optical  absorption,  but  the 
damage  can  be  limited  and  also  can  be  partially  healed,  often  by  incorporation  of 
oxygen.  Recent  work  has  centered  on  the  fluoride  of  samarium  (see  Fig.  1)  and  the 
nature  of  the  vacancies  left  by  the  bombardment. 

The  usefulness  of  ion-assisted  deposition  when  the  bombarding  species  is  a 
reactive  gas  was  mentioned  above.  Oxygen  has  long  been  used  in  the  optical  coating 
of  refractory  oxides.  A  background  of  oxygen  in  the  coating  chamber  improves  the 
degree  of  oxidation  of  the  films  and  decreases  their  optical  absorption.  Bombardment 
of  growing  metal  films  with  oxygen  ions  can  result  in  very  good  oxide  layers. 
Nitrogen  is  another  compound-forming  gas  which  lacks  sufficient  reactivity  to  have 
been  a  useful  background  gas  in  deposition.  The  nitrides  have,  therefore,  never  been 
used  successfully  in  classical  evaporation  processes,  although  some  success  has  'oeen 
reported  using  sputtering.  This  research  has  confirmed  that  nitrogen  ion 
bombardment  during  deposition  is  a  viable  technique  for  producing  nitrides.  Recent 
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Figure  1.  The  measured  refractive  index  of  SmFi  in  the  visible  and  near  UV. 


progress  in  this  area  has  been  in  the  study  of  the  aluminum-oxygen-nitrogen  system. 
A  range  of  compounds  can  be  produced,  with  refractive  indices  from  1.67 
(corresponding  to  pure  oxide)  to  1.8S  (corresponding  to  the  nitride),  by  evaporating 
metallic  aluminum  under  various  residual  atmospheric  conditions.  The  conditions 
finally  chosen  for  the  demonstration,  described  in  the  appended  manuscript,  involve  a 
background  gas  of  oxygen  together  with  nitrogen  ion  bombardment.  This  particular 
arrangement  has  the  benefit  of  avoiding  the  corrosive  effect  of  oxygen  on  the  hot 
cathode  of  the  Kaufman  sources,  enormously  increasing  their  lifetimes,  and  also 
increasing  the  reactivity  of  nitrogen  so  that  compounds  containing  it  can  be  produced. 
The  range  of  refractive  indices  as  a  function  of  composition  is  shown  both  for 
nitrogen  bombardment  and  bombardment  by  a  nitrogen-oxygen  mixture  (Figs.  2 
and  3).  As  a  test  of  the  process  and  a  demonstration  of  its  potential,  a  IS-period 
rugate  filter  has  been  constructed.  Each  period  was  made  up  of  10  subperiods  so 
that  the  entire  filter  consists  of  150  separate  layers  (Figs.  4  and  5).  The  variation  of 
properties  throughout  the  entire  filter  was  achieved  by  varying  the  background 
oxygen  pressure  and  keeping  the  ion-gun  conditions  and  the  aluminum  evaporation 
conditions  constant.  Nitrogen  was  the  bombarding  gas.  The  only  variable,  therefore, 
for  the  complete  operation  was  the  background  oxygen  pressure.  The  packing 
density  achieved  was  good  because  no  vacuum-to-air  shift  was  detected  It  should 
be  noted  that  the  oxygen  partial  pressure  was  manually  controlled,  which  limited  the 
complexity  of  the  coating  that  could  be  attempted.  More  complex  and  demanding 
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pressure  was  varied  from  0.4xl0‘*  to 
L4xW'*  mbar,  while  a  nitrogen  ion  beam 
was  bombarded  on  a  growing  Al  film.  Ion 
energy  was  1250  eV  and  current  density 
was  40  pAfcn?, 


Figure  4.  A  step- index  profile  of 
AlO^Ny  films  obtained  from  Fig.  I  for  a 
rugate  filter.  One  cycle  consists  of  10 
homogeneous  step  layers  and  has  a 
169  nm  physical  thickness.  A  sinusoidal 
curve  is  fitted  to  the  step-index  profile. 
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Figure  3.  Variation  of  refractive  index 


and  extinction  coefficient  of  AlO^Ny  thin 
550  nm.  A 


films  at  550  nm.  A  mixture  of 
nitrogen! oxygen  gases  was  varied 
(OiSfxlO-^  to  (8/0)xl0'^  mbar.  A 
backfilled  oxygen  pressure  was 
1.2x10'*  mbar,  while  the  mixed  ion  beam 
was  bombarded  on  a  growing  Al  film. 


Figure  5.  Reflectance  and  transmittance 
of  the  real  and  simulated  rugate  filters. 
Fifteen  cycles  were  employed. 


coatings  may  be  possible  using  automatic  control  techniques.  These  investigations  are 
now  centered  on  arriving  at  an  understanding  of  the  nature  of  the  aluminum-oxygen- 
nitrogen  material  and  its  properties. 

Modeling 

Thermo-Optical  and  Nonlinear  Effects  in  Optical  Multilayers.  The  construction  of  a 
thermo-optical  model  began  in  connection  with  thermally  induced  bistability,  but  the 
code  that  has  been  developed  is  quite  general.  It  is  able  to  investigate  the  thermal 
effects  of  laser  light  on  the  optical  properties  of  optical  coatings  made  of  metal  or 
dielectrics.  Work  on  the  code  has  continued  under  this  contract.  The  code  takes 
into  account  the  electric-field  repartition  inside  the  coating  to  calculate  the  amount  of 
energy  absorbed  by  the  coating.  This  absorption  leads  to  an  increase  in  temperature, 
which  modifies  the  refractive  indices  of  the  materials  involved  in  the  multilayer 
coating.  The  coupling  of  heat-flow  equations  with  optical  properties  permits  the 
monitoring  of  the  changes  in  the  electric  field,  the  transmittance,  and  the  reflectance 
during  laser  irradiation.  This  simulation  has  been  used  to  investigate  the  bistability 
of  laser-irradiated  Fabry-Perot  etalons  made  of  2^  or  ZnSe.’  Until  that  time,  the 
nonlinearity  included  in  the  model  was  only  temperature  dependent. 

A  recent  paper  published  in  Physical  Review  Letters’  triggered  interest  here  in 
nonlinear  effects  occurring  in  structures  made  of  electric-field-dependent  nonlinear 
materials.  In  the  cited  work,  a  multilayer  dielectric  mirror  was  subjected  to  laser 
irradiation  at  a  wavelength  close  to  its  stop  band.  By  allowing  the  refractive  index 
to  vary  with  the  internal  electric  field,  those  authors  were  able  to  demonstrate  the 
existence  of  several  transmitting  states  according  to  the  incident  power  and  the  degree 
of  nonlinearity. 

The  existing  code  used  here  for  thermal  simulation  was  modified  to  include  the 
nonlinear  electric-field  dependence  of  the  refractive  index,  and  some  very  interesting 
bistability  loops  have  been  computed  that  are  in  agreement  with  the  referenced 
paper.  This  first  study  demonstrates  the  possibility  of  switching  a  reflecting 
component  to  a  transmitting  state.  Rugate  filters  are  mirrors  with  very  narrow  stop 
bands,  making  them  very  similar  to  this  example. 

An  investigation  of  the  consequences  of  nonlinearities  in  the  refractive  indices 
of  the  materials  involved  in  a  multilayer  structure  was  carried  out  for  the  case  of  a 
very  narrow  reflection  band  mirror.  The  design  chosen  for  this  study  is  simply  a 
quarterwave  stack  of  formula  (HL)’®®H.  where  H  is  a  quarterwave  of  material  of 
refractive  index  1.51  and  L  represents  a  quarterwave  of  material  of  refractive  index 
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1.50.  The  wavelength  of  investigation  was  the  central  wavelength  of  the  stop  band 
of  the  mirror.  This  simple  design  was  chosen  because  the  corresponding  reflectance 
properties  are  close  the  the  properties  of  a  real  rugate  filter  (Fig.  6.) 

TRANSMITTANCE  (%)  (SL)''200  H 


WAVELENGTH 

Figure  6.  Transmittance  versus  wavelength  for  a  (HL)^H  mirror. 


In  the  computations,  both  materials  have  an  electric-field  dependence  of  the 
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E  -  eo(l  +  ^  )  (2) 

with 

X  -  XoEj  .  (3) 

Equation  (2)  is  used  in  conjunction  with  Eq.  (3)  to  integrate  both  the  effect  of 
the  nonlinearity  of  the  material  and  the  effect  of  laser  intensity  in  one  parameter. 

In  these  three  equations,  (q  is  the  dielectric  permittivity  of  the  material,  \  is  the 
electric-field  dependence  of  the  permittivity,  E  is  the  electric-field  modulus.  Eg  is  the 
incident  electric-field  modulus,  and  E„  is  the  electric-field  modulus  normalized  to  a 
1  W/m*  incident  power  density. 

As  the  power  or  the  nonlinearity  is  increased,  corresponding  to  an  increase  in 
the  parameter  X,  the  transmittance  goes  through  a  hysteresis  loop  (Fig.  7),  allowing  the 
device  to  become  transparent  once  a  threshold  value  is  passed,  and  to  remain 
transparent  as  the  power  is  reduced. 
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Figure  7.  Bistability  of  a  mirror  (HL)^H. 


This  bistable  behavior  is  very  likely  to  exist  in  rugate  filters  if  the  parameter  \ 
can  reach  a  value  higher  than  the  threshold.  It  is  also  important  to  realize  that  laser 
irradiation  will  cause  heating  by  residual  absorption  and  other  losses,  causing  a 
variation  in  the  refractive  index  that  results  from  its  temperature  dependence. 
Therefore,  even  if  the  X  parameter  does  not  reach  the  threshold  value,  the  influence 
of  the  heating  may  still  cause  a  bistable  behavior  if  the  thermal  diffusion  allows  it; 
this  may  well  be  the  case  in  very  thick  coatings. 

It  is  hoped  that  this  investigation  can  be  pursued  by  introducing  into  the  code 
refractive  index  profiles  that  are  closer  to  rugate  filter  designs,  to  study  their 
behavior  under  energetic  laser  fluxes,  in  terms  of  both  electric-field  and  temperature 
effects. 

Personnel 

H.  A.  Macleod 

B.  G.  Bovard 

C.  K.  Hwangbo 
M.  R.  Jacobson 
J.  Lehan 

L.  Lingg 
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Proposal 

IN-SITU  STRESS  MEASUREMENTS 
Principal  Investigator:  B.  G.  Bovard 
Introduction 

Thin  films  deposited  by  thermal  evaporation  suffer  from  high  levels  of  stress, 
which  can  be  tensile  or  compressive,  depending  on  the  materials  and  the  growth 
process  used.  For  many  materials,  a  high  level  of  stress  limits  the  maximum 
thickness  possible,  and  is  an  important  cause  of  mechanical  failure. 

Stress  is  usually  divided  into  two  components:  the  intrinsic  part,  resulting  from 
the  growth  process  itself,  and  the  thermal  part,  resulting  from  the  differing  expansion 
coefficients  of  the  substrate  and  film  materials. 

Project  Goal' 

Monitoring  and  controlling  stress  during-  deposition  are  important  steps  toward 
understanding  the  underlying  mechanisms  that  drive  stress  during  film  growth. 
Therefore,  we  propose  to  install  an  in-situ  stress  monitor  in  one  of  our  vacuum 
chambers.  The  system  will  be  based  on  the  cantilever  beam  technique,  in  which  the 
displacement  of  the  end  of  a  thin  substrate  is  measured  to  calculate  the  stress.’ 

A  number  of  methods  of  measuring  the  displacement  are  possible  (e.g., 
capacitance  measurements).  Optical  measurement  of  the  displacement  using  a  laser 
deflection  technique  is  the  primary  candidate,  sdthough  other  techniques  will  be 
considered. 

Several  materials  will  then  be  deposited,  including  dielectrics  such  as  fluorides, 
which  usually  suffer  from  very  high  tensile  stress,  and  the  process  monitored. 

Reference 

1.  R.  W.  Hoffman,  "Mechanical  properties  of  non-metallic  thin  films,"  in  Physics  of 
Non-Metalllc  Thin  Films,  C.  H.  S.  Dupuy  and  A.  Cachard.  eds.  (Plenum  Press, 
New  York,  1974). 
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TASK  3 

STRUCTURE  MODinCATION  BY  ION-ASSISTED  DEPOSITION 
(ION  BEAM  ANALYSIS  OF  THIN  FILMS) 

Principal  Investigators:  J.  A.  Leavitt  and  L.  C.  McIntyre,  Jr. 

Project  Goals 

The  goals  of  this  program  are  threefold:  1)  to  provide  accurate  characterization 
of  thin  films  by  MeV  ion-beam  analysis  (IBA)  for  Thin  Film  Center  scientists  and 
other  collaborators;  2)  to  develop  new  and  improve  existing  analysis  techniques;  and 
3)  to  expand  the  capabilities  of  the  IBA  facility  itself. 

^He'*'  analysis  ion  beams  with  energies  in  the  l.S-  to  S.O-MeV  range  from  a 
6-MV  Van  de  Graaff  accelerator  are  used  to  determine  film  thickness,  stoichiometry, 
concentration  profiles,  and  crystalline  structure,  by  means  of  Rutherford 
backscattering  spectroscopy  (RBS),  high-energy  backscattering,  and  channeling.  A 
new  accelerator  source  is  being  developed  to  depth  profile  hydrogen  in  thin  films, 
using  ‘*N++  analysis  beams  with  energies  between  6  and  10  MeV,  and  utilizing  the 
nuclear  reaction  ‘H(‘*N,a7)‘*C. 

Summary 

During  the  second  year  of  the  project,  the  IBA  facility  continued  to  provide 
acciirate  characterization  of  thin  films  for  Center  scientists  (215  films  analyzed  during 
11/1/87  to  9/19/88)  and  other  collaborators  (239  films  analyzed  during  the  same 
period). 

Efforts  to  improve  techniques  have  included:  1)  completion  of  attempts  to  use 
heavier-than-*He  analysis  beams  to  improve  mass  and  depth  resolution  in 
backscattering  spectrometry;  results  indicate  no  major  improvements  with  the  heavier 
beams;  2)  the  successful  depth  profiling  of  P  (below  the  1%  level)  in  SiOj,  using  a 
nuclear-reaction  analysis  technique  involving  observation  of  protons  ejected  by 
incident  *He++;  3)  the  acquisition  of  data  to  measure  non-Rutherford  cross  sections 
for  *He  incident  on  light  elements  (C  to  Si)  in  the  1.8-  to  5-MeV  energy  range,  so 
the  higher  energy  beams  can  be  used  to  reduce  the  overlap  of  signals  from  light 
elements.  Preprints  describing  all  three  efforts  are  available. 
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Because  of  unexpected  difficulty  with  extraction  of  the  desired  doubly  ionized 
species  N'*"*'  and  He'^'*'.  the  testing  and  installation  of  the  new  Penning  ion  source  has 
not  proceeded  as  rapidly  as  expected.  However,  usable  beams  of  these  species  were 
obtained  recently  on  both  60-Hz  and  400-Hz  lab  test  benches;  the  new  source 
assembly  will  be  installed  in  the  Van  de  Graaff  terminal  shortly. 

Work  Accomplished 

Film  CharcKterization 

During  the  period  11/1/87  to  9/19/88,  a  total  of  454  backscattering  analyses 
were  provided.  Complete  analyses  (with  uncertainties)  of  the  data  were  reported  in 
most  cases;  that  is,  film  stoichiometry  (<  ±1%),  thickness  in  atoms/cm‘  (±3%),  and 
impurity  profiles  were  reported.  In  the  majority  of  analyses,  beam  energies  greater 
than  3.5  MeV  were  used  to  reduce  the  overlap  of  signals  from  different  elements;^'’ 
frequently,  the  IBA  facility  deconvolution  programs’  were  applied  to  separate  the 
remaining  overlaps.  In  particular,  for  Center  scientists,  the  following  analyses  were 
performed: 

1)  175  films  for  the  Macleod  group.  The  majority  of  these  films  (87)  were 
fluorides  of  Sm,  Gd,  Mg.  La.  Ce.  Ca,  Lu.  Dy,  Pr.  Tm.  and  Eu.  measured  in 
connection  with  the  use  of  ion-assisted  deposition  for  modification  of  film  structure.**’ 
Measurements  on  37  AlOy  and  AlOxNy  films  were  provided  in  a  similar  study.® 
Other  films  analyzed  were  A1  (26).  Ag  (6).  TiN/TiW  (13),  HfOy  and  SiN  (2). 

2)  Twenty-seven  films  for  the  Gibson  group.  Twelve  ZnS  films  were  emalyzed 
for  their  investigation  of  atomic  layer  evaporation  deposition;  other  films  analyzed 
were  MgF  (6).  ZrO,  (5),  TiO,  (4). 

3)  Four  Langmuir-Blodgett  films  for  the  Stegeman  group. 

Techruque  Improvement 

The  following  improvements  have  been  undertaken: 

1)  The  investigation  into  the  use  of  heavier-than-*He  analysis  beams  to 
improve  mass  and  depth  resolution  has  been  completed.  The  results,  with  C,  N,  Ar. 
and  Kr  analysis  beams,  indicate  that  no  major  improvements  in  mass  and  depth 
resolution  occur  with  use  of  these  beams.*®  at  least,  in  the  analysis  beam  energy 
range  I  to  5  MeV. 

2)  Potassium  was  successfully  depth  profiled  (below  the  1%  level)  in  SiOj, 
using  the  nuclear-reaction  analysis  technique  involving  observation  of  ejected  protons 
for  an  incident  *He  analysis  beam;"  the  method  depends  on  narrow  resonances  in  the 


reaction  cross  section  for  ’‘P(*He,p)’*S.  This  technique  was  used  because  the 
backscaitering  signal  from  Si  overwhelms  that  from  P.  so  backscattering  cannot  be 
used  to  depth  profile  P  in  SiOj. 

3)  A  massive  project  to  measure  non- Rutherford  cross  sections  for  *He 
incident  on  light  elements  (C.  N,  O,  F.  Mg.  Al,  and  Si)  in  the  1.8-  to  5.0-MeV 
energy  range  has  almost  been  completed.  The  data  from  about  500  backscattering 
runs  are  being  analyzed  currently.  The  He-C  data  analysis  is  finished;^’  the  result  is 
shown  in  Fig.  1.  This  result,  and  similar  results  for  the  elements  mentioned  above, 
will  allow  use  of  higher  energy  analysis  beams  to  reduce  the  overlap  of 
backscattering  signals  from  different  light  elements.  Note  that  these  cross  sections 
must  be  measured;  they  cannot  be  calculated,  as  is  the  case  when  the  cross  sections 
are  Rutherford. 

Facility  Development 

A  new  Penning  source  assembly  is  being  developed  for  use  in  the  6-MV  Van 
de  Graaff  to  allow  use  of  the  nuclear  reaction  ‘H(**N.a7)**C  to  depth  profile 
hydrogen  in  thin  films.  Unfortunately,  the  isolated  resonance  in  the  reaction  cross 
section  intended  for  use  occurs  for  a  ‘*N  analysis  beam  energy  of  6.35  MeV.  Hence, 
**N++  beams  must  be  used  to  cover  the  6  to  9  MeV  region  necessary  to  do  the  depth 
profiling. 

The  remainder  of  the  parts  (velocity  selector  and  Einzel  lens)  needed  to  complete 
the  source  assembly  were  received  on  1/19/88.  The  entire  assembly  underwent 
60-Hz  bench  tests,  and  some  time  was  required  to  extract  clean,  doubly  ionized 
beams.  Eventually  it  was  determined  how  to  tune  the  assembly  for  doubly  ionized 
beams.  Typical  results  are  shown  in  Fig.  2;  note  the  presence  of  usable  beams  of 
N++  (2.5  it\)  and  He++  (4.2  p/K). 

The  support  platform  has  been  constructed  for  insertion  of  the  assembly  into  the 
Van  de  Graaff  terminal  and  400-Hz  bench  tests  are  being  run.  Results  similar  to 
those  shown  in  Fig.  2  have  been  obtained,  except  for  the  presence  of  an  additional 
beam  signal  of  unknown  origin. 

The  unexpected  difficulty  with  extraction  of  the  doubly  ionized  beams  has  caused 
the  testing  and  installation  of  the  new  source  to  fall  behind  schedule;  however,  the 
new  source  assembly  is  expected  to  be  in  the  terminal  shortly. 
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Figure  I.  Measured  relative  cross  sections  for  I7(P  backscattering  of  ^He**  ions  by 
carbon  nuclei  for  *He**  incident  at  energies  between  1.5  and  5.0  MeV;  for  relative  cross 
sections  equal  to  1.00,  the  scattering  is  Rutherford:  that  is.  Coulomb's  law  holds.  These 
cross  sections  are  non- Rutherford  for  ^He  energies  above  2  MeV;  previous  measuurements 
of  *He-C  scattering  in  this  energy  range  are  not  suitable  for  use  with  the  backscattering 
technique.  The  strong  non- Rutherford  resonance  at  4260  keV  is  particularly  useful  for 
enhancing  the  sensitivity  for  depth  profiling  carbon  in  heavy-element  thin  films. 
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Figure  2.  Ion  beams  of  nitrogen  and  helium  from  the  new  Penning  source  assembly 
obtained  during  bench  tests  of  this  source.  Note  the  presence  of  doubly  charged  beams: 
2.5  fiA  of  N**  and  4.2  pA  of  He**.  Successful  operation  of  this  source  in  the  6  MV  Van 
de  Graaff  terminal  will  provide  He**  analysis  beams  for  backscattering  measurements  with 
energies  of  at  least  10  MeV,  and  will  allow  depth  profiling  of  hydrogen  in  thin  films  using 
the  nuclear  reaction  ^Hf^^N , ay J^^C. 


Ion  Beam  Analysis  Seminar 

A  10-hour,  no-credit,  no-grade  seminar  on  MeV  IBA  techniques  was  taught  by 
J.  A.  Leavitt  during  the  spring  semester  of  1988.  The  purpose  of  the  seminar  was 
to  teach  collaborators  how  to  analyze  the  backscattering  data  produced  by  this  group. 
Seminar  attendees  were  from  the  UA  Departments  of  Physics.  Chemistry,  Materials 
Science,  Electrical  Engineering,  and  Optical  Sciences,  as  well  as  from  local  industry 
(IBM).  Initial  attendance  was  28  persons;  this  dropped  to  21  by  the  end  of  the 
seminar. 

Plans 

Plans  for  the  next  year  include:  1)  continuing  to  provide  thin  film 

characterization  for  collaborators;  2)  completing  the  400-Hz  Penning  assembly  bench 
test  and  installing  the  source  assembly  in  the  accelerator;  testing  the  source  and 
making  any  necessary  alterations;  depth  profiling  hydrogen  in  a  test  sample  (TiH); 
and  commencing  routine  analysis  of  H  in  thin  films;  3)  analyzing  the  non-Rutherford 
cross-section  data  for  He  on  N,  O,  F,  Mg.  Al.  and  Si  in  the  1.8-  to  S.O-MeV  energy 
region  and  publish  the  results;  measuring  non-Rutherford  cross  sections  for  He  on 
light  elements  in  the  5-  to  10-MeV  energy  region. 

Personnel 

J.  A.  Leavitt 

L.  C.  McIntyre,  Jr. 

P.  Stoss  (partially  supported  through  Characterization  Facility) 

M.  D.  Ashbaugh  (10/1/87-4/1/88) 

J.  G.  Oder  (4/1 /88-present) 

Collaborations 

Macleod  group  (UA  Optical  Sciences):  ion-assisted  modification  of  film  structure, 
fluorides  and  oxynitrides. 

Gibson  group  (UA  Optical  Sciences):  atomic  layer  deposition  of  ZnS. 

Stegeman  group  (UA  Optical  Sciences):  Langmuir-Blodgett  films. 

Falco  group  (UA  Physics):  high  superconducting  films;  MBE  superlattices. 
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TASK  4 

A  NOVEL  TECHNIQUE  FOR  QUANTIFYING 
THE  MECHANICAL  PROPERTIES  OF  THIN  FILMS 

Principal  Investigators:  G.  I.  Stegeman  and  C.  T.  Seaton 

Project  Goals 

The  principal  objective  of  this  project  is  to  develop  Brillouin  scattering  as  a 
technique  for  measuring  the  elastic  properties  of  interesting  thin  films  and  interfaces. 
A  secondary  objective  is  to  use  second  harmonic  generation  (SHG)  as  a  tool  for 
studying  thin  films.  The  actual  work  performed  can  be  subdivided  into  these  two 
areas: 

a)  elastic  properties  of  films  and  their  modification; 

b)  SHG  of  thin  films. 

The  SHG  work  was  primarily  was  supported  by  other  sources. 

Work  Accomplished 

Work  with  previous  collaborators  inside  and  outside  URIP  was  continued.  With 
U.  Gibson,  the  Brillouin  spectra  of  very  rough  ZnS  films  were  examined  and  second 
harmonic  generation  measured  on  specially  prepared  ZnS  films.  Finally,  study  of 
polymerized  Langmuir-Blodgett  films  with  J.  Swalen  of  IBM  Almaden  was 
successfully  concluded. 

In  addition,  two  visitors  had  a  large  impact  on  this  research,  presenting  some 
new  opportunites,  specifically  in  the  area  of  Langmuir-Blodgett  (L-B)  films.  These 
films  are  fabricated  one  monolayer  at  a  time  in  a  Langmuir  trough,  and  their  range 
of  applications  in  electronics  and  optics  is  increasing  dramatically.  David  Gleed,  a 
graduate  student  from  the  University  of  Exeter,  brought  a  series  of  L-B  films 
fabricated  at  Exeter  for  investigation  by  Brillouin  scattering.  Professor  Wolfgang 
Knoll  visited  for  three  months,  made  the  L-B  trough  here  operational  again,  and 
produced  films  which  were  investigated  both  by  Brillouin  scattering  and  SHG. 

The  program  has  lost  one  of  its  key  personnel,  namely  Dr.  B.  Hillebriiids,  who 
returned  to  Germany  to  finish  his  Habilitation  in  July.  Fortunately,  a  highly 
qualified  postdoctoral  replacement  was  found,  namely  John  Dutcher  from  the 
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University  of  British  Columbia,  who  was  selected  from  more  than  fifty  applicants. 
He  is  expected  to  arrive  about  October  15. 

> 

BrUlouin  Scattering  From  Rough  Surfaces 

Three  series  of  experiments  were  performed  that  probed  the  elastic  properties  of 
imperfect  surfaces. 

(  ZnS  films  were  deposited  by  the  Gibson  group  for  roughness  studies;  the 

Rayleigh  wave  penetrates  less  than  1000  A  into  the  material,  and  normally  has  its 
maximum  amplitude  right  at  the  surface,  hence  it  is  very  susceptible  to  surface 
quality.  In  fact,  it  was  found  experimentally  that  the  presence  of  a  Rayleigh  wave 
^  could  not  be  detected  by  light  scattering  in  a  sample  with  microcracks,  indicating  this 

mode  to  be  very  heavily  damped  by  the  very  rough  film  surface.  However.  Brillouin 
scattering  was  detected  from  Sezawa  modes,  which  fill  the  complete  volume  of  the 
film  and  thus  are  less  heavily  damped  by  a  rough  surface.  This  work  will  continue 
I  to  quantify  the  effects. 

The  first  step  in  understanding  how  a  "rough"  surface  affects  the  acoustic 
phonon  spectrum,  and  hence  the  mechanical  properties,  is  to  investigate  the  simplest 
possible  case  of  a  rough  surface,  namely  a  periodically  corrugated  surface.  Gratings 
I  with  periodicites  varying  from  0.2  to  0.5  ftm  were  etched  into  the  (001)  surface  of 

silicon.  Spectra  obtained  from  small-corrugation  gratings  show  evidence  for  the 
existence  of  two  Rayleigh  modes,  one  corresponding  to  that  of  a  flat  surface,  and  the 
other  frequency-renormalized  by  the  presence  of  the  grating.  The  existence  of  two 
I  modes  is  not  yet  understood.  Also  observed  was  zone  folding  at  a  wavevector 

corresponding  to  the  grating  periodicity,  as  well  as  a  longitudinal  resonance  at  the 
zone  boundary.  The  results  for  deep  gratings  show  strongly  frequency-renormalized 
modes  and  a  large  dispersion  gap  at  the  zone  boundary.  Since  all  of  the  observed 
I  features  are  not  understood,  this  work  will  continue. 

The  last  report  discussed  preliminary  Brillouin  scattering  investigations  of  the 
elastic  properties  of  saturated  ODF.  and  unpolymerized  and  polymerized  u>-ODF  and 
ODF:ODM  Langmuir-Blodgett  films.  Because  of  the  laminar  film  structure,  which 
^  consists  of  parallel  long  molecules,  aligned  perpendicularly  to  the  interface  and  tightly 

bonded  to  each  other,  the  films  are  highly  anisotropic  elastically  and  generally  show 
an  unusual  small  value  of  the  transverse  sound  velocity;  i.e.,  C44.  The  elastic 
constants  decreased  on  polymerization  for  all  of  the  films  except  DOF:ODM.  which 
I  exhibited  a  small  stiffening  of  C44,  contrary  to  widely  held  opinions  in  the  L-B 

community.  Subsequent  electron  microscope  examinations  of  the  films  showed  that 
most  of  the  films  showed  microcracks  which  occurred  during  the  polymerization 
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process.  This  correlates  with  the  large  reductions,  as  much  as  40%,  in  the  measured 
elastic  constants.  For  the  ODF:ODM  system,  there  is  negligible  film  contraction  on 
polymerization,  no  significant  microcrack  formation  occurs,  and  the  elastic  constants 
do  increase  on  polymerization.  These  results,  which  corroborate  the  ZnS 
measurements  discussed  above,  show  that  Brillouin  scattering  can  definitely  be  used  to 
detect  microcrack  formation. 

Buffer  Layer  Effects  Measured  by  Brillouin  Scattering 

Brillouin  scattering  measurements  were  performed  on  three  different  samples  of 
ZnS  films  on  Si  substrates  to  study  possible  effects  of  different  buffer  layers  between 
the  substrate  and  film  created  by  different  substrate  cleaning  procedures. 

The  three  samples  used  in  this  project  were  grown  by  physical  vapor  deposition 
after  the  Si  substrate  was  cleaned  by  three  different  methods.  Sample  1  was  cleaned 
by  the  standard  substrate-cleaning  procedure  with  deionized  water  and  a  detergent, 
followed  by  ion  milling  for  one  minute.  This  sample  served  as  a  reference,  since  it 
did  not  have  a  buffer  layer  even  though  its  surface  may  have  been  damaged  by  the 
ion  milling.  Sample  2  was  cleaned  by  the  standard  cleaning  procedure,  but  without 
the  ion  milling.  This  one  may  have  a  SiOj  buffer  layer  between  the  substrate  and 
the  film.  Sample  3  has  been  contaminated  by  pump  oil  and  therefore  probably  has  a 
thin  pump-oil  buffer  layer.  The  thicknesses  of  the  ZnS  films  were  310  A,  230  A  and 
80  A  respectively,  as  measured  with  a  WYKO  Non-Contact  Profiler. 

Tlie  phonon  dispersion  curves  were  calculated  for  both  an  isotropic  ZnS  film  on 
a  Si  substrate,  and  for  an  isotropic  2^  film  separated  from  the  Si  substrate  by  an 
isotropic  SiO]  film,  to  search  theoretically  for  any  modification  resulting  from  the 
buffer  layer.  The  calculated  curves  do  not  show  any  changes  within  the  resolution 
of  the  Brillouin  spectroscopy  technique.  This  was  confirmed  by  the  experimental 
data.  Therefore,  except  for  the  data  from  Sample  3,  which  needs  more  analysis,  no 
effect  was  observed  on  the  acoustic  phonon  spectrum  of  different  cleaning  procedures 
causing  a  possible  buffer  layer  between  the  film  and  the  substrate.  This  result 
confirmed  theory,  which  predicted  no  modification  of  the  dispersion  curves. 

Brillouin  Scattering  From  Langmuir- Blodgett  Films 

Two  additional  sets  of  experiments  were  performed  on  Langmuir-Blodgett  films 
with  the  assistance  of  the  visitors.  One  of  the  questions  addressed  was  whether  the 
elastic  properties  of  L-B  films  depend  on  the  conditions  of  the  water  surface  in  the 
Langmuir  trough.  A  number  of  films  of  22-tricosenoic  acid  were  prepared  with 
increasing  water  surface  pressures,  and  the  Brillouin  spectra  then  were  measured. 
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(22-tricosenoic  acid  is  being  developed  as  a  high-resolution  photoresist  material.)  It 
was  found  that  both  the  Rayleigh  wave  and  first  Sezawa  mode  velocities  increased 
with  increasing  surface  pressure,  indicating  a  corresponding  increase  in  the  shear 
elastic  constant.  Differences  in  the  Brillouin  spectra  were  also  observed  for  phonon 
wavevectors  parallel  (versus  perpendicular)  to  the  film  "pulling"  direction,  revealing 
in-plane  asymmetry  in  the  elastic  constants.  A  detailed  analysis  of  the  data  for  one 
surface  pressure  yielded  the  film  elastic  constants.  Although  the  values  of  the 
longitudinal  elastic  constants  Cn  and  Csj  were  about  60%  of  the  values  found 
previously  for  CdA,  the  shear  constant  c^s  was  about  2.5  times  larger.  Therefore, 
although  22-tricosenoic  acid  has  a  small  Cjs/Cu  ratio  (Css/Cjj  also)  of  approximately 
1:10,  it  is  not  as  anomolously  small  as  CdA,  which  was  between  1:20  and  1:40. 

A  second  set  of  experiments  was  carried  out  on  films  made  from  GDIO  L-B 
molecules.  From  symmetries  in  the  Brillouin  spectra,  it  was  concluded  that  the  films 
were  indeed  hexagonal  in  structure.  In  contrast  to  other  L-B  films,  the  c-axis  lies  in 
the  plane  of  the  film,  with  the  long  molecular  axis  aligned  along  the  dipping 
direction.  By  using  Brillouin  spectroscopy  to  probe  the  phonon  modes  in  the  a-c  and 
a-b  planes,  and  to  identify  spectra  attributable  to  longitudinal  guided  waves,  the 
complete  set  of  elastic  constants  was  determined  for  an  L-B  film  for  the  first  time. 
As  noted  for  22-tricosenoic  acid,  the  ratio  of  the  shear  to  the  longitudinal  elastic 
constants  was  approximately  1:10. 

Surf  ace- Enhanced  Brillouin  Scattering  From  Metal  Films 

The  Brillouin  scattered  signals  from  metal  surfaces  are  usually  very  weak,  and 
hours  typically  are  required  to  accumulate  a  spectrum  with  a  good  signal-to-noise 
ratio.  By  coupling  the  incident  light  to  surface  plasmons,  the  electromagnetic  field  at 
the  surface  can  be  enhanced,  leading  to  stronger  signals.  Further,  if  the  Brillouin 
scattered  light  is  taken  out  in  the  form  of  surface  plasmons,  the  scattering  cross 
section  is  enhanced  further. 

A  silver  film  500  A  thick  was  deposited  onto  the  base  of  a  glass  prism.  The 
incident  light  was  coupled  into  a  surface  plasmon  mode  (as  evidenced  by  a  decrease 
in  the  reflected  light).  The  light  backscattered  into  a  surface  plasmon  mode  was 
outcoupled  using  the  same  prism  and  analysed  with  a  tandem  interferometer.  When 
the  direction  of  the  incident  light  was  angle  tuned  through  the  surface  plasmon 
resonance,  the  Brillouin  signals  for  scattering  from  the  Rayleigh  wave  and  the  first 
Sezawa  mode  were  enhanced  by  a  factor  of  5o0  over  the  off-resonance  case.  The 
observed  enhancement  versus  angle  was  in  excellent  agreement  with  theory. 


Second  Harmomc  Generation  From  Thin  Films 

This  is  the  first  year  that  URIP  has  provided  some  nominal  support  for  this 
prc^ram.  In  the  past  SHC  has  been  shown  to  be  a  very  powerful  probe  of  the 
anisotropy  of  thin  films. 

The  anisotropy  of  ZnS  films  deposited  under  varying  conditions  was  probed 
using  SHG.  The  results  are  discussed  under  Task  1  report. 

With  the  aid  of  Wolfgang  Knoll,  a  new  and  very  powerful  application  of  SHG 
to  monolayer  films  deposited  by  the  Langmuir-Blodgett  technique  was  demonstrated. 
Previous  SHG  studies  have  derived  information  about  molecular  hyper-polarizabilities 
or  orientational  properties  of  single  molecules  or  parts  of  molecules.  It  was  shown 
that  structural  information  is  also  available  at  the  supra-molecular  level  of  single 
crystallites  containing  many  molecules. 

When  the  cyanine  dye  S120  is  spread  and  compressed  on  the  water  surface  of  a 
Langmuir  trough,  crystalline  aggregates  (called  J  or  Scheibe  aggregates)  are  formed  for 
pressures  above  a  certain  critical  pressure.  It  is  generally  believed  that  the 
aggregates  are  randomly  distributed  on  the  water  surface,  and  that  this  random 
distribution  is  transferred  to  the  substrate.  By  performing  SHG  from  such  a 
deposited  monolayer  and  rotating  the  sample  about  the  incident  laser  beam,  it  was 
shown  that  the  symmetry  of  the  aggregates  is  C}.  and  that  for  certain  deposition 
conditions,  the  aggregates  are  not  randomly  distributed  on  the  surface.  In  fact,  it 
was  shown  that  the  results  depend  on  the  net  molecular  flow  velocity  in  the  region 
of  the  trough  when  the  barriers  are  moved  to  compensate  for  loss  of  material  when  a 
monolayer  is  being  deposited.  When  a  tank  designed  for  zero  velocity  at  the 
deposition  point  is  used,  there  is  well-defined  net  orientation  for  the  crystallites.  On 
the  other  hand,  when  the  deposition  is  performed  under  conditions  in  which  there  is 
a  net  flow  velocity  for  the  molecules,  the  orientations  of  the  crystallites,  which 
individually  still  have  C2  symmetry,  are  random. 

Plans 

Use  of  Brillouin  scattering  to  quantify  and  interpret  the  elastic  properties  of  thin 
films  will  continue.  In  particular,  investigation  will  continue  into: 

1.  the  fascinating  elastic  properties  of  L-B  films; 

2.  the  effects  of  roughness  on  surface  mechanical  properties,  including  grating 
structures; 

3.  the  elastic  properties  of  films  fabricated  by  different  techniques. 
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Some  new  directions  are  expected  that  will  incorporate  the  ideas  and  interests  of 
the  new  postdoctoral  fellow. 

In  the  area  of  SHG  from  films,  there  are  many  interesting  possible  experiments, 
as  is  typical  during  the  evolution  of  a  new  technique.  These  will  be  pursued  as 
they  become  apparent.  Study  of  the  properties  of  thin  film  deposited  under 
different  conditions  in  order  to  optimixe  film  properties  and  study  of  L-B  films  will 
continue. 

Collaborations 

U.  J.  Gibson  (UA):  adhesion  properties  of  thin  films 

C.  M.  Falco  (UA):  preparation  of  superlattice  samples 

J.  Swalen  (IBM  San  Jose):  preparation  of  Langmuir-Blodgett  films 

J.  R.  Sambles  (University  of  Exeter):  preparation  of  Langmuir-Blodgett  films 

W.  Knoll  (Max  Planck  Mainz):  preparation  of  Langmuir-Blodgett  films 
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TASK  5 


A  CONTINUUM  MODEL  FOR  THE  MORPHOLOGY  OF  THIN  FILMS 
Principal  Investigator:  Seth  Lichter 
Project  Goal 

While  some  features  of  film  morphology  are  correlated  to  substrate  temperature 
by  the  structure  zone  model  of  Movchan  and  Demchishin.  the  ability  to  predict  film 
properties  (such  as  film  density,  rate  of  growth,  and  column  size)  from  deposition 
conditions  is  an  actively  pursued  goal  for  both  its  theoretical  interest  and  practical 
importance.  The  purpose  of  this  research  is  to  analyze  how  deposition  parameters 
affect  film  morphology.  This  is  done  by  first  identifying  the  physical  mechanisms 
which  control  adatom  motion  on  the  film.  Then,  a  mathematical  model  including 
these  effects  is  formulated. 

Work  Accomplished 

To  date,  the  work  accomplished  may  be  summarized  as  follows:  1)  the 
equations  describing  the  evolution  of  a  growing  film  have  been  formulated,  and 
2)  the  behavior  of  the  equations  in  some  limiting  cases  has  been  determined.  A  more 
extensive  description  of  how  the  film  responds  under  a  variety  of  deposition 
conditions  is  forthcoming,  pending  completion  of  the  numerical  coding  necessary  to 
determine  some  of  the  parameters  which  appear  in  the  equations.  These  results  will 
be  available  in  a  preprint  later  this  month. 

The  large-scale  properties  of  a  thin  film  are  obtained  by  averaging  over  an 
ensemble  of  individual  particles  passing  through  a  control  volume  on  the  surface. 
The  local  analysis  is  performed  in  a  coordinate-independent  fashion.  The  local 
balance  then  yields  deposition  equations  in  general  coordinates  that  have  the 
advantage  of  being  fully  nonlinear  for  the  macroscopic  scales.  A  simple  ballistic 
argument  yields  a  scattering  cross  section  that  determines  the  number  of  adatoms 
added  to  the  surface  layer  and  their  velocity  along  the  surface. 

The  basic  assumption  is  that  the  incident  adatoms  hit  the  substrate  and  a 
portion  of  them  join  a  loosely  bound  surface  layer.  The  adatoms  in  this  surface 
layer  then  eventually  stick,  leading  to  accretion  of  the  deposited  film. 
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The  formulation  incorporates  the  effects  of  incoming  momentum,  adatom  motion, 
including  convection  and  diffusion,  and  sticking.  The  incoming  beam  is  assumed  to 
be  unidirectional  and  uniform.  The  flux  rate  onto  a  particular  portion  of  the  surface 
is  then  dependent  upon  the  orientation  of  the  surface.  Incident  particles  will  hit  the 
surface  and  many  will  re-evaporate.  A  simple  ballistic  argument  determines  the 
portion  of  the  beam  flux  which  is  not  re-evaporated.  Adatoms  move  by  convection 
and  surface  diffusion  until  they  stick.  Bulk  diffusion  is  ignored.  Finally,  it  is 
assumed  that  adherence  to  the  surface  is  a  Poisson  process  and  that,  consequently,  the 
deposition  rate  is  proportional  to  the  adatom  concentration. 

Consider  a  control  volume  on  the  surface  whose  location  is  given  by  p.  In  the 
control  volume,  the  concentration  of  adatoms,  c,  will  be  augmented  by  a  portion  of 
the  incoming  flux  through  the  top  control  surface.  On  the  lower  control  surface,  a 
certain  proportion  of  the  adatoms,  a,  stick  per  unit  time.  Adatoms  also  leave  through 
the  sides  of  the  control  volume  because  of  their  non-zero  velocity. 

A  particle  that  is  deflected  along  the  surface  will  travel  with  a  fixed  velocity, 
V.  until  it  strikes  some  other  particle.  It  is  assumed  that  the  probability  of  the 
particle  being  deflected  a  second  time  is  a  Poisson  process  with  time  scale  I/7.  The 
average  length  of  a  deflection  is  then  given  by  V/7.  It  is  assumed  that  both  of  these 
scales  are  small  compared  to  the  characteristic  scales  of  the  surface;  that  is. 


The  density  of  coherent  adatoms,  c,  can  now  be  written  in  terms  of  the  ballistic 
cross  section  by  integrating  the  contributions  from  the  incident  flux  in  the 
neighborhood  of  a  particular  point,  p.  on  the  surface, 

f2ff  r«» 

c(p,  t)  -  J  d^fr  J  dr  J^(p  -  r^.  t  -  t’)e-'if‘'  .  (2) 

where  the  local  polar  coordinates  (r,  iji)  have  been  introduced  on  the  surface.  The 
flux  has  to  be  normalized  by  the  velocity  V  to  give  a  particle  density  per  unit 
length.  The  exponential  term  is  the  probability  density  that  a  particle  remains 
undeflected  after  a  time  t',  where  t'  is  the  transit  time  from  the  location  where  the 
particle  strikes  p  -  r\i  to  the  location  p,  or 
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The  integral  Eq.  (2)  can  be  simplified: 


» 


i 


C(P.  t)  -  ^  Jn  -  (J«>t  - 


(4) 


Two  local  properties  of  the  ballistic  cross  section  are  defined  here;  first,  the  total 
incident  flux, 


J 


a 


(5) 


and  second,  the  average  velocity  of  an  incident  particle: 


The  subscript  s  on  the  vector  operators  indicates  that  they  are  acting  along  the 
surface  defined  by  p.  Using  the  condition  in  Eq.  (1),  it  is  clear  that  the  leading- 
order  solution  to  the  concentration  is 

c(p,  t)  -  i  J„  ,  (7) 

which  can  be  interpreted  physically  as  the  incident  flux  multiplied  by  the  average 
time  to  deflection  of  an  adatom.  The  first-order  corrections  to  Eq.  (7)  can  now  be 
written  in  terms  of  c  instead  of  J„: 


where  the  last  two  terms  in  Eq.  (8)  are  the  leading-order  corrections  for  the  temporal 
fluctuation  and  the  mean  advection  of  adatom  concentration. 

Equation  (8)  describes  the  evolution  of  c  attributable  to  the  coherent  motion  of 
particles  caused  by  ballistic  scattering.  In  addition  to  this  motion,  a  contribution 
from  the  motion  of  the  particles  caused  by  random  scattering  on  the  surface  can  be 
included.  It  is  assumed  that  if  a  particle  in  the  surface  layer  is  deflected,  the 
orientation  of  its  velocity  will  be  determined  randomly.  Consequently,  the  motion  of 
the  particle  can  be  modeled  as  a  diffusive  process. 
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The  rate  at  which  particles  leave  the  coherent  adatom  concentration  is  given  by 
7C.  It  is  assumed  that  adatoms  become  part  of  an  incoherent  adatom  density,  d. 
Adatoms  are  removed  from  the  incoherent  state  by  sticking  to  the  surface;  this  is 
again  assumed  to  be  a  Poisson  process  with  a  mean  lifetime  of  q/7'.  The  equation 
of  motion  for  d  is  then  given  by 

dt  -  7C  -  7'd  +  tfV^d  .  (9) 

where  </  is  a  diffusion  constant  which  can  be  thought  of  as  the  product  of  the  mean 
speed  and  the  mean  free  path. 

As  adatoms  are  deposited,  the  surface  grows  outward  along  its  normal  in 
proportion  to  the  rate  of  deposition;  that  is. 

Pt  -  7'dn  .  (10) 

This  equation,  together  with  Eqs.  (8)  and  (9),  specifies  the  evolution  of  the 
surface.  To  summarize,  the  equations  governing  the  growth  of  the  interface  are 

Ct  -  J„  -  7C  -  V,  .  I^V^c]  , 

d|  -  7c  -  7'd  +  t'V^d  , 

Pt  -  7  dn  , 

where 

p  -  location  of  the  surface, 

n  -  surface  normal, 

c  -  coherent  adatom  volume  density, 
d  -  incoherent  adatom  volume  density, 

-  adsorbed  adatom  flux, 

Vg  -  average  adsorbed  adatom  velocity, 

1/7  -  coherent  adatom  deflection  rate, 

1/7'-  incoherent  adatom  sticking  rate, 

1/  m  incoherent  adatom  diffusion. 

Results  from  the  above  formulation  will  be  contained  in  the  attached  preprint  of 
an  article  to  be  submitted  to  Physical  Review  B. 

Checking  some  of  the  more  refined  results  of  the  model  with  experiments  being 
performed  by  J.  Ruffner  and  U.  Gibson  in  the  Optical  Sciences  Center  w£is  intended. 
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but  there  was  not  sufficient  time.  However,  in  the  preprint,  the  results  are 
compared  with  previous  experimental  data. 


Personnel 
S.  Lichter 

Collaboration 

A.  J.  Bemoff,  Department  of  Mathematics,  University  of  Arizona 
U.  J.  Gibson:  Experimental  confirmation  of  model  results 
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TASK  6 


STUDIES  OF  SURFACE  ROUGHNESS  HEALING 
Principal  Investigator:  C.  M.  Falco 
Project  Goal 

The  goal  of  this  task  is  to  develop  t^hniques  for  determining,  and  reducing,  the 
roughness  of  substrate  surfaces.  Scanning  tunneling  microscopy  (STM)  •'ind  optical 
interferometry  techniques  will  be  developed  to  study  the  roughness  of  surfaces  on 
length  scales  ranging  from  angstroms  to  millimeters.  In  combination  with  these 
studies,  sputtering  and  molecular  beam  epitaxy  (MBE)  techniques  will  be  developed  to 
deposit  "buffer  layers"  to  reduce  the  roughness. 

Summary 

The  roughness  at  the  interfaces  within  multilayer  coatings  has  a  negative  effect 
on  most  of  the  desirable  physical  properties  of  such  coatings.  Multilayers  currently 
being  produced  for  some  purposes  have  layer  thicknesses  on  the  order  of  1  nm  or 
less.  In  such  structures,  even  sub-nanometer  roughness  can  have  a  large  effect.  As 
a  result,  it  has  become  a  fairly  common  practice  to  pre-deposit  an  initial  thick  layer 
of  some  material  on  a  substrate  in  hope  of  healing  the  roughness,  before  beginning 
deposition  of  the  actual  thin  film.  However,  systematic  studies  of  the  effects  of  these 
"buffer  layers"  have  yet  to  be  conducted.  Of  great  interest  is  understanding  of  the 
precise  mechanism  (presumably  involving  preferential  surface  migration  and  nucleation 
in  the  microscopic  "valleys")  that  results  in  healing  of  the  roughness.  Study  of  the 
effect  of  such  buffer  layers  on  surface  roughness,  as  well  as  the  effect  of  other 
substrate  pre-treatments,  has  begun. 

Work  Accomplished 

Substrate  Preparation 

Si  films  have  been  deposited  by  sputtering  on  a  variety  of  substrates,  and  the 
effect  of  the  films  on  surface  roughness  characterized  with  the  techniques  described 
below.  The  substrates  studied  to  date  include  diamond-turned  aluminum,  float  glass. 
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and  polished  semiconductor-grade  Si  wafers.  These  substrate  materials  are  important 
for  future  thin-film  work  in  the  visible  and  XUV  optics,  electronics,  and  magnetism 
fields.  Si  wafers  and  float  glass  were  treated  with  three  commonly  used  cleaning 
procedures  to  determine  their  effects  on  roughness.  Simple  detergent  cleaning,  a 
three-step  acid  cleaning  procedure  and  a  22-step  acid  cleaning  procedure  were  used 
to  prepare  Si  wafers  for  growth  of  films  by  MBE. 

Characterization  of  Surface  Roughness 

Substrate  surfaces  were  characterized  with  a  scanning  tunneling  microscope 
(STM)  and  a  WYKO  phase-shifting  interferometer,  both  before  and  after  the 
deposition  of  Si  films  by  sputtering.  The  STM  used  for  this  initial  work  was  made 
available  through  a  collaboration  with  Dr.  Dereniak's  group. 

The  STM  and  the  WYKO  interferometer  have  comparable  vertical  resolutions: 
1  A  for  the  former  instrument  and  3  A  for  the  latter.  However,  their  horizontal 
resolutions  differ  by  a  factor  of  10*;  that  is.  approximately  1  A  for  the  STM  and  2 
/tm  for  the  WYKO.  Thus,  the  WYKO  measurement  determines  the  "waviness"  of  the 
surface,  while  the  STM,  with  its  atomic  scale  resolution,  measures  roughness  over 
spatial  dimensions  that  have  large  effects  on  the  physical  properties  of  many  thin 
films  and  multilayers. 

The  STM  measurements  made  to  date  were  complicated  slightly  by  the  fact  that 
the  instrument  used  operates  in  air.  and  by  the  instrument  not  always  being 
immediately  available  for  use.  Since  this  technique  requires  a  conductive  surface,  a 
thin  layer  of  Ag  (100  A)  was  deposited  over  the  Si  buffer  layer  to  prevent  the 
growth  of  an  insulating  oxide  layer  prior  to  the  STM  measurement.  This  also  was 
necessary  for  the  float  glass,  in  order  to  create  a  conductive  surface. 

STM  surface  profiles  of  a  diamond-turned  aluminum  flat  were  found  to  show  a 
small  periodic  ripple  in  the  surface  before  deposition,  presumably  an  artifact  of  the 
diamond-turning  process.  However,  measurements  of  the  auto-correlation  function  of 
the  surface  after  deposition  of  the  buffer  layer  showed  the  surface  to  have  become 
much  more  random.  Analyzing  many  such  measurements  leads  to  the  conclusion  that 
the  Si  buffer  layer  succeeds  in  randomizing  the  diamond-turned  surface,  but  has  little 
effect  on  the  root-mean-square  (RMS)  roughness.  Thus,  in  the  case  of  diamond- 
turned  aluminum,  there  may  be  some  advantage  to  depositing  a  Si  buffer  layer,  since 
a  random  surface  can  have  a  more  predictable  effect  on  the  properties  of  a  coating 
than  one  with  periodic  fluctuations.  For  the  other  two  substrate  materials  studied, 
polished  Si  and  float  glass,  the  deposition  of  an  Si  buffer  layer  did  not  improve 
either  the  RMS  roughness  or  the  randomness  of  the  surface.  However,  this 
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prelioainary  conclusion  should  be  re-examined  after  further  measurements  are  made 
on  Si  buffer  layers  deposited  under  other  conditions. 

WYKO  measurements  also  were  performed  on  the  substrates  before  and  after 
the  various  cleaning  procedures.  The  three-step  acid  etch  was  found  to  leave  a  clean 
but  otherwise  unmodified  surface.  The  detergent  method  left  rough  patches,  while 
the  22-step  acid  procedure  increased  the  roughness  of  the  polished  wafers  over  the 
entire  area  sampled. 

Plans 

One  major  goal  of  this  work  is  to  develop  techniques  to  determine  the  auto¬ 
correlation  of  surface  roughness  over  length  scales  continuously  spanning  the  range  of 
approximately  1  A  to  1  mm.  Initial  results  toward  this  goal  have  been  highly 
encouraging.  The  continuation  of  these  roughness-healing  studies  will  make  use  of 
improved  STM  hardware  and  software  to  bridge  the  gap  in  horizontal  resolution 
between  the  STM  and  the  WYKO  measurements.  During  the  coming  year,  a 
Nanoscope  II  STM  will  be  used  to  make  measurements  over  much  larger  areas  than 
was  possible  previously,  and  then  averaging  software  will  be  used  to  produce  an 
accurate  picUire  of  the  roughness  over  a  variety  of  intermediate  horizontal-length 
scales. 

A  systematic  study  of  roughness  as  a  function  of  surface  pretreatment  and 
sputter-deposition  parameters  will  be  performed.  A  search  will  be  made  for  an 
improved  MBE  substrate-cleaning  procedure  that  can  remove  chemically  bonded 
impurities,  such  as  carbon,  without  increasing  the  roughness  of  the  polished  surface. 

Personnel 

C.  M.  Falco 
Z.  Milanovic 
J.  M.  Slaughter 

Collaborations 

Dereniak:  STM  studies  of  substrates 
Gibson:  STM  studies  of  surfaces 
Leavitt:  RBS  analysis  of  films 
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TASK  7 


CHARACTERIZATION  FACILITY 
Principal  Investigator:  U.  J.  Gibson 
Project  Goal 

Set  up  a  central  facility  for  characterization  of  thin  films,  surfaces,  and 
substrates. 

Work  Accomplished 

The  Characterization  Facility  became  fully  operational  in  the  last  year,  with 
construction  of  a  dedicated  site,  installation  of  the  electron  microscopes,  acquisition  of 
sample  preparation  and  other  equipment,  and  mastery  of  their  operation.  The  facility 
is  seeing  increased  use  from  a  variety  of  the  URl  participants. 

Equipment 

TEM:  JEOL  2000FX,  with  1.2  A  resolution  and  micro-diffraction  capability. 

SEM:  JEOL  JSM-U3;  should  be  replaced  with  new  instrument  by  mid-December. 
WYKO  interferometric  profilometer 

OLYMPUS  analytical  microscope  with  phase  contrast  and  TV  monitor 
NANOSCOPE  I  scanning  tunneling  microscope  (STM) 

Sample  preparation  equipment  for  SEM  and  TEM 
Scatterometer/scanning  microscope  from  Task  8 

Experiments  in  Progress 

Use  of  crystalline  Si,  SiOj,  and  opal  for  polishing  damage  studies 
STM  imaging  of  polished  and  coated  substrates 
Development  of  microtome  techniques  for  thin-film  samples 


! 
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PUns 

Acquire  and  install  scanning  coils  and  a  TV  for  the  TEM 

I 

Acquire  and  install  a  cold  stage  for  the  ion  thinner 
Upgrade  the  STM  for  digital  data  acquisition 

Personnel 

I 

U.  J.  Gibson 

V.  Lindley 
M.  Kelley 
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TASK  8 


MICRO-SCANNING  SCATTEROMETER  AND  POLISHING  STUDIES 
Principal  Investigator:  R.  Parks 
Project  Goal 

Design  and  implemention  of  a  surface-scattering  measurement  instrument,  with 
good  lateral  resolution  and  a  greater  sensitivity  to  surface  rather  than  volume  defects, 
was  the  major  goal  of  this  Task.  Comparative  studies  of  polishing  procedures  were 
also  performed. 

Summary 

An  instrument  has  been  constructed  that  measures  scattered  light  through  a 
small  aperture  as  a  sample  is  scanned.  The  instrument  and  its  software  are 
functional,  and  the'  instrument  was  evaluated  by  W.  Wolfe,  whose  assessment  points 
out  that  the  instrument  is  too  unlike  conventional  scatterometers  to  be  calibrated 
against  them.  It  is,  however,  more  sensitive  to  local  variations  in  scattering 
efficiency  than  other  techniques  such  as  BRDF. 

In  tests  against  a  WYKO  interferometer,  the  instrument  appeared  less  sensitive 
to  surface-height  variations  than  the  WYKO  and  lateral  resolution  was  comparable. 

It  is  not  yet  clear  that  the  instrument  has  merits  above  existing  techniques  for 
characterizing  bare  surfaces.  A  study  to  determine  if  it  can  be  used  to  characterize 
thin  films  by  means  of  light  scattered  out  of  a  waveguide  may  be  undertaken. 

Work  Accomplished 
Scatterometer 

During  the  past  year,  the  micro-scanning  scatterometer  was  completed  as  a 
stand-alone  device,  and  software  for  its  control  and  image  processing  was  written  and 
documented.  Following  a  period  of  check  out  and  calibration,  the  device  was  turned 
over  to  W.  Wolfe's  group  for  calibration  against  sample  of  known  BRDF 
characteristics;  his  report  appears  as  an  appendix  later  in  this  document. 

Currently,  the  instrument  is  fitted  with  a  sample  holder  that  accepts  I -in. 
diameter  by  0.12S-in.  samples.  A  calibration  sample  is  used  to  center  the  scan  and 
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set  the  laser  Uluminator  and  photometer  to  a  constant  level  before  scanning  actual 
samples.  Samples  may  be  scanned  over  areas  1  cm‘  or  smaller,  in  user-defined  steps 
down  to  1  fim.  Sample  resolution  at  the  detector  is  currently  7  fim.  This  could  be 
reduced  by  a  factor  of  three  before  the  illumination  level  would  be  too  weak  to  be 
practical. 

During  the  check  out  and  calibration  phase,  a  10-/tm  slit  etched  in  chrome  was 
used  to  determine  the  instrument  profile.  Both  sides  of  the  slit  were  clearly  visible 
to  the  scatterometer  and  there  was  no  evidence  of  stray  light  giving  a  false  or  ghost 
image  of  a  second  slit.  The  instrument  was  alw  checked  for  fixed  pattern  noise 
caused  by  speckle  from  the  laser  illumination  system.  Again,  there  was  no  evidence 
of  fixed  pattern  noise. 

Measurement  of  a  Lambertian  surface  (as  determined  by  BRDF  with  a  2-cm 
spot  size),  revealed  intensity  variations  point-to-point  with  a  S:1  ratio.  This  is  the 
result  of  the  small  aperture  and  large  collection  angle  of  the  new  scatterometer, 
which  meastires  short  distance  variations  in  the  surface. 

Polishing  Studies 

Several  surfaces  with  different  degrees  of  roughness  and  damage  were  produced 
for  the  Thin  Film  Center  and  to  support  development  of  the  scatterometer. 

A  polishing  procedure  was  established  that  generates  a  surface  on  a  f used-silica 
substrate  with  a  roughness  of  about  0.2-nm  rms.  (This  is  the  measurement  limit  of 
the  WYKO  profilometer.)  Crystals  of  silicon  and  quartz  were  polished;  a  surface 
roughness  of  0.8-nm  rms  was  produced  on  silicon  substrates,  and  surfaces  on  quartz 
substrates  with  a  roughness  of  about  0.3-nm  rms  were  produced. 

A  comparison  of  Silvemail  polishing  compounds  is  underway  using  fused-silica 
substrates.  Preliminary  results  show  that  these  compounds  can  produce  very  smooth 
surfaces,  but  the  results  are  limited  by  the  measurement  capability  of  the 
profilometer. 

Some  experiments  were  conducted  with  electroless  nickel  to  obtain  metal 
substrates  for  further  polishing  studies. 

Personnel 

R.  E.  Parks 
R.  Sumner 
J.  Sasian 
C.  Kittrell 
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Development  Proposal 
SPECTRAL  ELLIPSOMETRY 
Principal  Investigator:  R.  Zanoni 
Project  Goal 

Funding  is  requested  for  development  of  a  research-grade  spectroscopic 
ellipsometer.  This  instrument  will  be  based  on  a  new  design  that  combines  multiple 
angle-of-incidence  measurements  with  a  spectroscopic  phase-modulation  ellipsometer 
(MA-SPME).  The  accuracy  of  commercial-grade  ellipsometers  is  at  least  an  order  of 
magnitude  smaller  than  that  of  a  PME  operating  at  one  wavelength.  With  the 
multiple  wavelength  and  multiple  angle-of-incidence  capabilities,  this  new  instrument 
will  have  increased  accuracy  and  additional  measurement  capabilities.  The  MA- 
SPME  will  be  able  to  measure  the  refractive  index  of  optical  surfaces  and  films,  film 
thickness,  inhomogeneous  films,  spectral  absorption,  birefringence,  crystal  symmetry, 
and  surface  roughness. 

Introduction 

The  research  programs  funded  by  URIP  require  basic  characterization  of  the 
linear  optical  properties  of  thin  films  such  as  refractive  index,  absorption,  dispersion, 
thickness,  and  birefringence  over  a  wide  range  of  optical  wavelengths.  These 
programs  include: 

*  Brillouin  spectroscopy 

*  Thin-film  growth  and  characterization 

*  Nonlinear  optics  of  thin  films: 

-  Organic  film  characterization 

-  Semiconductor  characterization 

-  Second  harmonic  generation 

*  Surface  studies: 

-  Roughness 

-  Optical  constants 

-  Effects  of  surface  polishing 
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For  example,  Brillouia  spectroscopy  requires  information  about  thin-film  optical 
symmetries,  refractive  index,  and  thickness  for  films  as  thin  as  several  Langmuir- 
Blodgett  monolayers  (1  monolayer  at  3  nm).  Second  harmonic  generation  requires  the 
measurement  of  index  of  refraction  and  thickness  for  films  over  a  wavelength  range 
of  O.S  to  2.0  fim.  Since  investigation  of  new  materials  is  becoming  an  increasingly 
important  area  in  optics,  the  ability  to  characterize  these  new  materials  quickly  and 
'  accurately  is  crucial  in  order  to  utilize  these  materials  within  existing  optics 

programs.  A  variety  of  methods  are  used  currently  to  make  these  measurements,  but 
a  single  work  station  that  can  routinely  characterize  new  thin-film  materials  would 
be  an  invaluable  addition  to  the  Center. 

EUipsometry 

Invented  nearly  100  years  ago,  ellipsometry  provides  a  convenient  method  for 
,  obtaining  optical  properties  by  measuring  the  intensity  of  polarized  light  reflected 

from  an  optical  surface.  The  interpretation  of  ellipsometric  data  is  not  a  simple  task, 
and  thus  the  development  of  ellipsometry  has  paralleled  the  development  of 
computers.  Within  the  past  25  years,  new  variations  of  ellipsometry  have  been 
implemented  and  tested.  It  is  now  possible,  with  the  appropriate  instrumentation,  to 
make  measurements  on  bulk  or  thin-film  samples  that  may  be  anisotropic, 
inhomogeneous,  microscopically  rough,  covered  with  adsorbed  layers,  or  have  finite 
interface  widths. 

I  Conventional  null  ellipsometry  (CNE)  used  nulling  methods  to  find  the 

polarization  states  that,  on  reflection  from  a  sample,  extinguish  the  light  incident  on  a 
detector.  The  quadratic  dependence  of  the  intensity  on  instrument  settings  near  null, 
and  the  necessity  of  operating  the  detector  at  very  low  light  levels,  limit  null 
I  instruments  to  intense  spectral  line  sources.  Optical  characterization  requires  accurate 

spectral  data,  and  this  is  the  domain  of  spectral  ellipsometry  (SE).  Spectral 
ellipsometry  is  dominated  by  photometric  instruments  in  which  one  of  the  elements  is 
rotated  or  modulated  to  vary  systematically  the  flux  reaching  the  detector.  In 
»  addition  to  being  easily  automated,  these  instruments  operate  at  high  optical  efficiency 

and  are  well  suited  for  use  with  weak  continuum  sources. 

There  are  three  classes  of  spectral  ellipsometer;  rotating  analyzer  (RAE), 
rotating  polarizer  (RPE),  and  phase  modulated  (PME).  The  RAE  and  RPE  have  the 
►  advantage  of  simplicity,  since  the  only  wavelength-dependent  element  is  the  sample 

itself.  Systematic  errors  are  less  likely  with  the  RAE  than  with  the  RPE  because 
polarization-insensitive  detectors  are  available,  while  polarization-free  sources  are  not. 
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The  PME  is  not  affected  by  source-polarization  or  detector-polarization  sensitivity,  is 
capable  of  high-speed  operation,  and  has  no  mechanically  moving  parts.  So.  the 
spectroscopic  PME  is  the  instrument  of  choice  for  characterizing  thin-film  materials. 


EUipsometric  Accuracy 

EUipsometers  measure  the  complex  reflectivity  p  -  Rp/R,  -  tan^e*^^  of  a 
reflecting  sample,  where  Rp  and  R,  are  the  p-  and  s-polarized  light  reflectivities. 
The  ellipsometric  angles,  and  A,  are  related  to  the  polarizer  angle  and  the  analyzer 
angle  in  CNE.  These  ellipsometric  angles,  together  with  the  angle  of  incidence, 
determine  the  ultimate  accuracy  for  film  thickness  and  refractive  index  in  a  given 
ellipsometric  measurement.  The  accuracy  of  a  thickness  measurement  made  by 

ellipsometry  as  a  function  of  wavelength  is  shown  in  Fig.  1.  The  accuracy  of  the 
instrument  increases  to  1%  for  short  wavelengths  (as  250  nm).  provided  that  the 
ellipsometric  angles  0,  A,  and  0  are  accurate  to  ±  0.001°.  The  accuracies  of 

thickness  and  refractive  index  measurements  at  \  -  0.6328  pm  for  a  commercial- 
grade  instrument  versus  the  proposed  research  instrument  are  listed  in  Table  1.  (The 
accuracies  of  the  ellipsometric  angles  for  the  commercial  instrument  are:  0  -  ±  0.04°, 
A  -  ±  0.08°,  and  0  -  ±  1/60°.  These  values  were  obtained  from  the  literature  for  a 
commercial  instrument  made  by  SOPRA  (ESG-2).) 


Table  1. 

Measurement  accuracy 
ellipsometers. 

at  X  -  0.6328  pm 

for  research 

and  commercial 

Film 

Thickness 

(nm) 

Research-grade 
instrument  (±  0.001°) 

%(5tr™,/t)  %(5n™,/n) 

Commerciiil-grade 
instrument  (SOPRA) 
%(5tr„,s/t)  %(«n,„,/n) 

2.5 

17.0 

0.89 

460.0 

24.0 

5.0 

4.3 

0.29 

120.0 

6.0 

10.0 

1.1 

0.20 

30.0 

1.5 

50.0 

0.72 

0.19 

2.1 

0.08 

100.0 

0.015 

0.19 

0.75 

0.08 

The  film  model  assumed  in  this  case  is  a  polymer  film  (n  -  1.55)  on  a  fused- 
silica  substrate  (n  -  1.46).  The  research-grade  instrument  is  clearly  at  least  an  order 
of  magnitude  better  than  the  commercial  ellipsometer.  Further,  the  errors  shown  can 
be  reduced  an  order  of  magnitude  by  decreasing  the  wavelength,  as  shown  in  Fig.  1. 
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Figure  1.  Ellipsometric  accuracy  for  an  ellipsometer  with  60  =  5A  =  =  ±  0.001°: 

(a)  percent  error  in  thickness  (ht/t)  versus  wavelength;  (b)  percent  error  in  refractive 
index  (hnfn)  versus  wavelength. 


Program  Description 

Construction  of  a  research-grade,  automated  spectroscopic  ellipsometer  is 
proposed.  The  instrument  design  will  be  based  on  a  phase-modulated  ellipsometer. 
However,  unlike  a  commercial  instrument,  the  polarizer  position  and  analyzer  angle 
will  be  fixed.  The  incident  polarization  will  be  modulated  by  an  acoustic  shear 
wave  transducer.  The  accuracy  of  the  ellipsometric  angles  can  be  expected  to  be 
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±  O-OOl".  The  instrument  will  be  mounted  on  a  computer-controlled  goniometer  to 
allow  multiple  angle-of-incidence  ellipsometric  measurements.  This  will  improve  the 
overall  accuracy  of  the  ellipsometric  measurements,  as  well  as  providing  additional 
information  to  obtain  birefringence  data.  A  tungsten-halogen  light  source  will  be 
used,  together  with  a  monochromator,  to  scan  a  range  of  wavelengths  from  200  to 
800  nm.  Provision  for  IR  wavelengths  from  800  to  2000  nm  will  be  added  later. 

The  program  will  procceed  in  two  phases: 

1)  In  the  first  year,  the  basic  ellipsometer,  without  spectral  capability,  will  be 
made  operational.  A  multi-line,  HeNe  source  will  be  used  to  simplify  the 
development  of  the  ellipsometer;  in  a  parallel  effort,  the  spectral  source  and 
monochromator  will  be  assembled. 

2)  The  goniometer  and  spectral  capability  will  be  added  in  the  second  year, 
completing  the  MA-SPME. 

Development  of  this  instrument  will  allow  single-station  characterization  of  thin 
films  at  the  Center,  as  well  as  providing  insight  into  the  ultimate  accuracy  of  MA- 
SPME  instruments. 
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XPS  and  RHEED  studies  of  Si(lOO)  Substrate 
Preparation  for  Epitaxial  2LnS  Growth 
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Ursula  J.  Gibson 

Optical  Sciences  Center 
University  of  Arizona 
Tucson.  AZ  85721 


Abstract 

We  report  on  angle  resolved  photoemission  of  Si  (100)  substrates  prepared  for 
growth  of  ZnS  by  atomic  layer  or  molecular  beam  epitaxy.  In  particular,  we 
investigate  the  growth  habit  of  the  first  monolayer  of  zinc,  and  find  qualitative 
differences  in  the  behavior  as  a  function  of  substrate  cleaning  procedure,  and 
pretreatment  with  HjS  gas.  RHEED  studies  did  not  indicate  a  significant  lowering  of 
the  temperature  required  to  remove  SiOj  from  the  Si  wafer. 


53 


INTRODUCTION 


A  great  deal  of  attention  has  been  paid  in  recent  years  to  the  deposition  of 
epitaxial  thin  films  of  zinc  chalcogenides.  These  compounds  have  already  found 
significant  application  as  electroluminescent  displays,  and  their  large  direct  bangap 
has  led  to  their  pursuit  for  short  wavelength  light  emitting  diodes  and  lasers. 

Several  reports  have  been  made  of  epitaxial  growth  of  ZnS  on  Si(lOO),  by 
MBE,^  from  compound  source  materials.  By  contrast,  in  ALE.  once  growth  of  ZnS 
is  established,  sequential  pulses  of  Zn  and  S,  or  of  compounds  containing  these 
elements,  are  brought  to  the  substrate  with  a  time  lag  between  them.^  This  allows  a 
layer  by  layer  growth  of  the  material,  and  has  been  reported  to  assist  in  the 
exclusion  of  some  impurities.*  The  formation  of  the  first  layer  of  ZnS  is  obvioxisly 
critical  for  epitaxial  growth,  whether  from  an  elemental  or  compound  source  material. 
A  recent  report  on  MBE  growth  of  ZnS  on  Si  suggested  that  an  HjS  treatment  of  the 
wafer,  prior  to  deposition  of  ZnS,  enhanced  the  resulting  film  quality.*  This  work 
was  initiated  to  clarify  the  effect  of  HjS  exposure  on  a  Si  wafer,  and  to  determine 
its  effect  on  the  growth  habit  of  initial  layers  of  Zn  deposited  onto  the  surface. 


EXPERIMENTAL 

Silicon  (100)  substrates  were  cleaned  using  a  Pirhana/RCA-I/RCA- 
2/H2O/HF/H2O  procedure.*  The  final  H2O  rinse  lasted  3  or  S  minutes,  allowing 
formation  of  an  oxide  on  the  silicon  surface.  The  wafer  pieces  were  then  mounted 
with  gallium  onto  stainless  steel  one  inch  Perkin-Elmer  (P^  XPS  sample  stubs,  and 
transferred  to  the  growth  chamber  via  a  PE  sample  introduction  rod. 

The  H2S  exposures  and  Zn  depositions  were  performed  in  a  UHV  chamber 
(shown  in  Fig.  1)  pumped  by  a  combination  of  a  turbopump  and  Ti  sublimation 
pump.  A  liquid  nitrogen  cooled  source  shroud  contained  the  Zn  effusion  cell,  and 
the  H2S  was  bled  in  through  a  needle  valve,  and  a  nozzle  that  directed  the  gas 
stream  toward  the  substrate.  Base  pressures  were  typically  5x10"*“  Torr  prior  to 
heating  of  the  stub  and  introduction  of  the  H2S  gas.  The  samples  were  heated  by 
electron  bombardment  to  a  maximum  temperature  of  700-770°C.  Sample  temperatures 
were  measured  approximately  using  a  thermocouple  mounted  on  the  manipulator  stage 
and  calibrated  against  the  stub  temperature.  In  each  experiment,  the  substrate  was 
heated  until  a  RHEED  pattern  was  observed. 

Experiments  were  performed  to  determine  whether  the  RHEED  pattern  was 
observable  at  lower  temperatures  in  the '  presence  of  HjS,  and  to  measure  surface 
impurites  after  cleaning  and  heating  with  and  without  H2S.  and  to  determine  the 
growth  habit  of  thin  Zn  layers  deposited  on  substrates  prepared  these  two  ways. 

For  the  RHEED  pattern  onset  studies,  the  sample  was  loaded  onto  the 
manipulator,  and  the  stub  temperature  was  raised  in  5(PC  increments.  At  each 
temperature,  the  manipulator  was  adjusted  in  search  of  a  pattern.  This  procedure 
was  followed  both  with  and  without  a  backpressure  of  1x10"^  Torr  of  HjS  gas.  The 
procedure  for  the  depositions  was  as  follows:  After  cleaning  and  loading,  the 
samples  were  heated  to  700®C  for  approximately  30  minutes.  The  samples  were  then 
cooled  to  300®C,  which  required  approximately  30  minutes.  If  HjS  was  used,  it  was 
admitted  to  a  pressure  of  1x10“’’  Torr  for  30  minutes,  with  the  substrate  held  at 
300®C.  These  substrates  will  be  referred  to  as  Si(A).  If  no  H2S  was  admitted  (Si(B) 
sample),  the  substrate  was  held  at  300®C  to  keep  heating  times  equal.  The  samples 
were  then  either  transferred  in  an  unpumped  vacuum  vessel  to  the  PHI  5100  for 
XPS  analysis,  or  after  pumping  out  the  HjS,  deposited  upon  with  Zn  at  the  rate  of 
1  A  per  second,  to  a  thickness  of  100  A  indicated  on  the  quartz  crystal  monitor. 


The  difference  in  temperature  between  the  crystal  (22**)  and  the  substrate  (SOO^Q  led 
to  an  expected  deposit  thickness  of  10-20  A. 


RESULTS  &  DISCUSSION 

The  onset  temperature  for  the  RHEED  patterns  was  not  measurably  affected  by 
the  presence  of  during  heating,  however,  both  the  quantity  and  bonding  of 
carbon  and  oxygen  impurities  showed  the  influence  of  the  HjS  treatment.  The 
strongest  effect  was  seen  in  the  growth  habit  of  the  Zn  layers  after  treatment,  where 
H^S  was  seen  to  promote  layer-like  growth.  The  presence  of  an  interfacial  layer  of 
S  is  believed  to  be  responsible  for  this  behavior. 

The  angle  resolved  data  were  analyzed  using  the  model  of  C.  S.  Fadley^  for  the 
formation  of  thin  layers  on  a  substrate.  In  his  model,  the  ratio  of  the  film  to 
substrate  signals  is  described  by 

_  1^^  1  -  exp(-t/A(E‘)sinfl) 

(1  -  7)  +  yexpC-t/AIEi'lsind) 

where 

R  -  ratio  of  film  to  substrate  signal, 

K  -  geometrical  factor. 

7  -  percentage  of  the  covered  substrate. 

t  -  film  thickness, 

ACE*)  -  mean  free  path  of  an  electron  with  kinetic  energy  E'. 

6  m  electron  take  off  angle. 

This  model  results  in  a  smoothly  decreasing  ratio  of  the  film  to  substrate  signal  ratio 
as  the  detector  is  moved  from  a  glancing  angle  geometry  to  90**  with  respect  to  the 
substrate  plane.  For  the  case  of  thin  layers,  the  coverage  of  the  substrate  and  film 
thickness  may  be  determined  by  performing  a  fit  to  the  experimentally  observed  data. 

XPS  survey  and  angle  resolved  analysis,  seen  in  Fig.  2,  after  chemically 
cleaning  the  silicon  substrate  indicated  that  there  were  3-4  monolayers  (ML)  of 
oxygen  and  2-3  ML  of  carbon  contamination.  Increasing  the  length  of  each  of  the 
second  HjO  rinses  from  3  to  S  minutes  affected  the  quantity  and  bonding  of  the 
impurities:  the  amount  of  oxygen  doubled,  the  carbon  increas^  slightly,  and  the  C- 
O  (286  eV)*  bonding  was  favored  over  0-Si  (533  eV).*  C-Si  (283  eV),®~*®  C-C,  and 
C-H  (284  eV)“  bonds. 

After  heating  the  wafer  at  700*’C  for  30  minutes,  dissociation  of  C-0  bonds  and 
formation  of  C-Si  bonds  were  observed  (shown  in  Fig.  3).  A  decrease  of  oxygen  and 
I  carbon  impurities  occurred  for  Si(A)  wafers  exposed  to  HjS,  but  not  for  Si(B)  wafers 

which  did  not  undergo  the  H2S  treatment.  Also,  a  larger  percentage  of  the  total 
oxygen  bonded  to  Si(B)  than  to  Si(A). 

The  graphs  in  Fig.  4  show  the  XPS  Zn(A)  signal  was  significantly  weaker  than 
that  from  Z^(B).  indicating  more  zinc  condensed  onto  the  Si(B)  wafers.  Figure  5 
shows  graphs  of  relative  intensities  Zn(A)/Si(A)  and  Zn(B)/Si(B)  vs.  take  off  angle. 

(  The  first.  Zn(A)/Si(A)  smoothly  decreased  which  is  characteristic  of  a  thin  layer  with 

approximately  50  coverage  and  8  A  thickness.  The  second  curve,  Zn(B)/Si(B)  does 
not  decrease  smoothly.  A  high  zinc  aspect  ratio,  (i.e.  three-dimensional  growth), 
could  account  for  the  rapid  decrease  in  substrate  signal  at  lower  angles. 

In  summary,  we  have  observed  that  H2S  pretreatment  of  Si  wafers  affects  the 
surface  impurities,  and  modifies  the  growth  of  thin  Zn  layers.  Removal  of  O  and  C 
)  contamination  is  believed  to  be  responsible  for  an  improved  interface.  Bonding  of  S 

to  the  surface  of  the  wafers  encourages  two-dimensional  growth  in  thin  Zn  layers, 
and  may  be  a  useful  technique  for  quenching  antisite  defects  in  ZnS  grown  from  a 
molecular  source. 
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Fig.  1  Vacuum  system  for  MBE  and  ALE  experiments 
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Fig.  4  XPS  survey  of  Zn  film  on  a)Si(A)  and  b)Si(B). 
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Abstract 

Optical,  electrical,  and  inicrostrucmral  effects  of  Ar  ion  bombardment  and  Ar  incor¬ 
poration  on  thermally  evaporated  Ag  and  A1  thin  films  are  investigated.  The  results 
show  that  as  the  momentum  supplied  to  the  growing  films  by  the  bombarding  ions 
per  arriving  metal  atom  increases,  refractive  index  at  632.8  nm  increases  and  extinc¬ 
tion  coefficient  decreases,  lattice  spacing  expands,  grain  size  increases,  electrical  resis¬ 
tivity  increases,  and  trapped  Ar  increases  slightly.  In  Ag  films,  stress  reverses  from 
tensile  to  compressive  and  in  A1  films  compressive  stress  increases.  In  the  A1  films 
the  change  in  optical  constants  is  due  to  the  presence  of  voids.  The  reversal  of 
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stress  from  tensile  to  compressive  in  Ag  films  requires  a  threshold  level  momentum. 
The  increase  in  electrical  resistivity  is  related  to  the  decrease  in  grain  size  and  incre¬ 
ase  in  trapped  Ar  in  both  types  of  films.  Many  of  these  properties  correlate  well 
with  the  momentum  transferred,  suggesting  that  the  momentum  is  an  important  physi¬ 
cal  parameter  in  describing  the  influence  of  ion  beams  on  growing  thin  films  and 
determining  the  characteristics  of  lAO  thin  metal  films. 
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loa-Assisted  Deposition  of  Thermnlly  Evaporsted  Ag  and  A1  Films 


Chang  Kwon  Hwangbo.  Linda  J.  Lingg.  John  P.  Lehan.  H.  Angus  Macleod, 

John.  L.  Makous.  and  Sang  Yeol  Kim 

John.  L.  Makous  is  with  University  of  Arizona,  Department  of  Physics.  Tucson.  AZ 
85721;  Sang  Yeol  ICim  is  with  Pennsylvania  State  University.  Materials  Research  Lab¬ 
oratory,  University  Park.  PA  16802;  the  other  authors  are  with  University  of  Ari¬ 
zona.  Optical  Sciences  Center.  Tucson.  AZ  85721. 


I.  Introduction 

Ion-assisted  deposition  (lAD)  is  known  to  produce  thin  films  with  characteristics 
that  conventionally  deposited  thin  films  do  not  have.‘'*'^  For  instance,  the  I  AD  of 
dielectric  thin  films  has  shown  that  it  increases  the  packing  density  and  refractive 
index,  reduces  the  moisture  adsorption,  and  improves  the  adhesion  under  certain  mat¬ 
erial  dependent  conditions.*’*’*'^’*’’  lAD  can  also  alter  the  stoichiometry  and  crystal¬ 
linity  of  dielectric  thin  films. For  metal  thin  films,  lAD  enhanced  adhesion  in 
Au.  Ag  and  A1  films  on  glass  substrates.  ”’**’‘*’*’  I  AD  of  Cr  thin  films  resulted  in 
an  increase  in  reflectance  and  reversal  of  stress  from  tensile  to  compressive  when  the 
ion  flux  exceeded  a  critical  value. ‘*  In  addition,  the  non-bulklike  optical  density, 
reduction  of  optical  scatter,  and  control  of  microstructure  in  Cu  films  under  ion  bom¬ 
bardment  were  reported.*^’**’** 
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The  correlation  of  optical  properties  and  microstructure  of  lAD  thin  films  with 
ion  beam  parameters  is  very  important  because  it  provides  insights  for  understanding 
and  producing  these  films.  It  is  known  that  the  ion  energy  and  current  density  are 
two  ion  beam  parameters  which  alter  the  microstructure  of  naturally  growing  thin 
films  and  their  optical  properties.  It  is  not  clear,  however,  what  physical  parameter 
(energy,  momentum  or  some  ftmction  of  the  two)  of  the  ion  beam  affects  growth  of 
thin  films.  ^  Parmigiani  et  al.*  related  the  normalized  energy,  defined  as  the  energy 
delivered  by  the  bombarding  ions  per  arriving  metal  atom,  to  the  optical,  electrical, 
and  microstructural  changes  of  metal  films  in  dual  ion  beam  sputtering.^'  On  the 
other  hand.  Windischmann”  derived  a  model,  based  on  the  linear  cascade  theory  of 
forward  scattering,  such  that  the  intrinsic  stress  of  ion  beam  sputtered  films  depends 
on  the  imparted  momentum  rather  than  energy. 

In  this  paper,  we  investigate  the  optical,  electrical,  and  microstructural  effects  of 
At  ion  bombardment  and  Ar  incorporation  on  thermally  evaporated  Ag  and  A1  thin 
films.  These  effects  and  incorporated  Ar  are  presented  in  terms  of  the  momentum 
supplied  to  the  growing  thin  film  by  the  bombarding  ions  per  arriving  metal  atom. 
We  define  the  momentum  (P)  normalized  by  the  Ar  mass  as  7>/E,  where  7  is  the  (Ar 
ion)-to-<metai  atom)  arrival  rate  at  the  substrate  and  E  is  the  Ar  ion  energy. 

II.  Deposition  Conditions 

Ag  was  evaporated  from  a  Mo  boat  and  Ai  was  e-beam  evaporated  from  an  in- 
termetallic  crucible  in  a  Balzers  760  box  coater.  as  shown  in  Fig.  1.  The  base  pres¬ 
sure  was  IxlO^  mbar.  The  chamber  temperature  was  ambient  and  its  variation 
during  a  single  layer  deposition  was  less  than  S'*  C.  The  deposition  rates  of  Ag  and 
Al  films  were  0.5  nm/sec  and  1.0  nm/sec,  respectively,  as  measured  by  a  quartz  cry¬ 
stal  oscillator.  Prisms,  glass  slides  and  graphite  slabs  were  used  as  substrates  for  the 
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various  post-deposition  measurements.  Just  prior  to  deposition,  substrates  were 
sputter-cleaned  in  vacuo  by  the  Ar  ion  beam  for  S  to  60  seconds,  depending  upon  the 
type  of  substrates  used. 

We  employed  a  3-cm  aperture  Kaufman  hot-cathode  ion  source  with  Ar  gas  as 
the  bombarding  species.  Ar  gas  was  admitted  through  the  ion  gun  until  a  pressure 
of  8x1 0~^  mbar  in  the  chamber  was  obtained.  The  Ar  ion  energies  investigated  were 
200.  300.  and  500  eV  while  the  Ar  ion  current  density  was  varied  from  0  to  20 
fiA/cm*.  These  parameters  resulted  in  the  mass  normalized  momentum  (P)  being 
varied  from  0  to  0.48  VW  for  Ag  films  and  0  to  0.46  >/eV  for  A1  films,  respec¬ 
tively.  The  non-IAO  (conventional)  thin  film  has  zero  P  by  definition. 


III.  Optical  Measurements 

A.  Surface  Plasmon  Resonance 

Optical  constants  at  632.8  nm  wavelength  and  thickness  of  Ag  and  A1  films  were 
obtained  using  the  surface  plasmon  resonance  (SPR)  technique.”  In  the  Kretschmann 
configuration  of  the  SPR.  a  thin  metal  film  is  deposited  on  the  hypotenuse  of  a  right 
angle  prism.”  When  the  angle  of  incidence  of  a  TM  wave  inside  the  prism  is  larger 
than  the  critical  angle,  the  incident  wave  is  totally  reflected  and  an  evanescent  wave 
is  absorbed  by  the  metal  film.  At  a  certain  incidence  angle  and  film  thickness,  the 
component  of  the  wave  vector  of  the  electric  field  parallel  to  the  film  plane  equals  to 
that  of  surface  plasmon  propagating  along  the  metal-air  interface  and  is  thus  coupled 
to  this  surface  plasmon  mode.  When  this  occurs,  the  energy  of  the  incident  wave  is 
absorbed  by  the  surface  plasmon  and  no  wave  is  reflected  by  the  metal  film.  This 
results  in  the  reflectance  showing  a  sharp  resonant  minimum  as  the  incidence  angle 


67 


..  :  '"-  ’■■'>■  - - ”-: - — . :■  ■ ; . ■ '  ■  ’  ' 

-  6  - 

varies  above  the  critical  angle.  The  resonance  angle  and  thickness  which  yields  zero 
reflectance  depend  on  the  optical  constants  of  the  metal  film.”  Optimization  software 
employing  admittance  matching  between  the  incident  medium  and  metal  film  was 
used  to  calculate  the  resonance  angle  and  thickness.”  Optimum  thicknesses  of  about 
50  nm  and  15  nm  were  chosen  for  Ag  and  A1  films,  respectively.  From  the  mea- 
^  sured  reflectance  with  respect  to  the  incident  angle,  the  optical  constants  and  thick¬ 

ness  were  calculated  using  a  non-linear  least  squares  method.”  Reflectance  calcu¬ 
lated  from  the  optical  constants  and  thickness  found  by  the  program  was  in  good 
agreement  with  the  experimental  reflectance  profile.  Measurements  were  performed 
within  5  hours  after  the  films  were  taken  out  of  the  vacuum  chamber  because  the 
SPR  technique  is  sensitive  to  the  oxidation  of  the  metal-air  interface.” 

^  Refractive  index  (n)  and  extinction  coefficient  (k)  of  Ag  films  at  632.8  nm  are 

plotted  in  Fig.  2  in  terms  of  P.  As  P  increases,  n  increases  and  k  decreases.  Ref¬ 
lectance  of  100  nm  thickness  Ag  films  at  normal  incidence  is  calculated  from  these 
I  optical  constants  and  listed  in  Table  1.  At  the  highest  P  (.48  vW),  reflectance  drops 

by  2.5  %  compared  to  the  non-IAO  film.  In  other  words.  LAO  Ag  films  become 
lossy  when  bombarded  with  Ar  ions. 

I  The  n  and  k  of  A1  Hlms  at  632.8  nm  shown  in  Fig.  3  are  scattered  from  a  line. 

Reflectance  of  A1  films  of  100  nm  thickness  at  normal  incidence  is  also  listed  in 
Table  1.  When  P  is  small,  reflectance  does  not  change  much  from  that  of  non-lAD 
(  film.  However,  when  P  exceeds  .23  reflectance  of  lAD  films  decreases. 

B.  Spectroscopic  Ellipsometry 


Recently,  spectroscopic  ellipsometry  (SE)  has  been  used  as  a  nondestructive  tool  to 
depth-profile  multilayer  structures  and  void  fractions  in  thin  films,  in  addition  to  its 
traditional  role,  measuring  the  optical  constants,  dispersion,  and  thickness  of  thin 
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filois.  The  void  fractioa  in  Au  and  amorphous  I  AD  Ce  thin  films  and  depth-profil¬ 
ing  Si,  ZnS,  MgO,  and  ThF^  thin  films  have  recently  been  studied.”*”’*'’”*”  In  this 
work  SE  was  used  to  measure  the  optical  constants,  dispersion,  and  thickness  of  metal 
thin  films  on  glass  substrates  at  an  angle  of  incidence  of  70*  in  the  spectral  range 
1. 5-4.5  eV  in  steps  of  .05  eV.  The  void  fraction  in  A1  films  and  its  oxide  layer 
thickness  were  calculated  by  applying  the  Bruggeman  effective  medium  approximation 
and  linear  regressional  analysis  to  the  SE  data.”*” 

Pseudo-refractive  index  (Op)  and  pseudo-extinction  coefficient  (kp),  the  refractive 
index  and  extinction  coefficient  of  the  film-substrate  assembly,  for  some  Ag  films  are 
plotted  in  Fig.  4  (a)  and  (b)  in  terms  of  photon  energy,  np  of  the  lAD  Ag  films  is 
larger  than  that  of  the  non-IAO  Ag  film  below  the  plasma  frequency  near  3.8  eV  (at 
which  np  equals  to  kp).”  But  in  the  2.5-3.0  eV  spectral  region,  np  of  the  lAD  films 
is  almost  the  same  as  the  non-IAD  film.  In  the  interband  transitions  above  the 
plasma  frequency.  Op  of  the  lAO  films  is  lower  than  that  of  the  non-IAD  film,  kp 
of  the  lAD  Ag  films  is  smaller  than  that  of  the  non-IAD  film  except  in  the  vicinity 
of  the  plasma  frequency.  Also,  the  high  P  film  shows  larger  np  and  smaller  kp  than 
the  low  P  film,  which  is  consistent  with  the  SPR  results.  Reflectance  of  Ag  films  at 
normal  incidence  is  calculated  from  np  and  kp  and  plotted  in  Fig.  S.  Reflectance  of 
these  lAD  films  drops  by  about  1  %  in  the  intraband  transitions,  while  it  drops  by 
about  3%  in  the  interband  transitions.  Around  the  plasma  frequency,  surprisingly, 
the  reflectance  (interband  absorption)  increases  (decreases)  by  about  3%  for  lAD  films. 

Op  and  kp  of  A1  films  are  plotted  in  Fig.  6  (a)  and  (b).  Near  the  l.S  eV  inter¬ 
band  absorption.  Op  of  the  lAD  Al  films  is  smaller  than  that  of  the  non-IAO  A1 
film.  In  the  rest  of  the  measured  spectral  region  (from  about  1.8  eV  to  4.5  eV).  the 
high  P  film  shows  larger  Op  than  the  non-IAD  film  while  the  small  P  film  behaves 
like  the  non-IAD  film,  kp  of  the  lAD  films  is  smaller  than  that  of  the  non-IAO 
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film  in  the  whole  measured  spectral  region.  Reflectance  calculated  from  Op  and  kp 
at  normal  incidence  for  100  nm  thick  films  is  plotted  in  Fig.  7.  Reflectance  of  the 
lAD  films  drops  by  about  2  %  except  at  the  1.5eV  interband  absorption. 

Optical  properties  of  metal  thin  films  in  the  spectral  range  of  intraband  transi¬ 
tions  can  be  described  phenomenologically  by  the  Drude  free-eiectron  model.”  At 
near  IR  frequencies  (ca»l/r)”,  the  real  and  imaginary  parts  of  complex  dielectric 
constants  can  be  written  as  follows: 

a  WpV*^  -  <0  -  Ep*/E*  (1) 

and 

*  cjpVCwH  (2) 

where  represents  core  interband  transitions.  Up  is  the  plasma  angular  frequency 

for  intraband  transitions,  r  is  the  intraband  relaxation  time,  and  &>  is  the  angular  fre¬ 
quency  of  the  light.  Ep  and  E  are  the  energies  of  the  plasma  frequency  and  the  in¬ 
cident  wave,  respectively.  The  intraband  plasma  frequency  of  the  conduction  elec¬ 
trons.  Qp,  is  given  by 

Up-  >/4irNe*/m.  (3) 

where  N  is  the  effective  number  density  of  electrons  participating  in  the  intraband 

transitions,  and  m  and  e  are  the  free-electron  mass  and  charge,  respectively.  Also, 
the  d.c.  conductivity  (a)  in  the  Drude  free-electron  model  is  expressed  by 

<T  m  «p*T/4ir.  (4) 

Using  eq.  (1).  we  can  determine  Up  of  a  metal  film  from  the  slope  of  a  plot  of  vs 
(1/E^.  From  Wp  and  eq.  (3).  N  can  be  obtained. 

The  plots  of  -<i  vs  (1/E*)  calculated  from  the  pseudo-dielectric  constants  for  Ag 
and  A1  films  in  the  intraband  transitions  are  shown  in  Fig.  8  (a)  and  (b),  respec¬ 
tively.  Ep  of  the  lAD  Ag  film  at  P  -  .48  -y/eV  is  7.6  eV,  while  that  of  the  non- 
lAD  film  is  7.7  eV.  For  the  A1  films,  Ep  drops  from  10.7  eV  of  the  non-IAD  film 
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to  10.1  eV  of  the  lAO  filin  at  P  •  .46  >/W.  The  slope  for  the  lAO  films  decreases 
as  P  increases.  The  decrease  of  Ep  means  a  reduction  in  N,  which  has  to  do  with 
the  change  of  microstructure  of  the  metal  film.”*^  The  fraction  of  voids  in  the  A1 
films  is  calculated  by  using  the  effective  medium  approximation  and  linear  regres- 
sional  analysis^'**  and  the  results  are  shown  in  Table  2.  Because  the  A1  films  had 
been  exposed  to  the  atmosphere  for  a  few  months  before  the  SE  measurement  was 
performed,  a  thin  AljOj  layer  on  top  of  the  A1  Him  was  taken  into  account.^  Since 
the  100  nm  thickness  of  A1  film  is  much  larger  than  the  skin  depth  of  A1  (about  10 
nm),  the  incident  wave  does  not  see  the  glass  substrate.  Hence,  bulk  A1  plus  voids  is 
assumed  as  the  substrate.  In  Rg.  9,  the  measured  dielectric  constants  of  the  non- 
lAO  A1  nim  are  in  good  agreement  with  those  calculated  from  the  data  modeled  in 
Table  2.  The  fraction  of  voids  increases  as  P  increases,  while  the  thickness  of  the 
oxide  layer  does  not  seem  to  correlate  with  P. 

Analogous  to  the  reduction  of  interband  absorption  in  the  lAO  Ag  films,  the  1.5 
eV  interband  absorption  of  the  lAD  A1  films  decreases.  Interestingly,  the  1.5  eV  in¬ 
terband  absorption  in  the  A1  films  tends  to  disappear  as  P  increases.  Similar  obser¬ 
vations  for  A1  films  on  cold  substrates  were  reported  by  Bemland.’^  In  his  work,  as 
the  substrate  temperature  decreased  from  198  K  to  25  K.  the  1.5  eV  interband 
absorption  decreased  and  finally  disappeared  at  25  K.  Possible  explanations  were  a 
lack  of  {200)  translational  symmetry  at  25  K  and  an  intermediate  state  between  a  cry¬ 
stalline  and  an  amorphous  state  in  the  A1  film  on  the  cold  substrate.”'”  When  we 
investigated  our  A1  films  with  X-ray  diffraction  (XRD).  we  did  not  observe  a  (200) 
orientation  for  both  non-lAO  and  LAD  A1  films.  So  these  explanations  do  not  seem 
applicable  to  our  films.  If  we  assume  the  reflectance  drop  for  the  intraband  transi¬ 
tions  due  to  the  fraction  of  voids  is  extended  to  the  interband  transitions,  the  reduc¬ 
tion  of  the  1.5  eV  interband  absorption  in  the  lAD  A1  films  should  also  be  due  to 
these  voids.  Generally,  the  presence  of  voids  in  a  thin  film  means  that  the  charac- 
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teristics  of  the  film  deviate  from  the  bulk.^***  Since  the  1.5  eV  interband  absorption 
is  a  characteristic  of  the  bulk,  the  presence  of  voids  may  influence  it  in  such  a  way 
that  the  1.5  eV  interband  absorption  decreases  as  the  void  fraction  increases. 

Given  that  the  fraction  of  voids  correlates  well  with  P.  we  believe  that  the  pres¬ 
ence  of  voids  in  the  Him  induced  by  Ar  ions  is  an  important  factor  which  affects 
the  effective  number  density  of  electrons  in  the  intraband  transitions  and  thus  the 
optical  constants. 


IV.  X-Ray  Diffraction  Analysis 

A  Bragg-Brentano  X-ray  diffractometer  (XRD)  was  used  to  investigate  the 
microstructure  of  polycrystalline  Ag  and  A1  films  such  as  the  preferred  orientation, 
lattice  spacing  and  grain  size.  In  addition,  the  strain  and  average  stress  in  the  film 
plane  were  calculated  from  the  change  in  lattice  spacing. 

For  Ag  films,  the  [Ill]  diffraction  profile  showed  the  strongest  preferred  orienta¬ 
tion.  The  intensity  of  [200]  orientation  is  about  1/7  of  the  intensity  of  the  [111] 
orientation.  Thus  we  will  consider  only  the  [III]  orientation. 

The  situation  for  A1  films  was  more  complicated.  Thin  A1  films  (15  nm.  the 
same  thickness  at  which  the  optical  constants  of  A1  films  were  measured  using  the 
SPR  technique)  showed  no  diffraction  peaks.  Since  these  thin  A1  films  had  been 
exposed  to  the  atmosphere  for  a  while  before  the  XRD  measurements,  an  oxide  layer 
comprising  an  average  of  20  to  30  atomic  %  of  the  film  (as  measured  by  the  Ruther¬ 
ford  Backscattering  Spectrometry)  formed  at  the  air-film  interface.  We  believe  that 
the  thin  oxide  layer  disturbed  the  structure  of  very  thin  AI  films  and  resulted  in  no 
diffraction  profiles.  For  100  nm  thick  Al  films,  the  oxide  layer  comprised  only  an 
average  of  2  to  3  atomic  %  of  the  film  and  so  a  diffraction  profile  was  observed,  but 
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oaiy  th«  [111]  orientation. 

Lattice  spacing  (d„i)  perpendicular  to  the  film  surface  was  calculated  by  the 
Bragg  diffraction  formula.  The  variation  of  dui  for  Ag  films  is  shown  in  Table  3  in 
terms  of  P.  When  P  is  small,  dm  does  not  change  significantly  from  that  of  the 
non-lAD  Ag  film.  When  P  exceeds  .2  vW.  however,  dm  increases  and  levels  off 
at  2.360  A  which  is  slightly  larger  than  the  value  for  Ag  powder  given  by  JCPOS** 
(2.3S9  A).  The  results  for  A1  films  are  shown  in  Table  4.  In  contrast  to  Ag  films, 
dm  of  lAD  A1  films  keeps  expanding  slightly  as  P  increases. 

The  change  in  lattice  spacing  is  related  to  the  strain  (e^)  perpendicular  to  the 
film  plane  as  follows:^ 

-  (d-dj/d,. 

where  d  and  d,  are  the  strained  and  unstrained  lattice  spacings.  respectively. 

If  we  assume  that  the  film  plane  is  elastically  isotropic  and  follow  the  developments 
of  Vook  and  Wiu^  and  Huang  et  al.**.  the  average  stress  (c^y)  and  the  strain  (e,)  in 
the  film  plane  can  be  written  by 

(r.y  -  e,/(2S„') 
and 

+  Su"), 

where  Su'.  Sia',  and  S31'  are  the  elastic  constants  in  the  [111]  orientation  calculated 
by  rotating  the  elastic  constants  for  the  [001]  orientation.**  The  resulting  4  and  e, 
for  Ag  and  A1  films  are  listed  in  Tables  3  and  4.  respectively.  The  average  stress 
in  the  film  plane  is  plotted  in  Fig.  10  as  a  function  of  P.  Non-IAD  and  lightly 
bombarded  lAD  Ag  films  exhibit  tensile  stress,  while  lAD  Ag  films  above  P  -  .2 
•«/eV  show  a  constant  compressive  stress  value  independent  of  P.  Similar  stress  rev¬ 
ersal  has  been  reported  in  lAD  Nb  films,**  dual  ion  beam  sputtered  Ag  films,**  and 
in  other  sputtered  metal  films, *^'**  In  contrast  to  Ag  films,  the  stress  of  lAD  A1 
films  is  compressive  throughout  the  entire  P  range  and  increases  slightly  with  P. 
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It  appears  that  the  compressive  stress  of  lAO  Al  films  correlates  well  with  P 
while  the  stress  reversal  from  tensile  to  compressive  for  lAD  Ag  films  requires  a 
threshold  value  of  P.  Recently.  Windischmann”  proposed  a  simple  scaling  law 
which  relates  the  compressive  stress  for  polycrystalline  thin  films  prepared  by  ion 
beam  sputtering  with  the  momentum.  Also.  Hoffman  and  Gaerttner**  demonstrated 
that  the  important  parameter  for  reversal  of  stress  from  tensile  to  compressive  is  rel¬ 
ated  to  the  momentum  transfer  rather  than  energy  by  investigating  evaporated  Cr 
films  bombarded  with  inert  gases.  Our  results  for  stress  analyses  of  lAO  Ag  and  Al 
films  appear  to  support  Windischmann's  scaling  law  and  be  in  good  agreement  with 
Hoffman  and  Caerttner'  results. 

In  addition  to  strain  and  stress.  XRO  allows  us  to  examine  the  grain  size  in  our 
films.  Grain  size  (Lui)  is  given  by  the  Scherrer  equation.** 

Liu  *  liK/PcosB, 

where  k  is  a  constant  (taken  to  be  1),  X  is  0.1542  nm  (Cu  Ka  radiation),  9  is  the 
Bragg  diffraction  angle,  and  is  the  angular  width  of  the  pure  diffraction  profile  on 
26  scale.  We  corrected  for  instrumental  broadening  by  assuming  Gaussian  profiles 
and  applying  the  John's  correction  method.**  As  shown  in  Fig.  11.  Lm  of  lAD  Ag 
films  is  smaller  than  that  of  the  non-IAD  film  and  it  does  not  change  significantly 
within  the  experimental  error  as  P  increases.  On  the  other  hand.  Lui  of  lAD  Al 
films  decreases  monotonically  as  P  increases  and  at  P  >  .46  >/eV  Lm  of  the  lAD  Al 
film  was  only  about  67  %  of  that  of  the  non-lAD  film. 


V.  Electrical  Resistivity  Measurement 


The  electrical  resistivity  (p)  of  Ag  and  Al  films  at  ambient  temperature  was  mea¬ 
sured  by  the  van  der  Pauw  four-point  method.**  With  the  appropriate  correction 
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factors,  p  of  a  flat  metal  film  of  arbitrary  shape  can  be  measured  if  the  four  contacts 
are  at  the  circumference  of  the  homogeneous  film.  Resistivity  for  Ag  and  A1  films 
are  shown  in  Fig.  12  as  a  function  of  P.  p  of  lAD  films  increases  as  P  increases. 

It  is  known  that  the  electrical  resistivity  of  a  thin  metal  film  is  larger  than  that 
of  the  bulk  and  decreases  as  the  thickness  of  the  film  increases.  Fuchs  developed  a 
free-electron  model  to  explain  this  phenomenon  by  considering  the  size  effects  of  a 
thin  film.’*  In  his  theory,  the  thin  film  is  assumed  to  be  a  homogeneous  slab  and 
electrons  are  scattered  not  only  in  the  bulk,  but  also  by  the  smooth  film  surfaces. 
Real  polycrystalline  films,  however,  have  different  structures  from  the  corresponding 
bulk.  Modified  models,  based  on  Fuchs'  theory  and  taking  into  account  surface 
roughness’’,  grain  boundaries”,  and  charge  carriers  localization’*,  were  later  proposed. 
These  also  are  not  complete  because  the  polycrystalline  films  studied  have  defects, 
dislocations,  and  impurities.  With  this  in  mind.  Sambles”  pointed  out  that  extreme 
care  must  be  taken  in  analyzing  and  interprming  the  data  of  resistivity  measurement. 

Since  Ar  ion  bombardment  on  the  growing  Ag  and  A1  films  introduces  micros- 
tructural  changes,  the  pure  Fuchs'  theory  is  not  valid  for  our  films.  The  RMS  sur¬ 
face  roughness  of  our  films  was  about  1  nm  as  measured  by  a  non-contact  surface 
profiler  (WYKO  TOPO  2D)  and  no  significant  difference  in  the  surface  roughness 
between  non-IAO  and  lAO  films  was  observed.  We  can  thus  exclude  surface  rough¬ 
ness  as  a  factor  in  the  differences  in  p  betw««a  lAD  and  non-IAD  films.  Also,  sur¬ 
face  scattering  due  to  the  film  thickness  effects  is  relatively  negligible  because  the 
grain  size  of  Ag  and  A1  films  is  smaller  than  the  thicknesses  (  about  50  nm  and  100 
nm  for  Ag  and  A1  films,  respectively)  and  so  volume  scattering  dominates. 

When  Parmigiani  et  al.”  applied  Mayadas'”  grain  boundary  scattering  model  to 
ion  bombarded  Ag  films  prepared  by  dual  ion  beam  sputtering,  it  was  found  that  the 
reflection  coefTicient  (Rg)  at  the  grain  boundary  depended  on  the  deposition  condi- 
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tions.  Since  Mayadas'  model  assumes  a  constant  for  films  of  the  same  metal,  the 
different  for  each  film  was  attributed  to  the  different  potential  barrier  at  the 
grain  boundaries  due  to  the  structural  changes. 

In  this  paper,  we  applied  Vancea's  model**  that  the  reduced  effective  number 
density  of  electrons  (N)  participating  in  the  d.c.  conductivity  ir  a  cause  for  the  incre¬ 
ase  of  resistivity  of  metal  thin  films.  Vancea's  model  seems  to  be  a  more  important 
factor  for  the  increase  of  p  of  our  films  than  grain  boundary  scattering  because  we 
observed  the  reduction  of  N  in  the  analysis  of  the  real  part  of  complex  dielectric 
constants  which  is  related  to  the  decrease  of  the  d.c.  conductivity  in  the  free-electron 
model  by  equation  (4).  Also,  in  Vancea's  model,  each  grain  acts  as  a  potential  well 
with  the  boundaries  as  barriers.  Free  electrons  which  pass  through  the  potential  bar¬ 
rier  contribute  to  the  d.c.  conductivity.  As  the  grain  size  decreases,  the  number  of 
barriers  increases,  more  electrons  are  trapped  at  grain  boundaries,  and  the  number  of 
effective  electrons  participating  in  the  d.c.  conductivity  decreases.  This  model  thus 
relates  the  reduction  of  N  and  decrease  in  grain  size  to  the  increase  of  p.  Since  the 
reduction  of  N  implies  an  increase  in  void  fraction  as  discussed  in  the  SE  section, 
the  presence  of  a  larger  fraction  of  voids  in  lAD  films  seems  to  contribute  to  the  in¬ 
crease  of  electric  resistivity  via  Vancea's  model. 

VI.  Rutherford  Backscattering  Spectrometry  Analysis 


The  stoichiometry  of  our  films  was  measured  by  Rutherford  Backscattering  Spec¬ 
trometry  (RBS).  Incident  singly  ionized  helium  (*He*)  ions  with  1.892  MeV  energy 
are  backscattered  by  a  thin  film  sample.  The  number  of  scattered  ions  per  unit  solid 
angle  per  incident  ion  is  measured  and  compared  with  the  known  Rutherford  scatter¬ 
ing  cross  section.  A  simple  calculation  then  yields  the  absolute  number  per  unit  area 
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of  each  atomic  species  in  the  film.’*  The  film  stoichiometry  thus  can  be  measured 
very  accurately.  Graphite  substrates  were  used  because  the  Si  and  O  peaks  in  a 
glass  substrate  would  overlap  the  A1  and  O  peaks  in  the  films  and  thus  reduce  the 
accuracy. 

The  atomic  %  of  Ar  in  Ag  and  A1  films  is  plotted  in  Fig.  13  in  terms  of  P. 
Even  though  the  amount  of  trapped  Ar  is  less  than  2  atomic  %.  it  responds  well  to  P 
in  both  Ag  and  A1  films.  The  average  atomic  %  of  oxygen  in  Ag  films  is  less  than 
.1%  and  therefore  negligible.  A1  films,  on  the  other  hand,  have  most  of  their  oxygen 
at  the  surface  with  the  variation  of  average  atomic  %  of  oxygen  from  run  to  run 
being  less  than  .2  %.  Therefore,  it  seems  that  Ar  ion  bombardment  on  A1  films  does 
not  influence  the  natural  oxidation  of  the  A1  surface  as  discussed  in  modeling  of  A1 
films  in  the  SE  section. 

It  is  known  that  the  trapped  Ar  atoms  in  films  under  ion  bombardment  reside  at 
grain  boundaries.’*  According  to  Vancea's  model.  Ar  atoms  trapped  at  grain  boun¬ 
daries  increase  the  height  of  potential  barriers,  more  electrons  are  localized  at  grain 
boundaries,  and  the  electrical  resistivity  increases.  However,  we  can  speculate  that 
since  Ar  atoms  in  films  are  forced  into  the  film  by  the  ion  gun,  they  all  do  not  nec¬ 
essarily  reside  at  grain  boundaries  and  some  may  reside  within  grains.  If  this  hypo¬ 
thesis  is  correct,  the  small  amount  of  Ar  in  grains  cannot  be  neglected  as  a  factor  in 
the  change  of  electrical  resistivity  and  stress.  Therefore,  in  either  case,  the  inclusion 
of  Ar  in  films  seems  to  affect  electrical  and  structural  properties.  Since  the  electri¬ 
cal  resistivity  and  Ar  concentration  increase  as  P  increases,  we  can  expect  the  electr¬ 
ical  resistivity  to  increase  with  Ar  concentration. 
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The  following  are  observed  in  lAO  Ag  and  A1  films  as  the  momentum  trans¬ 
ferred  to  the  growing  film  by  the  Ar  ion  beam  per  arriving  metal  atom  increases: 

(1)  In  the  intraband  transitions,  the  refractive  index  increases  and  extinction  coeffi¬ 
cient  decreases,  resulting  in  a  decrease  in  reflectance.  The  interband  absorption 
(reflectance),  however,  decreases  (increases).  The  changes  in  optical  constants  are  due 
to  the  reduction  of  the  effective  number  density  of  free  electrons  in  the  intraband 
transitions  and  the  increase  in  the  void  fraction  of  the  films. 

(2)  The  lattice  spacing  perpendicular  to  the  film  surface  expands.  The  stress  changes 
from  tensile  to  compressive  for  Ag  films.  The  compressive  stress  of  A1  films  incre¬ 
ases  slightly.  The  grain  size  decreases  for  both  types  of  films.  The  dependence  of 
stress  on  P  supports  Windischmann's  scaling  law  which  relates  stress  to  the  momen¬ 
tum. 

(3)  The  electrical  resistivity  increases.  This  is  attributable  to  the  reduced  effective 
number  density  of  free  electrons  rather  than  to  grain  boundary  scattering.  The  dec¬ 
rease  of  grain  size  increases  the  number  of  potential  barriers  which  results  in  the 
decrease  of  effective  number  density  of  free  electrons.  This  is  in  good  agreement 
with  the  observed  reduction  of  the  effective  number  density  of  free  electrons  in  the 
intraband  transitions  resulting  from  an  increase  in  the  fraction  of  voids. 

(4)  The  amount  of  trapped  Ar  increases  slightly.  Whether  Ar  atoms  are  at  grain 
boundaries,  in  grains,  or  in  voids,  they  cannot  be  neglected  as  a  factor  which  affects 
the  optical,  electrical,  and  structural  properties. 

From  the  results  of  this  study,  lAD  is  likely  to  make  Ag  and  A1  metal  thin  films 
less-bulklike  by  changing  their  microstructure.  Also,  many  of  their  optical,  electrical, 
and  microstructural  properties  correlate  well  with  the  momentum  transferred  to  the 
growing  film  by  the  Ar  ions  per  arriving  metal  atom,  suggesting  that  the  momentum 
is  an  important  physical  parameter  in  describing  the  influence  of  ion  beams  on  grow¬ 
ing  thin  films  and  determining  the  characteristics  of  lAD  thin  metal  films. 
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Table  1.  Reflectance  of  Ag  and  A1  films  of  100  nm  thickness  at  632.8  nm  at 
normal  incidence  is  calculated  from  the  optical  constants  obtained  from  the  SPR  mea> 
surements. 


Ag  nims 

A1  films 

P  (vW) 

R  (%) 

P  (VeV) 

R  (%) 

0 

98.3 

0 

88.8 

0.18 

97.1 

0.06 

88.7 

0.22 

97.0 

0.23 

88.6 

0.29 

96.7 

0.36 

87.2 

0.30 

96.7 

0.46 

84.6 

0.37 

96.5 

0.48 

95.8 

Table  2.  The  thickness  (t)  of  an  aluminum  oxide  layer  and  the  fraction  (f) 
in  A1  films  are  simulated.  The  thickness  of  film  is  100  nm. 


P  (VeV) 

t  (nm) 

f  {%) 

0 

3.32 

12.4 

0.12 

3.47 

18.2 

0.46 

2.44 

20.7 

of  voids 
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Table  3.  The  lattice  spacing  (dm),  strain  (c,)  perpendicular  the  film  surface,  and 
strain  (<,)  in  the  film  plane  are  listed.  The  thickness  of  Ag  films  is  50  nm. 

Ag  films 


P  (^/eV) 

dm  (A) 

(10-3) 

(10-3) 

0 

2.354 

-2.1200 

2.40 

0.18 

2.354 

-2.1200 

2.40 

0.22 

2.358 

-0.4240 

0.48 

0.29 

2.360 

0.4239 

-4.80 

0.30 

2.360 

0.4239 

-4.80 

0.37 

2.360 

0.4239 

-4.80 

0.48 

2.360 

0.4239 

-4.80 

Powder 

2.359 
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Table  4.  The  lattice  spacing  (d,„),  strain  (Cj)  perpendicular  the  film  surface,  and 
strain  («,)  in  the  film  plane  are  listed.  The  thickness  of  A1  films  is  100  nm. 


A1  films 


P  (>/eV) 

(A) 

0 

2.343 

0.06 

2.344 

0.12 

2.344 

0.29 

2.345 

0.46 

2.346 

Powder 

2.338 

(10-3) 

<,  (10-3) 

2.14 

2.01 

2.57 

2.42 

2.57 

2.42 

2.99 

2.82 

3.42 

3.22 
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Fig.  1.  The  experimental  configuration  of  a  Balzers  760  coating  plant. 

Fig.  2.  Refractive  index  and  extinction  coefficient  at  632.8  nm  of  Ag  films  mea¬ 
sured  by  the  SPR  plotted  as  a  function  of  P. 

Fig.  3.  Refractive  index  and  extinction  coefficient  at  632.8  nm  of  A1  films  mea¬ 
sured  by  the  SPR  plotted  as  a  function  of  P. 

Fig.  4.  (a)  Pseudo-refractive  index  and  (b)  pseudo-extinction  coefficient  of  Ag  films 

measured  by  the  SE  plotted  as  a  function  of  incident  photon  energy. 

Fig.  5.  Reflectance  of  Ag  films  calculated  from  the  pseudo-refractive  index  and 
pseudo-extinction  coefficient. 

Fig.  6.  (a)  Pseudo-refractive  index  and  (b)  pseudo-extinction  coefficient  of  A1  films 

measured  by  the  SE  plotted  as  a  function  of  incident  photon  energy. 

Fig.  7.  Reflectance  of  A1  films  calculated  from  the  pseudo-refractive  index  and 
pseudo-extinction  coefficient. 

Fig.  8.  plotted  as  a  function  of  (1/E*)  for  (a)  Ag  films  and  (b)  A1  films. 

Fig.  9.  Measured  pseudo-refractive  index  and  pseudo-extinction  coefficient  compared 
with  those  calculated  from  the  model  in  Table  2. 

Fig.  10.  Stress  of  Ag  and  A1  films  as  a  function  of  P. 
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Fig.  11.  Grain  size  of  Ag  and  A1  films  as  a  function  of  P. 


Fig.  12.  Resistivity  of  Ag  and  AI  films  as  a  function  of  P. 


Fig.  13.  Ar  concentration  in  Ag  and  Al  films  as  a  function  of  P. 
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Characterization  of  optical  coatings  with  backscattering  spectrometry 


J.  A.  Leavitt.  L.  C.  McIntyre.  Jr..  M.  D.  Ashbaugh.  B.  Dezfouly-Arjomandy. 
and  J.  G.  Oder 


We  present  three  examples  to  illustrate  the  use  of  backscattering  spectrometry 
(utilizing  helium  analysis  beams  in  the  1. 5-5.0  MeV  energy  range)  to 
determine  film  stoichiometry,  areal  density,  and  impurity  levels  in  optical 
coatings:  1)  a  magnesium  fluoride  coating.  2)  a  tungsten/silicon  multilayer  soft 
x-ray  mirror,  and  3)  a  tri-layer  optical  data-storage  film. 


The  authors  are  with  University  of  Arizona,  Department  of  Physics,  Tucson, 
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Backscattering  spectrometry,*  utilizing  analysis  beams  of  *He  ions  in  the 
energy  range  1.5-S.O  MeV,  is  an  underused  technique  for  determining  film 
thickness,  stoichiometry,  and  impurity  levels.  The  basic  technique  is 

illustrated  in  Fig.  1.  Analysis  ions  rebound  from  struck  target  (film)  ions  to 
impinge  on  the  detector  (subtends  solid  angle  Q)  with  energy  characteristic  of 
the  mass  of  the  struck  particle;  further,  the  He  ions  lose  energy  passing  into 
and  out  of  the  film  material.  Energy  analysis  of  the  backscattered  He  ions 
from  a  film  of  uniform  composition  yields  a  nearly  flat-topped  "peak"  for 
each  element  present  in  the  film.  The  film  elements  may  be  identified  by  the 
energies  (E|c,  Eig)  high-energy  sides  of  the  peaks.  Areal  density.  (Nt)}. 

in  atoms/cm^  may  be  determined  for  each  element  from  knowledge  of  solid 
angle  Q.  cross  section  and  integrated  peak  count  Aj  for  Q  incident  analysis 
ions.  Concentration  variations  appear  as  variations  in  the  "flatness"  of  the 
peak  tops.  Note  that  the  average  stoichiometric  ratio.  N<./Ni,.  depends  only  on 
the  ratio  of  integrated  peak  counts,  Aq/A^,  and  knowledge  of  the  scattering 
cross  section  ratio,  Cross  sections  for  scattering  of  2  MeV  *He  on  all 

elements  are  nearly  Rutherford  (Coulomb's  law  holds),  so  Rutherford 
backscattering  (RBS)  is  an  absolute  method  that  does  not  require  the  use  of 
standards.  Average  stoichiometric  ratios  may  be  determined  as  accurately  as 
0.1%  by  acquisition  of  sufficient  data,  while  areal  densities  are  uncertain  by 
about  3%.  Conversion  of  areal  density  (in  atoms/cm^)  to  physical  film 
thickness  (in  A)  requires  knowledge  of  film  density.  Conversely,  if  the 
physical  thickness  is  known,  areal  density  measurement  will  provide  film 
density.  Low-mass  substrates  (such  as  carbon)  are  preferred  so  that  the 
substrate  signal  does  not  interfer  with  the  film  signals. 

We  present  three  examples  to  illustrate  the  p^wer  of  the  technique. 
All  backscattering  data  were  acquired  at  the  Ion  Beam  Analysis  Facility^  at 
the  University  of  Arizona;  a  6  MV  Van  de  Graaff  accelerator  was  the  source 
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of  the  *He  beams. 

The  first  example  is  a  single  layer  of  MgF},  with  impurities,  laid  down 
with  Ar  ion-assisted  deposition  (lAD),  as  part  of  a  study’  of  the  effect  of  I  AD 
on  film  packing  density  and  moisture  resistance.  The  RBS  data  are  shown  in 
Fig.  2;  stoichiometry  and  elemental  areal  densities  are  given  in  Table  1.  Note 
the  high  accuracy  of  the  average  stoichiometry,  the  incorporation  of  Ar  in  the 
film  and  the  replacement  of  F  by  the  ever-present  O.  No  other  impurities 
are  present  in  signficant  amounts.  Measurement  of  physical  film  thickness  (in 
A)  was  made  with  a  WYKO  optical  profiler  and  combined  with  the  RBS  areal 
density  measurement  to  produce  the  film  density.  The  study  showed  that  Ar 
LAD  does  produce  MgF2  coatings  with  increased  packing  density  and  moisture 
resistance. 

The  second  example  is  a  tungsten/silicon  multilayer  soft  x-ray  mirror* 
consisting  of  seven  layers  (S6  A  each)  of  Si  alternating  with  eight  layers  (28.5 
A  each)  of  W  deposited  on  a  carbon  substrate  by  magnetically  confined  dc 
sputtering  with  Ar.  The  sputtering  rates  had  previously  been  calibrated  by 
depositing  films  of  sufficient  thickness  (3000-6000  A)  so  that  their  physical 
thickness  could  be  accurately  measured  with  a  Mirau  optical  interferometer. 
The  RBS  spectrum  is  shown  in  Fig.  3.  Since  the  RBS  depth  resolution  is 
~100  A,  signals  from  the  individual  W  and  Si  layers  are  not  resolved.  The 
film  may  be  treated  as  a  compound  film  of  uniform  composition  for  analysis 
of  W,  Si,  and  Ar;  the  stoichiometry  and  elemental  areal  densities  are  given  in 
Table  2.  Division  of  the  areal  density  (per  layer)  by  the  above  physical 
thicknesses  yields  film  densities  of  2.1  g/cm’  for  Si  and  17.3  g/cm*  for  W. 
The  Si  density  is  in  the  range  for  amorphous  Si,  while  the  W  density  is  10% 
below  that  for  bulk.  The  oxygen  impurity  in  the  film  is  confined  to  a 
surface  layer  less  than  100  A  thick;  if  WOj  of  bulk  density  is  assumed,  then 
the  WOj  layer  is  29  A  thick.  More  detailed  characterization  of  similar  W/Si 
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mtiltilayers  has  involved  utilization  of  low-angle  x-ray  diffraction  data  with 
the  RBS  data  in  an  iterative  procedure  to  determine  layer  thicknesses,  layer 
densities,  and  interface  roughness  for  a  normal-incidence  reflector  of  ~200  A 
radiation.**’ 

The  third  example  is  a  tri-layer  optical  data-storage  film  composed  of 
a  900  A  magneto-optical  rare-earth  transition-metal  (RE-TM)  alloy  film 
sandwiched  between  passivating  layers  (3300  A,  11(X)  A)  of  AijO,. 
Backscattering  spectra  for  incident  *He  analysis  beams  of  energies  1.9  and  3.8 
MeV  are  shown  in  Fig.  4;  areal  densities  and  stoichiometries  for  the  three 
layers  as  determined  from  the  3.8  MeV  data  are  given  in  Table  3.  Note  that 
accurate  determination  of  stoichiometry  and  areal  densities  from  the  1.9  MeV 
data  is  probably  not  possible  due  to  extensive  overlap  of  signals  from 
different  elements  and  different  layers.  The  complexity  of  this  spectnun  gives 
an  indication  of  why  ordinary  2  MeV  *He  RBS  has  not  been  extensively 
applied  to  analysis  of  these  multilayers.  The  3.8  MeV  data  of  Fig.  4  show 
non-overlapping  signals  from  all  elements  in  all  three  layers  except  the  Ar 
signal  in  the  tail  of  the  TM  signal  (this  overlap  is  eliminated  in  4.7  MeV 
data,  not  shown).  It  is,  of  course,  not  necessary  to  completely  eliminate 
overlapping  signals  from  different  elements  in  different  layers.  In  certain 
simple  cases,^  reliable  results  can  be  obtained  by  using  least-squares-fits  to 
deconvolute  the  overlapping  signals.  Use  of  *He  analysis  beams  in  the  high 
energy  (2-5  MeV)  region  does  require  use  of  measured’  cross  sections  for 
scattering  of  *He  by  light  elements  such  as  C,  N,  O,  etc.,  since  these  cross 
sections  are  not  Rutherford  in  this  energy  region.  The  uncertainties  listed  in 
Table  3  are  typical;  the  accuracy  of  the  average  stoichiometric  ratio  could  be 
improved  somewhat  by  taking  more  <'ata.  This  example  illustrates  the  utility 
of  the  higher  energy  *He  analysis  beams. 
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In  sununary.  we  have  presented  three  examples  that  illustrate  the  use 
of  backscattering  spectrometry  for  determining  areal  density,  stoichiometry,  and 
impurity  levels  in  optical  coatings. 

This  work  is  supported  by  the  Air  Force  Office  of  Scientific  Research 
through  the  University  Research  Initiative  Program  and  by  the  Optical  Data 
Storage  Center  at  the  University  of  Arizona.  We  are  grateful  to  P.  Stoss  for 
providing  the  *He  beams,  G.  Van  Zijll,  and  Z.-M.  Yang  for  assistance  with 
acquisition  and  analysis  of  data,  and  to  C.  Falco,  F.  Fernandez.  A.  Macleod, 
and  M.  Messerly  for  providing  samples  and  stimulating  discussions. 
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Table  1.  Results  of  RBS  Analysis  of  the  Mg?,  Optical  Coating. 


Element 

Nt* 

Atomic  % 

Mg 

17.8  ±  0.6 

34.1  t  0.4 

F 

31.2  ±  1.0 

59.9  t  0.4 

0 

2.2  t  0.1 

4.2  ±  0.2 

Ar 

0.93  ±  0.04 

1.80  t  0.10 

W 

0.08  ±  0.01 

<  0.1 

*Nt  -  Atoms/cm^  x  10“ 
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Table  2.  Results  of  RBS  Analysis  of  the  W/Si  Multilayer. 


Element 

Nt* 

Atomic  % 

W 

12.7  ±  0.4 

41.4  ±  0.3 

Si 

17.7  ±  0.6 

57.5  ±  0.3 

Ar 

0.32  ±  0.04 

1.00  ±  0.10 

0^ 

0.16  ±  0.02 

— 

*Nt  -  Atoms/cm‘  x  10‘*. 
■•'o  in  surface  layer  only. 
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Table  3.  Results  of  RBS  Analysis  of  the  Optical  Data  Storage  Tri-Layer. 


Layer 

Element 

Nt* 

Atomic  % 

Ratio 

I 

AJ 

14.8  ±  0.4 

38.9  ±  0.6 

0 

22.3  ±  0.8 

58.6  ±  0.7 

No/Nai  -  1.51  ±0.04 

Ar 

0.9  t  0.2 

2.5  ±  0.3 

2 

RE 

1.13  ±  0.03 

17.4  ±  0.6 

TM 

4.56  ±0.13 

70.3  ±  0.6 

Nre/(Nre+Nti^) 

-  0.199  ±  0.001 

Ar 

0.80  ±0.17 

12.3  ±  1.8 

3 

AJ 

4.8  ±  0.2 

38.2  ±  0.9 

O 

7.5  ±  0.3 

59.8  ±  1.0 

No/Nai  -  1.57  ±  0.06 

Ar 

0.25  ±  0.10 

2.0  ±  0.6 

*Nt  -  Atoms/cm*  x  I0‘*. 
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nOURE  CAPTIONS 

Fig.  1.  Production  of  the  backscattering  spectrum  for  a  two-element 

compound  film  on  a  low-mass  substrate  and  determination  of 
areal  density.  Nt,  and  stoichiometric  ratio.  N^/Nt,. 

Fig.  2.  The  1.9  MeV  *He  backscattering  spectrum  of  a  620  A 

magnesium  fluoride  coating  on  a  carbon  substrate.  The  film 
was  manufactured  by  electron-beam  evaporation  with  Ar  ion- 
assisted  deposition.  Data  acquisition  time  was  five  minutes. 

Fig.  3.  The  1.9  MeV  backscattering  spectrum  of  a  620  A 

tungsten/silicon  multilayer  deposited  on  a  carbon  substrate  by 

dc-sputtering  with  argon.  The  same  data  are  shown  with 

different  vertical  scales  in  a)  and  b).  Data  acquisition  time  was 
five  minutes. 

Fig.  4.  Backscattering  spectra  of  a  5300  A  rare-earth  transition-metal 

optical  data-storage  tri-layer.  Ai20,/(RE-TM)/Af20,.  deposited  on 
a  carbon  substrate  by  Ar  sputtering,  for  incident  *He  analysis 
beams  of  energies  a)  1.9  MeV  and  b)  3.8  MeV.  Data 

acquisition  time  was  less  than  ten  minutes  per  spectrum.  The 
numbers  in  parentheses  are  layer  numbers,  with  layer  1  at  the 
surface. 
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ABSTRACT 

A  continuum  theory  for  the  deposition  and  growth  of  solid  is  presented. 
The  theory  is  developed  in  a  coordinate  independent  manner  and  so  incorporates 
the  fully  nonlinear  physics.  The  evolution  of  the  film  is  modelled  in  three  steps. 
First,  the  2idsorption  of  atoms  in  the  incident  beam  is  modelled  as  a  ballistic 
process.  Second,  the  random  motion  of  the  adatoms  is  treated  as  a  diffusive 
process.  Finally,  sticking  of  adatoms  to  the  film  occurs  as  a  Poisson  process.  The 
resulting  system  of  differential  equations  is  examined  in  several  parameter  limits. 
The  diffusively  dominated  limit  appears  similar  to  zone  I  of  the  str\icture-zone 
model.  Generically  the  surface  slope  develops  discontinuities;  these  ’kinks’  play 
the  role  of  grain  boundaries.  In  the  ballistically  dominated  case  these  kinks  may 
be  advected  along  the  surface  giving  rise  to  columnar  like  microstructure. 
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1.  INTRODUCTION 

While  some  features  of  film  morphology  are  correlated  to  substrate  temperature  by 
Movchan  and  Demchishin’s  (1969)  structxire  zone  model,  prediction  of  film  properties  (such 
as  film  density,  rate  of  growth,  column  size)  &om  deposition  conditions  is  an  actively 
pursued  goal  for  both  theoretical  interest  and  technological  importance.  Past  research  has 
concentrated  on  explaining  the  microscopic  dynamics  of  individual  adatoms.  Computer 
simulations  serially  deposit  spheres  or  disks  onto  previously  deposited  particles  to  simulate 
the  arriving  beam  flux.  The  adatoms  then  move  over  the  surface  according  to  some  hopping 
algorithm.  Continuum-like  models  also  explicitly  include  microscopic  variables  such  as 
adatom  size.  (See  Bartholomeusz  et  al.  (1987)  for  a  review  of  previous  research.) 

In  this  paper  the  properties  on  the  large  scale  are  obtained  by  averaging  over  an  en¬ 
semble  of  individual  particles  passing  through  a  control  volume  on  the  surface.  The  local 
analysis  is  performed  in  a  coordinate  independent  fashion;  consequently  the  local  balance 
then  yields  deposition  equations  in  general  coordinates  that  have  the  advantage  of  being 
fully  nonlinear  on  the  macroscopic  scale.  A  simple  ballistic  argument  yields  a  scattering 
cross-section  which  determines  the  number  of  adatoms  that  are  added  to  the  surface  layer 
and  their  velocity  along  the  srirface.  The  present  paper  attempts  to  track  an  atom  &om  the 
beam  until  it  is  fixed  to  the  film:  it  incorporates  both  a  model  of  the  initial  collision  with 
the  interface  and  the  adatom’s  mobility  on  the  surface.  It  is  hoped  that  this  more  compre¬ 
hensive  viewpoint  will  eventually  explain  a  number  of  morphological  features.  In  this  paper 
some  basic  analysis  of  pertinent  limiting  cases  is  considered;  a  comprehensive  categorizar 
tion  of  behavior  over  the  full  parameter  range  is  left  for  future  efforts.  In  Section  II  the 
physical  mechanisms  are  identified  and  incorporated  into  a  mathematical  formulation.  An 
analysis  of  this  model  is  presented  in  Section  III.  Results  and  comparison  with  experiment 
is  in  Section  IV.  In  the  appendix  details  of  the  ballistic  cross-section  are  presented. 
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n.  FORMULATION 

In  this  section  a  local  analysis  of  the  adatom,  dynamics  is  used  to  deduce  the  deposi¬ 
tion  equations.  The  formulation  incorporates  the  effects  of  incoming  momentum,  adatom 
motion  including  convection  and  diffusion,  and  sticking.  The  incoming  beam  is  assiimed 
to  be  unidirectional  and  uniform.  The  incident  atoms  hit  the  substrate,  and  a  portion  of 
them  join  a  loosely  boimd  surface  layer.  A  ballistic  argument  determines  the  portion  of  the 
beam  flux  which  is  not  reevaporated.  The  flux  rate  onto  a  particular  portion  of  the  surface 
is  then  dependent  upon  the  surface’s  orientation  (Dirks  and  Leamy,  1977).  The  adatoms 
in  this  surface  layer  then  move  by  convection  and  surface  diffusion  until  they  eventually 
stick  leading  to  accretion  of  the  deposited  film.  Bulk  diffusion  is  ignored.  Finally,  adhering 
to  the  siuface  is  modelled  as  a  Poisson  process  and  that  conseqriently  the  deposition  rate 
is  proportional  to  the  ada.tom  concentration. 

A  length  and  time  scale  are  chosen  such  that  the  variation  on  the  substrate  is  initially 
order  unity.  Consider  a  location  on  the  surface,  and  let  the  normal  at  this  point  be 
denoted  by  n  (Figure  la).  The  incident  adatom  flux,  J,  will  be  assximed  to  be  in  a  direction, 
J  and  of  intensity  /, 

(2.1) 

Moreover,  the  magnitude  of  the  incident  velocity,  V,  will  also  be  assumed  to  be  constant. 

As  the  atoms  strike  the  surface,  some  of  them  will  be  reflected  off  the  surface,  and 
some  will  be  adsorbed  into  the  surface  layer.  In  the  appendix  this  process  is  examined  in 
detail;  here  only  two  properties  of  the  adatoms  adsorbed  at  each  point  need  to  be  known 
in  detail,  the  angular  flux  density  of  adatoms,  2md  the  mean  velocity  of  the  adsorbed 
adatoms  V^.  Define  a  polar  angle  xjf  in  the  tangent  plane  and  a  corresponding  unit  vector 
Then  is  the  flux  density  per  unit  radian  and  is  the  mean  velocity  of  the  adatoms 
leaving  in  the  ^  direction  after  collision. 

The  adatom  travels  along  the  surface  in  the  direction  xmtil  it  collides  with  another 
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adatom  or  soir*  local  stirface  feature.  The  occurence  of  this  second  collision  is  a  Poisson 
distributed  process  with  timescale  I/7.  The  length  travelled  before  the  second  collision 
is  then  of  order  V/t-  The  case  when  both  these  scales  are  small  in  comparison  to  the 
characteristic  scales  of  the  surface  is  considered,  that  is 

1  V 

(2.2) 

7  7 

Before  the  adatom  undergoes  a  second  coUirion,  its  velocity  is  determined  by  the 
scattering  angle;  consequently,  this  population  of  adatoms  will  be  called  coherent.  After  the 
adatom  has  imdergone  a  second  collision,  it  will  be  assumed  that  it  is  scattered  randomly. 
This  population  will  be  referred  to  as  incoherent. 

The  equation  of  motion  for  the  density  of  coherent  adatoms,  c,  can  now  be  computed 
by  considering  the  flux  balance  in  a  small  region.  A,  of  the  surface.  Denote  the  edge  as  <7 
and  the  normal  to  the  edge  in  the  surface  as  rf  (flgure  lb).  Then  the  flux  balance  yields 

I  jf  c  --  7  c  id  +  </^(^ .  (2.3) 

Here  the  left  hand  side  of  (2.3)  is  the  time  rate  of  change  of  the  volume  of  adatoms  in 
A.  The  terms  on  the  right  hand  side  are  the  flux  through  the  top  of  A,  the  loss  due  to 
deflection,  and  the  flxix  through  the  edge  of  A,  respectively.  The  adsorbed  flux  density,  7„, 
at  a  given  point  on  the  sxirface,  p,  can  be  computed  by  summing  the  fluxes  scattered  on 
the  surface  through  all  angles, 

j„(^=  rdti^Mp),  (2.4) 

Jq 

and  the  incident  flux  adsorbed  through  the  upper  surface  is  then  given  by  integrating  Jn 
over  A.  The  second  term  represents  the  loss  of  flux  through  deflection.  This  is  modelled  eis 
a  Poisson  process;  consequently  the  time  rate  of  loss  of  coherent  adatoms  is  just  a  constant 
fraction,  7,  of  the  population.  The  last  term  in  (2.3)  represents  the  flux  of  adatoms  through 
the  edge,  er,  computed  by  integrating  the  stirface  flux  due  to  adatoms  moving  in  a  particular 
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orientation  through  a  boundary  element.  The  flux  of  adatoms  moving  in  the  ip  direction 
at  a  point  on  the  botmdary  p(a)  can  be  computed  as 


f  drJ^(p(<T)-rip,t-t')e 
Ja 


(2.5) 


where  the  integrand  is  the  contribution  to  the  flux  from  a  point  a  distance  r  in  the  — 
direction.  The  exponential  term  is  the  probability  that  a  particle  remains  in  the  coherent 
population  after  a  time  t'.  The  particles  which  strike  the  surface  at  a  point,  p  —  rip,  will 
arrive  at  the  location  p  at  time  t'  if 

r  =  V^t‘.  (2.6) 

Using  the  restrictions  (2.2)  the  integral  (2.5)  can  be  evaluated  to  leading  order  as 

^*(M)  «  =  -Vf(?('r),t)J„(P(,r),t),  (2.7) 

Jo  7 

Applying  the  divergence  theorem  to  equation  (2.3),  and  letting  A  shrink  to  a  small  area 
yields  the  flux  balance, 

Ct  =  Jn  -  7C  -  iv.  .  (faJn),  (2.8) 

where  the  average  velocity  Va  has  been  defined  as 


CdiPJ,^V,l,iP 


Jn 


(2.9) 


The  subscript  s  on  the  vector  operators  is  to  indicate  that  they  zu:e  acting  along  the  surface 
defined  by  p.  The  quantities  /„  and  Va  are  determined  by  the  particular  model  used  for 
ballistic  scattering.  These  quantities  are  computed  for  various  cases  in  the  appendix. 
Applying  (2.2)  shows  that  at  leading  order  the  solution  to  (2.8)  is 


c(P7<)  =  -Jn 
7 


(2.10) 


which  can  be  interperted  physically  as  the  incident  flux  times  the  average  time  to  deflection 
of  an  adatom. 

"  1,6 


Equation  (2.8)  describes  the  evolution  of  c  due  to  the  coherent  motion  of  ballistic 
I  scattering.  The  incoherent  motion  of  adatoms  caused  by  random  scattering  on  the  surface 

must  also  be  included.  It  will  be  asstuned  that  after  an  adatom  in  the  surface  layer  is 
deflected  that  the  orientation  of  its  velocity  is  determined  randomly.  Consequently,  the 
^  motion  of  the  particle  can  be  modelled  as  a  difEusive  process. 

The  rate  at  which  particles  leave  the  coherent  adatom  concentration  is  7c;  they  then 
become  part  of  an  incoherent  adatom  density  d.  Adatoms  are  removed  &om  the  incoherent 
1^  state  by  sticking  to  the  surface;  this  is  again  assumed  to  be  a  Poisson  process  with  a  mean 

lifetime  I/7'.  The  equation  of  motion  for  d  is  now  given  by 

^  dt  =  7c  -  7'd  +  i/V^d,  (2.11) 

where  is  a  diffusion  constant. 

As  adatoms  are  deposited,  the  surface  grows  outward  along  its  normal  in  proportion 
to  the  rate  of  deposition, 

Pt  =!  7'dn.  (2.12) 

I  This  equation  together  with  (2.8)  and  (2.11)  specify  the  evolution  of  the  surface. 

To  summarise,  the  adatom  motion  on  the  surface  is  due  to  a  component  which  depends 
on  the  direction  of  the  incident  flux  plus  a  randomly  oriented  diffusive  component.  The  flux 
^  rate  of  particles  which  ultimately  stick  on  the  surface  is  given  by  Jn  and  will  be  determined 

in  the  appendix  by  considering  the  collision  of  an  incident  particle  with  the  surface.  The 
model  equations  governing  the  process  of  film  growth  are 


Ct  =  7n-7C--V,.(f.7n), 

7 

(2.13a) 

d,=7c-7'd+./V2d, 

(2.136) 

i 

J 

1 

[ 

Pt  =  y'dh, 

(2.13c) 
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where 

p  s:  location  of  the  surface 
n  =  surface  normal 
c  =  coherent  adatom  volume  density 
d  =s  incoherent  adatom  volmne  density 
Jn  —  adsorbed  adatom  flux 
Va  s:  average  adsorbed  adatom  velocity 
1/7  =  coherent  adatom  d^ection  rate 
1/7' s  incoherent  adatom  sticking  rate 
u  =  incoherent  SMlatom  diffusion  constant 

HI.  ANALYSIS 


A.  Conservation  of  integrated  flux 

In  this  section  it  is  shown  how  the  integrated  adsorbed  incident  flux  relates  to  the 
integrated  flux  in  the  coherent  and  incoherent  adatom  populations  and  the  increase  in 
volume  of  the  film.  By  considering  the  total  incident  flux  onto  the  surface  it  is  possible 
to  determine  an  expression  for  the  rate  of  growth  of  the  film.  Consider  a  domain,  A,  with 
periodic  or  no-flux  boundary  conditions.  Define  the  volume  of  previously  deposited  film 
below  as  V.  Similarly,  define  the  volume  of  incoherent  and  coherent  adatoms  on  the 


surface  as 


c  d5  =  C, 


(3.1a) 


J  ddS^V, 

respectively,  where  dS  is  an  element  of  stirface  area. 


(3.15) 


Integrating  equation  (2.13a)  over  the  surface  relates  the  total  adsorbed  incident  fliix, 


=  jjndS, 
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to  the  evolution  of  C, 


Cf\‘'rC=^J.  (3.3) 

Similarly  integrating  (2.13b)  relates  the  evolution  of  C  and  P, 

Vt  +  =  7C.  (3.4) 

Finally  the  time  rate  of  change  of  the  volume  of  the  film  can  be  related  to  the  sticking 
of  incoherent  adatoms  by  dotting  (2.13c)  into  the  surface  normal  and  integrating  over  the 
surface, 

Vt  =  7'C.  (3.5) 

The  conservation  of  total  mass  can  now  be  obtained  by  adding  (3.3, 4, 5)  together, 

dt(C  +  V  +  V)=^J.  (3.6) 

Equations  (3.3-3.6)  demonstrate  how  the  incident  flux  cascades  through  the  various  stages 
of  the  modeL  They  allow  the  mean  growth  rate  of  the  thin  film  to  be  determined  as  a 
function  of  the  integrated  adsorbed  flux. 

B.  DifiPusion  Dominated  Dynamics 

In  this  section  the  case  when  diffusion  dominates  the  evolution  of  the  interface  will  be 
considered.  The  mean  diffusion  length,  of  an  incoherent  adatom  can  be  approximated 
by  considering  (2.13b), 

(,.T, 

Suppose  the  surface  microstructure  varies  on  a  scale  that  is  small  compared  to  the  diffusion 
length,  that  is  >>  1.  In  this  case  diffusion  will  be  the  dominant  surface  process.  Equation 
(2.13b)  will  relax  towards  an  equilibrium  where  V,d  is  zero,  that  is  a  state  where  the 
concentration  of  d  is  uniform  on  the  surface.  This  implies 


I#'  .i-l  .1 


Solving  (3.4)  now  yields 

d{t) «  4  /* 

^ — oo 

or  applying  (2.10), 

d(t)  «  4  /* 

OO 

In  the  limit  when  JtjJ  ">  7'  this  reduces  to 


which  is  simply  the  average  adsorbed  flux. 
If  we  define  a  new  time  scale, 


equation  (2.13c)  takes  a  particularly  simple  form 

^  A 

Pr  -  n. 


(3.9) 


(3.10) 


(3.11) 


(3.12) 


(3.13) 


Physically,  the  surface  grows  at  a  constant  rate  in  the  direction  of  its  normal.  In  the  limit 
defined  by  (3.11)  the  rate  of  growth  is  simply  the  average  adsorbed  flux  rate.  Equation 
(3.13)  is  an  eikonal  equation  that  commonly  arises  when  considering  the  rays  of  a  wave 
equation  and  can  be  solved  implicitly  with  the  method  of  characteristics  (Whitham,  1974). 

It  will  be  convenient  to  introduce  a  coordinate  system  for  equations  (2.13);  define  the 
surface  elevation,  p,  measxired  above  some  plane  with  cartesiein  coordinates  Xi  =  (x,l/), 
then 

P=ix,y,p{x,y,t)).  (3.14) 


Equation  (2.13c),  which  gives  a  Lagrangian  description  of  the  interface,  can  be  converted 
to  2U1  Eulerian  description.  The  normal  at  a  point  on  the  surface  is  given  by 


n  = 


(i+pi+rf)''’ 


(~Pr>  —Py?  !)■ 


(3.15) 
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If  we  denote  by  a  dot  the  time  derivative  in  the  Lagrangian  coordinate,  the  chajTi  rule 
yields 

Pr  =  P  +  P*x  +  Pj,y,  (3.16) 

which  after  using  (3.13)  yields 

p^  =  (l+p|+pJ)^/2_  (3.17) 

Following  Whitham  (1974)  an  implicit  solution  to  the  initial  value  problem  (3.17)  can 
be  foimd.  It  will  be  convenient  in  what  follows  to  use  index  notation  with  the  summation 
convention;  define  the  gradient  di  =  (5*,5y).  Suppose  initially, 


p(xj)  =  P(x,y)  at  r  =  0. 


(3.18) 


Then  at  a  later  time  t  the  solution  at  a  point  =s  (x,  p)  can  be  found  in  terms  of  P(x,  y). 


p(i,^,r)*P(x,y)  + 


r 

(i  +  dtPdtPyi^' 


where  x,-  is  defined  implicitly  as  a  function  of  x,-  and  r. 


(3.19) 


rdiP 

(l  +  afcP5feP)V2' 


(3.20) 


The  implicit  solution  (3.20)  allows  some  explicit  special  solutions;  in  particular,  a  plane 
will  grow  with  a  unit  velocity  normal  to  its  surface  md  a  spherical  cap  will  remain  sphericsd 
with  a  constantly  increasing  radius. 

Because  (3.20)  is  an  implicit  solution  corresponding  to  characteristics,  it  admits  the 
possibility  of  multivalued  solutions.  To  determine  if  and  when  this  happens  it  is  necessary 
to  decide  if  (3.20)  provides  a  continuous  bijection  from  x,-  to  ij.  This  will  be  the  case  if 
the  Jacobian,  Q,  between  the  two  coordinates  is  adways  non-zero.  Using  (3.20)  to  compute 
the  Jacobian  yields  the  determinant 


<3  = 


dxj 

dxj 


(3.21) 
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where  the  Sii  is  the  Kronecker  delta  function  and 


. ,  diiP{i+ dkPdkPy^^  -  diP  dkP  dikP 
'  (1  +  dkPdkPyn 


(3.22) 


To  understand  (3.22)  requires  a  basic  knowledge  of  differential  geometry  (Do  Canno,  1976). 
The  determinant  and  trace  of  Mij  are  the  Gaussian  curvature  and  twice  the  mean  curvature 
respectively.  This  implies  that  the  eigenvalues  of  Mij  are  exactly  the  principal  curvatures, 
Ki,  K2  of  the  stnface.  Geometrically  1/ki,  I//C2  correspond  to  the  maximiiTn  and  minimum 
radius  of  the  osculating  circle  at  a  point  on  P.  From  (3.21)  it  follows  that  Q  will  vanish 
when 

{*<}  \«i  >^2/ 

Note  that  if  ki,K2  are  never  greater  than  zero,  the  initial  condition  is  convex,  and  the 
solution  never  breakdowns.  However,  if  the  initial  condition  is  non-convex  anywhere,  at 
least  one  of  the  curvatmres  will  be  positive  and  the  solution  will  eventually  break  down. 

The  form  of  this  breakdown  can  be  deduced  by  notin  ’  "t  the  gradient  of  (3.17)  is 
a  first  order  semi*linear  hyperbolic  system  for  dip.  Conse  '  y  ‘  e  gradient  will  develop 
shocks  (discontinuities)  and  p  will  develop  ’kinks’  (disconimoii  ;  the  derivative  of  p). 
In  two  dimensions  these  kinks  will  occtir  at  isolated  points;  in  thr  *  dimensions  they  will 
occur  on  curves  on  the  sxirface.  A  description  of  the  motion  of  these  kinks  is  needed  for 
(3.17)  to  provide  a  complete  description  of  the  dynamics.  Fortunately,  phy3ic2d  insight 
provides  a  unique  solution  for  the  motion  of  the  kinks.  The  equations  aie  invariant  under 
multiplying  Xj,  p  and  r  by  a  constant;  from  this  it  follows  that  the  velocity  is  dependent 
only  on  the  slope  of  the  tangent  plane  on  each  side  of  the  kink.  The  isotropy  of  space  now 
implies  that  the  kink  must  propagate  along  the  angle  bisector  of  the  dihedral  angle  formed 
by  the  tangent  planes  at  the  kink  (figure  2).  This  rule  together  with  the  characteristic 
description  given  by  (3.19,20)  is  sufficient  to  describe  the  dynamics.  A  sample  numericsJ 
evolution  of  a  two  dimensional  initial  condition  is  shown  in  figure  3. 

In  three  dimensions  other  phenomena  can  occur;  the  kinks  correspond  to  boimdaries 
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between  grains.  If  the  evolutionirom  an  initial  condition  corresponding  to  a  rough  substrate 
is  considered,  various  bifurcations  occur.  The  initial  appearance  of  the  kink  will  still 
occur  at  a  point  corresponding  to  the  highest  curvature  of  a  concavity.  After  formation 
the  kink  will  become  a  line  of  discontinuity.  Pairs  of  lines  of  kink.s  can  collide  creating 
triple  points  where  three  grains  meet  (figure  4).  These  triple  points  can  also  occiir  by  the 
appearance  of  a  new  discontinmty  in  curvatxue  transverse  to  am  existing  line  of  kinks.  The 
evolution  of  this  triple  conjunction  is  again  easily  prescribed;  it  must  propagate  along  the 
lirn*  given  by  the  intersection  of  the  angle  bisectors  of  the  three  dihedral  angles  created  by 
the  tangent  planes  at  the  point.  In  thi«  way  an  initially  variable  topogra^y  will  lead  to  a 
system  of  kinks  coimected  at  triple  points,  resembling  the  networks  of  grains  seen  in  thin 
film  growth. 

C.  Ballistic  Dominated  Dynamics 

A  second  special  case  of  interest  is  when  ballistic  effects  dominate  the  effect  of  dif> 
fusion.  The  limit  considered  here  is  when  the  diffusion  length  is  much  smaller  than  the 
characteristic  length  in  the  problem,  that  is  ^  1.  In  addition,  it  will  be  assiimed  that 
the  deposition  time  is  short  compared  to  the  the  surface  growth  rate,  that  is  I/7'  1.  In 

this  case,  at  leziding  order 

■y'd  faycfti  Jn.  (3-24) 

Equation  (2.13c)  now  reduces  to 

pi  =  7„n.  (3.25) 

It  will  be  convenient  to  consider  a  two  dimensional  evolution  for  equation  (3.25). 
Following  (3.14)  define  the  height  of  the  surface  as 

P=  (x,y,p(x,t)),  (3.26) 

where  the  dependence  of  p  on  y  has  been  suppressed.  In  addition,  it  is  more  illuminating  to 
introduce  defined  as  the  angle  the  surface  normal  makes  with  the  direction  of  incoming 
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flux.  It  will  be  assumed  that  J  x  —z,  which  implies 


=  (3.27) 

Equation  (3.27)  and  the  restriction  that  4  is  positive  when  the  slope  of  p  is  positive  specify 
4  in  terms  of  p. 

The  assumption  J  =  —  z  is  not  at  all  restrictive;  the  isotropy  of  our  formulation 
guarantees  that  looking  at  vertical  flux  on  a  slanted  surface  is  identical  to  looking  at  flux 
incident  from  the  same  angle  on  a  flat  surface.  Using  (3.27),  p  can  be  eliminated  in  favor 
of  ^  in  equation  (3.25), 

4i  =  C03^{4)dx{Jn3ec{4)),  (3.28) 

From  the  ballistic  results  in  the  appendix,  it  is  clear  that  depends  only  on  (^,  so  equation 
(3.28)  can  be  rewritten  as 

4t  +  *  0.  (3.29) 

where 

=  -{{^n)<t,cos{4)  +  Jn3in{4))  (3.30) 

Equation  (3.29)  is  a  wave  equation.  Any  initial  condition  with  a  variation  in  4  ^nd 
U  should  cause  the  formation  of  shocks.  Note  that  a  shock  corresponds  to  a  discontinuity 
in  4i  from  (3.27)  it  follows  that  this  corresponds  to  a  kink  in  the  surface  of  the  film  as 
described  by  p. 

Figure  5  shows  U  graphed  for  various  ballistic  models.  Physical  symmetry  implies 
that  the  velocity  is  an  antisymmetric  function  of  4-  is  cxirious  to  note  that  in  the  three 
dimensional  hard  sphere  case  that  U  is  zero,  and  equation  (3.29)  corresponds  to  growth  at 
a  constant  rate  towards  the  source.  For  the  other  three  dimensional  cases  note  that  U  is 
negative  for  ^  >  0,  and  decreases  to  a  minimum  and  tends  to  zero  as  4  tends  to  ir/2.  For  the 
two  dimensional  case  hard  sphere  U  changes  sign;  this  will  lead  to  qualitatively  different 
dynamics.  This  should  be  taken  as  a  caution  ageiirtst  using  two  dimensional  ballistic  models. 


•*»:  v'v; 


When  equation  (3.29)  becomes  invalid  due  to  shock  formation,  it  is  essentially  because 
of  the  formation  of  high  gradients  in  <f>.  Because  derivatives  become  large,  higher  order 
terms  in  the  equation  become  important.  This  will  lead  to  a  regularization  of  the  wave 
equation. 

The  next  higher  order  contributions  due  to  corrections  to  (3.24)  can  easily  be  com¬ 
puted.  Substituting  (3.24)  into  (2.20)  3rields 


yd  +  vV\d  -dt-ct-  i V. .  (f.7„) 

which  resubstituting  (3.24),  and  using  (3.27)  and  (3.29)  yields 

cos{<f>) 


(3.31) 


«  A  +  [(i + i)  vij.u  -  j„) 


<^x  -  ^C03((^)dj:  {cOs{(f>)dj;  ( J„))  , 


(3.32) 


The  second  and  third  terms  on  the  right  hand  side  represent  corrections  to  the  approximar 
tion  made  by  (3.24).  By  balancing  the  leading  order  term  Jn  against  each  of  these  terms 
the  shock  width  produced  by  each  of  the  regularizations  can  be  deduced.  Note  that  the 
second  term  has  a  single  x  derivative;  it  corresponds  to  a  dissipative  regularization.  This 
term  becomes  important  on  the  shock  lengthscaie  L,  which  is  dimensionally  approximated 
by 

+  (3-33) 

Note  that  in  general  this  scale  is  on  the  order  the  ballistic  scattering  distance;  conse¬ 
quently  the  arguments  used  to  derive  (2.13a)  are  no  longer  valid.  If  this  term  leads  to  the 
regularization,  a  detailed  examination  of  the  ballistic  effects  in  this  region  will  be  needed. 

If  the  third  term  in  (3.32)  leads  to  the  regularization  of  the  shock,  the  shock  width 
L'  will  be 

/ 


which  is  exactly  the  diffusion  length!  Note  that  the  term  has  two  x  derivatives  correspond¬ 
ing  to  a  dispersive  regularization. 
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To  tinderstand  qualitatively  the  dynamics  of  these  two  regularization  the  special  case 
^  is  nearly  constant  is  considered.  Let 

<i>{x,  t)  =:  (pQ  +  €^ix,  t)  (3.31) 

where  e  -C  1.  Consider  a  change  of  variables  into  a  frame  moving  at  the  phase  velocity 
U{4>q),  and  on  a  slow  time  scale, 


=  r  - 


(3.35a) 


r  =  e<, 


Then  equation  (3.29)  can  be  expanded  as 


(3.356) 


+  F{^q)^xx  +  C(^o)^ix*  +  . . .  =  0.  (3.36) 


where 


nM = [(i + i)  , 


G(^o)  ==  — 7COs’(^o)(>^n)ffr' 
*7 


(3.37a) 

(3.376) 


In  the  case  when  L,  ^  Ld  the  $xz  term  will  be  the  dominant  regularization.  In  this 
case  equation  (3.36)  is  Burger’s  equation.  When  F(^o)  is  negative,  the  shock  will  have  a 
width  of  0(X,/c),  and  characteristics  will  flow  uniformly  into  the  shock.  This  leads  to  the 
decay  of  the  shock  and  for  large  time  the  solution  tends  to  a  imiform  constant.  When  F(<^o) 
is  positive  the  dissipation  is  acting  in  a  negative  sense,  and  the  equation  is  ill-posed.  Any 
small  disttirbance  will  tend  to  focxis  and  blow  up  in  a  finite  amoimt  of  time.  Presumably 
this  blow  up  might  saturate  if  higher  order  terms  are  included.  However,  any  model  with 
this  asymptotic  limit  presumably  exhibits  highly  chzwjtic  dynamics,  similar  to  what  is  seen 
in  the  Kuramoto-Sivashinsky  equation. 

In  the  case  when  L,  C  Ld  the  $***  term  will  be  the  dominant  regularization.  In 
this  case  equation  (3.36)  is  the  Korteweg-deVries  (KdV)  equation.  KdV  is  well  known  for 
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exhibiting  soliton  type  solutions;  lump-like  objects  which  are  exponentially  localised  and 
translate  without  change.  The  width  of  the  soliton  would  be  0(Ld/€).  Note  that  in  terms 
of  the  surface  elevation  p  a  soliton  would  be  a  step  separating  two  regions  of  constant  slope 
(figure  6).  As  c  increases,  presumable  this  would  correspond  to  a  narrow  region  of  width 
Ld  separating  two  flat  regions,  and  translating  with  speed  of  U(4>q)- 

IV.  RESULTS  AND  CONCLUSIONS 

In  this  paper  a  formalism  is  described  for  deriving  a  continuum  model  of  thin  film 
growth.  The  deposition  is  broken  into  three  processes:  a  ballistic  scattering  model  is  used 
to  determine  adsorption  onto  the  surface;  the  random  walk  of  the  adatoms  on  the  surface 
is  described  by  a  difiusion  equation;  and  sticking  of  the  adatoms  to  the  surface  is  described 
as  a  Poisson  process.  Because  the  model  is  formulated  in  a  coordinate  independent  fashion, 
the  local  dynamics,  which  are  basically  linear,  lead  to  a  fully  nonlinear  model  of  the  large 
scale  evolution.  This  supersedes  the  earlier  work  of  Lichter  and  Chen  (1986)  who  only 
considered  linear  perturbations  to  a  flat  film. 

The  proportion  of  incident  atoms  adsorbed  can  be  computed  from  the  ballistic  models 
described  in  the  appendix.  This  in  turn  allows  the  film  growth  rate  to  be  computed  and 
compared  to  experiment. 

In  the  case  when  sxirface  diffusion  is  the  dominant  process,  an  eikonal  equation  has 
been  derived  to  describe  the  growth  of  the  film.  The  film  locally  grows  at  a  constant 
rate  normal  to  its  surface.  The  eikonal  equation  has  a  special  solution  of  note;  a  sphere 
or  spherical  cap  will  grow  in  a  self  similar  fashion  with  a  constantly  growing  radius.  This 
means  a  local  bulge  in  the  substrate  leads  to  a  cone  like  solution.  In  addition  an  imperfection 
seeding  the  medium  may  also  cause  a  conical  growth  pattern. 

The  regime  in  which  this  behavior  is  expected  is  given  by  the  mean  diffusion  length, 
Ld  1.  This  parameter  can  be  expected  to  depend  exponentially  on  the  inverse  of  the 
temperature  (Burton  et  al,  1951);  this  implies  that  low  substrate  temperattires  will  cor- 
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Tespoad  to  a  large  diffusion  length.  The  cono'shaped  microstnictxire  with  spherical  caps 
corresponds  to  characteristic  morphology  of  zone  1  of  Movchan  and  Demchishin’s  structure 
zone  model. 

In  addition,  the  eikonal  eq\iations  predict  the  formation  of  kink  like  discontinuity  of 
slopes.  A  network  of  such  kinks  is  expected  for  an  initially  rough  substrate.  This  network 
of  grain  boundaries  seems  characteristic  of  experimental  observations  (Bartholomeusz  et 
al.,  1987). 

The  model  is  also  examined  in  the  case  when  ballistic  advection  is  much  greater 
than  diffusion.  Here  the  characteristic  equation  describing  the  film’s  growth  is  a  first-order 
wave  equation.  For  three  dimensional  hard  sphere  scattering,  the  phase  velocity  of  the 
wave  equation  is  zero;  consequently  the  solution  will  always  grow  in  a  self-similar  fashion 
towards  the  incident  flux.  However,  energy  loss  to  the  atomic  lattice  is  expected  (Lennard- 
Jones,  1937)  after  the  substrate  has  been  coated  with  the  incident  species.  This  energy 
loss  is  modelled  by  using  a  ballistic  cross-section  with  a  spherical  particle  with  a  finite 
effective  mass.  This  leads  to  a  non-uniform  phase  velocity  for  the  governing  wave  equation. 
This  non-uniformity  will  lead  to  shock  formation  which  once  again  manifests  itself  as  a 
discontinuity  in  the  surface  slope.  These  shocks  will  propagate  across  the  film;  their  tracks 
will  leave  columnar  like  deposition. 

The  details  of  the  shock  structure  depend  upon  the  higher  order  corrections  to  the 
wave  equation.  Two  cantfidates  for  regularizing  the  wave  equation  have  been  identified. 
When  surfeure  diffusion  is  small  compared  to  the  coherent  scattering  distance  of  an  eulatom, 
the  regularization  is  dissipative,  emd  a  Burger’s  type  equation  is  an  appropriate  model  for 
the  evolution.  In  this  case  shocks  are  localised  and  propagate  (WHiitham,  1974).  When  the 
diffusion  length  is  longer  than  the  coherent  scattering  length,  the  appropriate  regularization 
is  dispersive.  Here  the  canonical  governing  equation  is  Korteweg-de Vries.  This  allows  the 
possibilty  of  soliton  like  solutions  corresponding  to  a  step-like  surface  geometry.  These 
equations  undoubtable  exhibit  a  rich  and  varied  dynamics  depending  on  the  parameter 
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regimes  examined  and  consequently  they  are  left  for  future  study. 

Note  that  in  the  difEiision  dominated  limit,  a  flat  surface  will  grow  in  a  direction  normal 
to  its  surface  independent  of  the  direction  of  the  incident  flux.  In  the  ballistically  dominated 
limit,  when  three  dimensional  hard  sphere  scattering  is  considered,  the  surface  grows  at 
a  constant  rate  in  the  direction  of  the  incident  fltix  independent  of  the  surface  normal.  It 
seems  logical  to  assume  that  when  these  two  eflects  are  comparable  the  surface  will  grow 
at  some  angle  between  the  normal  and  the  direction  of  incident  flux.  Experiments  observe 
that  the  angle  of  columnar  growth  in  zone  2  of  Movchan  and  Demchishin’s  structure  zone 
model  is  intermediate  between  these  two  extremes,  and  follows  the  empirical  relationship 
known  as  the  tangent  rule  (Dirks  and  Leamy,  1977).  Our  model  suggests  that  this  is  effect 
may  be  due  to  a  balance  of  surface  diffusion  and  ballistic  effects. 

APPENDIX.  BALLISTIC  COMPUTATION  OF  FLUX  CROSS-SECTIONS 

In  this  appendix  a  formalism  for  using  the  ballistics  of  a  two  sphere  interaction  to 
derive  a  flux  and  velocity  cross-section  on  the  surface  of  the  film  is  presented.  Consider  an 
incident  spherical  particle  striking  a  surface  with  a  rough  atomic  scale.  Assume  that  the 
orientation  of  the  siuf2u:e  on  the  atomic  scale  is  beisically  random.  This  case  corresponds 
to  the  incident  particle  being  scattered  from  a  spherical  body  because  the  orientation  of 
the  normal  to  a  sphere  is  isotropic  in  space.  In  addition,  the  film  will  have  a  certain  degree 
of  elasticity.  This  will  be  modelled  by  defining  a  mass  ratio,  M,  between  the  mass  of 
the  incoming  adatom  smd  the  effective  mass  of  the  particle  struck  in  the  film.  The  limit 
M  -*  oo  correponds  to  scattering  off  a  fixed  hard  sphere,  whereas  M  =s  1  corresponds  to 
striking  an  tmboimd  adatom.  The  study  is  restricted  to  M  between  these  limits. 

Using  conservation  of  momentum  and  energy  allows  us  to  relate  the  exit  velocity  Ve 
and  scattering  angle,  as  a  function  of  the  incident  velocity,  V,  and  impact  angle,  ^ 
(figure  7).  The  velocity  ratio  can  be  computed  as 

Ve  (iV/^  -  -  cos{^') 


1-HM 


(Al) 
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The  relation  between  the  impact  angle  and  scattering  angle  can  only  be  written  implicitly, 

cos(2^)-^ 

The  incident  flux  is  assumed  to  be  arriving  vertically,  that  is 

J  =  -z.  (A3) 

The  scattering  from  a  unit  sphere  at  the  origin  will  be  computed  Define  spherical  polar 
coordinates,  where  a,j3  are  the  azimuthal  angle  measured  &om  x  zmd  the  declination 
from  z  respectively.  The  projection  along  the  z-  axis  is  the  unit  circle.  Define  the  radial 
component,  r,  in  the  x-y  plane, 

r  =  (x^  +  y*)*/*  =  sin/?,  (.44) 

then  the  differential  flux  density,  dJ,  is  given  by  the  flux  times  the  area  element  in  the 
plane  divided  by  the  total  area  of  the  projection. 


dJ=  —  rdr  da.  (AS) 

Notice  that  the  integral  over  the  entire  cross*section  simply  yields  J.  Substituting  (A4) 
into  (A5)  yields  the  cross-section  as  a  function  of  a,  /? 

J 

dJ  =  —sin{^)cos(0)  da  d^.  (A6) 

JT 

To  compute  the  flux  scattered  out  at  various  eingles,  the  scattering  angle  must  be 
known;  aasiune  that  the  particle  leaves  in  a  direction  a\  /?'  .  For  the  geometry  considered 
here  the  aizimuthal  angle  remains  the  same  (a’  =  a)  zmd  equation  (A2)  relates  the  scatter¬ 
ing  angle  to  the  impact  angle  The  differential  cross-section  for  the  scattered  particles 


can  now  be  written  zis 


dJ  =  ^sm(/?)cos(^) 
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where  it  is  understood  that  /9  is  a  function  of  P*. 

Finally,  it  is  desirable  to  compute  the  differential  cross  section  in  a  spherical  coordinate 
system  oriented  along  the  normal  to  the  surface,  n.  Let  the  angle  between  the  normal  and 
z  be  denoted  by  without  loss  of  generality  it  can  be  assumed  that  the  normal  lies  in  the 
x-z  plane  and  that  its  projection  on  x  is  negative.  Define  an  angle  of  declination  from  the 
normal,  Xi  ^  azimuthal  angle,  which  is  measured  from  the  x  direction.  The  two 


coordinate  systems  are  then  related  by 

cos(x)  =  cos(j3)cos(^)  —  sin(<^)sin(j3)cos(a),  (A8a) 

3in(x)cos(rj))  =  cos{^)sin{<j>)  +  cos{<i>)sin{0)cos{oL),  (i486) 

sin{x)sin{tj})  =  sm(/9)sxn(a).  (i48c) 

Equations  (AS)  can  now  be  used  to  convert  (A7)  into  an  expression  in  terms  of  x> 

(^)  (A9) 

where  ^  and  0'  must  be  eliminated  in  terms  x- 


Finally,  the  angular  cross  section  can  be  obtained  by  integrating  over  the  range  of  x 
which  are  deflected  onto  the  surface.  In  addition,  dJ  is  the  differential  cross-section  per 
unit  area  in  the  x-y  plane.  In  the  paper  the  incident  flux  per  unit  area  on  the  surface  of  the 
film  is  desired;  the  ratio  of  these  areas,  which  is  given  by  the  dot  product  of  the  normal  to 
the  surface  and  z,  is  simply  cos{(f>).  Consequently 

Jv-  =  cosid)  f  dJ.  (AlO) 

dx 

In  the  remainder  of  this  appendix  the  cross  sections  for  various  particular  cases  are 
computed.  Two  types  of  surface  trapping  are  considered;  tangent  trapping,  where  only  the 
adatoms  deflected  nearly  tangent  to  the  surface  will  not  be  revaporated  and  will  ultimately 
stick,  and  adsorption  trapping  where  any  z^latom  with  a  velocity  directed  into  the  surface 
after  scattering  is  assumed  to  be  entrapped  by  the  surface  adatom  layer. 
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This  corresponds  to  the  case  to  an  infinite  mass  ratio  (M  =  oo).  From  (Al,A2)  it 
follows  that 

Ve  =  V  0'  =  20.  (All) 

Substituting  into  (A7)  yields 

dJ  =  ^3in(0‘)  da*  d0*  =  ^dS*,  (A12) 

4ir  47r 

where  dS'  is  the  element  of  surface  area.  This  implies  that  the  adatoms  aire  scattered  with 
a  constant  density  per  solid  angle.  Moreover,  the  cross-section  is  independent  of  surface 
orientation 

To  evaluate  the  integral  (AlO),  first  consider  the  tangent  trapping  case.  Note  that 
scattering  with  X  =  f  corresponds  to  leaving  tangent  to  the  surface.  Define  a  small  angle 
A,  corresponding  to  the  angular  cross-section  aroimd  x  —  f  that  are  trapped  by  the 
surface.  In  this  case  can  be  approximated  by 


(M3) 


This  allows  the  computation  of 


and  (2.4)  yields 


Jjj,  =  ^^cos{<j>) 
4?r 


Jn  =  —  COS(<^). 


(A14) 


(A15) 


In  the  absorption  case  all  adatoms  with  y  <  x  <  be  added  to  the  surface  layer. 

In  this  case, 

(A16) 


and 


Ji,  =  -—cos{<f>), 
47r 


Jn  = 


(A17) 
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2.  Two-dimensional  Hard  Sphere  Scattpring:  Jla-nprent  Traoning 
If  the  two  dimensional  hard  sphere  case  is  considered,  isotropic  scattering  is  no  longer 
seen;  this  should  be  taken  as  an  omen  that  two  dimensional  ballistic  simulations  may  have 
drastically  different  results  &om  their  three  dimensional  analogies. 

In  two  dimensions  the  differential  flux  density  is  given  by 

dj  =  ^dr  =:  cos(j3)d0.  (^18) 

Substituting  (All)  yields 

dJ  =  jco3(^)d0'.  (A19) 

Note  that  x  =  In  this  case  is  replaced  by  J+  and  J_  defined  as  the  density  of 

adatoms  scattered  forward  and  backward  repectively.  Only  the  tangent  trapping  case  will 
be  computed  here;  for  it 

J±  *  ±~cos(^  q:  |)cos(<^).  (A20) 

The  integrated  flux  is  defined  as 

/A  6 

Jn  =  J^  +  J_=  .—cos(-)co3{<f>).  (A21) 

3.  Adatom  Scattering:  Tangent  Trapping 

In  this  case  an  adatom  scattering  off  a  particle  of  equal  mass  is  considered;  that  is 
M  =  1.  Considering  relationship  (A2)  shows  that 

=  0  +  (A22) 

Consequently,  ^  <  j3'  <  ir  and  for  this  range  the  velocity  is  given  by. 


=  -C03{^'). 


Substituting  into  (A7)  yields 

dJ  =  ——co3{0')3in{x)dxd'ijj. 

TT 


(A23) 


(A24) 
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Judidoxis  application  of  equations  (AS)  yields 


cos{^')  =  cos{x)cos{d)  +  sin{x)co3{‘4>)sin{9) 


(A25) 


and  applying  (AS)  yields 


2JA 

=: - cos(il;)co3{<f>)\3in{<l>)\ 


(A26) 


To  evaluate  the  integrated  flux  (A26)  is  integrated  over  ^  for  >  f ,  to  3rield 

7„  =  2^^jsm(^)lcos<^. 

IT 


(A27) 


4.Adatoin  Scattering:  Adsorption  Trapping 
In  the  case  of  a  adsorption  trapping,  the  density  given  by  (A24)  must  be  integrated 
over  the  sector  where  /?',X  >  f  •  This  yields  after  modest  calculation 


s  Jco3{<i>) 


22M+£fM  +  «^{2*  +  «„(2«)'  . 

2  2ir 


(A2S) 


C 


r 


22 


134 


iJ-l. 


BIBLIOGRAPHY 


B.J.  Baxtholomeusz,  K-H.  Muller,  &  M.R.  Jacobson,  “Computer  Simulation  of  the  Nucle- 
ation  and  Growth  of  Optical  Coatings,”  in  Modeling  of  Optical  Thin  Filnu,  Proceedings 
of  SPIE,  20-21  August  1987,  M.R.  Jacobson,  ed.  821  (1987). 

W.K.  Burton,  N.  Cabrera,  &  F.  C.  Frank,  “The  Growth  of  Crystals  and  the  Equilibrium 
Structure  of  Their  Surfaces”,  Phil.  Trans.  Roy.  Soc.,  A23,  299-358  (1951). 

A.G.  Dirks  and  H.J.  Leamy,  “Coltminar  Microstructure  in  Vapour  Deposited  Thin  Filins,” 
Thin  Solid  Films  47  219-233  (1977). 

M.P.  Do  Carmo,  Differential  Geometry  of  Curves  and  Surfaces,  Prentice-Hall,  New  Jersey, 
1976. 


S.  Lichter  &  J.  Chen,  “Model  for  Columnar  Microstructure  of  Thin  Solid  Films”,  Phys. 
Rev.  Lett.  56,  Vol.  13,  1396-1399  (1986). 

J.  E.  Leonard- Jones,  “The  Interaction  of  Atoms  and  Molecules  with  Solid  Surfaces;  XI  - 
The  Disperal  of  Energy  from  an  Activated  Link”,  Proc.  Roy.  Soc.  A163,  127-131  (1937). 

B.A.  Movchan  &  A.V.  Demchishin,  “Study  of  the  Structure  and  Properties  of  Thick  Vac¬ 
uum  Condensates  of  Nickel,  Titanium,  Aluminum  Oxide  and  Zirconium  Dioxide”,  Fiz. 
metal,  metalloved.  28  (1969)  83-90. 

G.  B.  Whitham,  Linear  and  Nonlinear  Waves,  Wiley-Intersdence,  New  York,  1974. 


1 


135 


FIGURE  CAPTIONS 


\ 


Figure  la.  Schematic  of  the  film  surface.  A  point  on  the  surface  is  denoted  by  p.  The 
normal  to  the  surface  is  n.  The  incident  flux  J  is  assumed  uniform  and  unidirectional. 


Figure  lb.  Control  volume  on  the  film  surface.  A  flux  balance  through  the  top  of  the  control 
surface  A  and  edge  <t  is  shown.  The  normal  to  A  is  denoted  n  whereas  the  normal  along  the 
edge  <T  is  denoted  fj.  Fltix  and  velocity  along  the  surface  and  V'^  axe  defined  as  functions 
of  location  and  polar  angle  on  the  tangent  plane.  At  each  point  on  (7  a  directional  flux 
is  defined.  The  flux  in  a  direction  ^  through  a  is  given  by  (0  • 


Figure  2.  Propagation  of  a  kink  along  the  angle  bisector. 


Figure  3.  Sample  evolution  for  the  eikonal  eqxxation.  The  evolution  of  an  initial  condition 
corresponding  to  a  sinusoidal  substrate  is  shown  at  evenly  spaced  time  intervals.  The 
concavity  leads  to  the  formation  of  a  kink  type  discontinuity.  The  kink  propagates  along 
the  angle  bisector  of  the  tangent  plane  (cf.  figure  2).  Note  that  surface  roughness  decreases 
with  increasing  time. 


Figure  4.  Triple  points.  In  three  dimensions  three  kinks  can  originate  from  a  single  point. 
Near  the  triple  point  the  surface  can  be  approximated  by  the  tangent  planes  to  the  three 
grains  between  the  three  kinks.  The  motion  of  the  point  is  along  the  angle  bisector  of 
the  dihedral  angles  formed  by  the  tangent  planes.  The  networks  of  kinks  connecting  triple 
points  form  a  grain-like  structure  which  is  similar  to  that  observed  experimentally. 
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Figure  5.  Velocity  17  as  a  function  of  (f>  for  various  models.  Note  that  the  three-dimensional 
hard  sphere  model  is  identically  zero.  The  two-dimensional  hard  sphere  model  exhibits  a 
change  of  sign  not  seen  in  any  of  the  three  dimensional  models.  Presumably  the  solutions 
for  finite  mass  ratios  (M)  greater  than  unity  lie  between  the  M  =  1  and  M  =  oo  cases 
and  also  exhibit  a  single  maximum. 


Figure  6.  Soliton  propagation.  Equation  (3.36)  has  soliton  like  solutions  in  the  dispersively 
regularised  limit.  The  soliton  in  4  corresponds  to  a  step  in  surface  elevation  of  width  Li 
between  regions  of  uniform  slope.  The  step,  like  the  soliton  will  propagate  with  a  speed 
U{<l>o). 


Figure  7.  Schematic  of  a  scattering  collision.  The  incident  particle,  mi ,  with  initial  velocity 
V  strikes  the  stationary  particle  at  an  impact  angle  It  is  scattered  through  an  angle 
with  an  exit  velocity  Ve-  The  ballistic  cross*section  can  be  examined,  and  and  Ve 
can  be  computed  as  a  function  of  V,  and  the  mass  ratio  M  =  m2 /mi. 
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Analysis  and  Results  of  Existing  Instrumentation  -  W.  Wolfe 

One  aspect  of  the  project  has  been  the  development  of  instrumentation  to 
measure  the  performance  and  properties  of  thin  films.  Although  there  is  no  dearth 
of  capability  to  measure  the  spectral  distribution  of  both  the  reflectivity  and 
transmissivity  of  thin  films,  there  is  a  paucity  of  measurement  capability  for  the 
assessment  of  the  scattering  characteristics  of  these  films.  This  section  addresses 
accomplishments  and  ideas  concerning  the  three-dimensional  properties  of  scattering 
from  thin  films:  scattering  as  a  function  of  position  on  the  surface  as  well  as  depth 
effects. 

Current  Instrumentation 

First,  the  instrument  developed  provides  a  good,  relative  assessment  of  the 
scattering  characteristics  of  a  surface,  whether  it  is  a  simple  surface  or  one  with 
coatings.  Light  illuminates  the  surface  at  an  incidence  angle  of  about  80^.  It  thus 
provides  a  rather  shallow,  broad  spot  on  the  surface  of  the  sample;  there  is  relatively 
little  depth  penetration.  The  light  is  from  a  HeNe  laser  and  is  well  collimated. 
Viewing  is  normal  to  the  surface,  with  a  microscope  objective  of  numerical  aperture 
0.4  and  focal  spot  7  /im-  The  detector  behind  this  optic  is  a  silicon  photodiode. 
The  geometry  is  shown  in  Fig.  1.  The  spot  can  be  scanned  over  a  total  field  of  1x1 
mm^  by  controlling  the  motion  of  the  sample,  which  is  situated  on  an  x-y  stage,  using 
a  computer.  The  illuminated  spot  is  actually  much  larger  than  the  viewing  spot. 

What  does  this  system  actually  measure  and  what  information  does  it  provide? 

The  viewing  system  senses  all  the  light  scattered  toward  it  in  its  solid-angle 
field  of  view.  This  is  a  combination  of  the  light  scattered  from  the  surface  and  the 
illuminated  sub-surface  "slice"  that  is  within  the  field  of  the  view  of  the  receiver. 
Although  the  interpretation  is  not  straightforward,  the  patterns  of  relative  intensity 
are  very  nice. 

Shown  in  Fig.  2  are  the  results  of  a  spatial  scan  on  a  section  of  Halon.  Halon 
is  a  material  that  is  used  as  a  photometric  reference  standard  in  the  visible  spectral 
region.  It  is  almost  isotropic  and  has  a  hemispherical  reflectivity  of  almost  one.  It 
accomplishes  these  two  feats  by  being  rough  (with  respect  to  wavelength)  and  non- 
absorbing.  The  pattern  shown  indicates  a  measured  intensity  ratio  of  about  S:l  over 
the  scanned  field  of  0. 1x0.1  mm.  This  spatial  pattern  is  not  inconsistent  with  a 
Lambertian  reflector.  On  a  scale  of  7  pm  the  surface  is  inhomogeneous.  The 
surface  is  rough  with  respect  to  the  wavelength,  and  it  does  reflect  more  or  less  in  a 
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particular  direction  for  each  small  area.  The  BRDF  curve,  shown  in  Fig.  3,  reveals 
that  these  variations  must  average  to  give  a  uniform  pattern. 

Shown  in  Fig.  4  is  the  pattern  of  a  highly  reflecting,  specular  surface.  It  too 
varies  from  spot  to  spot.  The  scanning  scatterometer  also  provides  a  reasonable 
representation  of  its  homogeneity  (or  lack  thereof).  A  BRDF  plot  is  given  in  Fig.  5. 

It  is  reasonable  to  ask  whether  a  BRDF  measurement  could  be  used  to  calibrate 
the  scanning  scatterometer  (SS).  TTie  answer  is.  "Yes.  but.  ..."  We  compare  the 
SS  to  the  AZSCAT  (Arizona  Scatterometer).  Both  can  be  made  to  have  th^.  same 
incidence  angle  and  wavelength.  However,  the  AZSCAT  collects  radiation  from  a 
1-cm  diameter  spot  and  l-jiisr  field,  whereas  the  SS  collects  from  a  7-^m  diameter 
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spot  and  0.3-sr  field.  To  generate  equivalence,  the  SS  results  would  have  to  be 
integrated  over  more  than  its  whole  field  of  view,  and  the  AZSCAT  results  would 
have  to  be  integrated  over  one-third  its  total  range.  This  is  more  integration  than 
has  been  done  in  the  United  States  since  the  Civil  Warl 

So,  AZSCAT  measures  angle-resolved  scattering  over  relatively  large  spots  and 
the  SS  measures  angle-integrated  scattering  from  small  spots.  Each  is  useful;  each  is 
different;  they  are  complementary  probes. 

Of  course  it  is  interesting  to  ask  whether  there  is  a  way  to  calibrate,  regularize, 
rationalize  the  two  instruments  so  that  a  calibrated  BRDF  spatial  scan  could  be  made. 
It  seems  clear  that  one  should  make  the  spot  of  the  AZSCAT  as  small  as  that  of  the 
SS  or  make  the  angle  of  the  SS  as  small  as  that  of  AZSCAT— or  at  least  make  the 
two  approach  each  other.  There  are  two  issues  involved.  One  is  radiometric 
sensitivity;  the  other  is  instrumental  capability. 

The  fundamental  equation  for  the  noise  equivalent  BRDF  in  this  situation  is 


where 


R» 


NEBRDF  - 


A„B 


EiAoA,cos0i0,D’ 


R  -  distance  from  the  sample  to  the  detector 
Ajj  -  detector  area 
B  ••  noise  bandwidth 
Ej  -  sample  irradiance 
A,  -  observation  area 
Af  -  receiver  area 
-  angle  of  incidence 
0,  -  angle  of  scatter 
D*  -  specific  detectivity 

With  appropriate  implementation,  the  bandwidth  B  can  be  1  Hz  and  the  detector 
area  A^  can  be  20  fim  on  a  side  (4x10'^  cm^).  A  reasonable  value  for  the  distance 
R  is  20  cm  and  for  the  laser  power  5  mW.  The  D*  of  a  good  silicon  detector  is 
about  10‘*  cmHz‘/*W"‘.  So  most  of  the  values  are  fixed,  chosen  or  determinable. 
For  these  values 


NEBRDF 

I 


(400  cm^X20xlO~*  cmXl  Hz) 

(Ej Ag ArCos0jCos0j)x  1 0‘*  cmHz‘/*W“‘ 

0.8x1 0~‘^  cm^W 
EjAgArCoseiSj 
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If  the  viewing  is  done  normally,  then  dj  -  d,.  If  it  is  60°  then  cosd  -  0.5.  The 
irradiance  will  be  the  power  divided  by  the  beam  area  divided  by  the  cosine  of  the 
incidence  angle.  Thus 

Ab(1.6xl0-*°  cm°W) 

NEBRDF  -  — — n-^r— j - 

Pi^oAr 

If  one  now  requires  that  the  NEBRDF  is  lO^sr"*.  then  must  be  I0"*cm'^. 

The  product  of  the  observation  area  and  the  collecting  area  is  a  constant.  One  can 
back  off  some  and  require  that 

10-^ 


A,Q  - 


400 


0.25x10^  cm^r 


The  system  throughput  is  a  constant.  For  10  ;im  and  0.3  sr,  this  is  3xl0~^  or 
0.3x10^.  These  are  approximately  the  parameters  of  the  SS.  For  the  AZSCAT.  the 
parameters  are  approximately  1  cm  and  10“*sr.  The  AZSCAT  could  reduce  its 
observation  area,  but  not  by  much.  So  neither  instrument  can  be  improved  by 
much— except  by  greatly  increasing  the  observation  time. 

It  would  seem  from  this  analysis  that  there  is  no  way  to  properly  compare  the 
two  instruments.  The  SS  measures  the  scattering  from  a  small  spot  over  a  large 
angular  area.  The  AZSCAT  measures  the  scattering  from  a  large  spot  over  a  small 
angular  area.  For  BRDFs  on  the  order  of  10“°sr"*  in  the  visible,  neither  can  be 
improved  markedly.  These  are  the  BRDFs  of  interest.  For  larger  ones,  comparisons 
can  be  made  and  measured.  For  instance,  a  BRDF  of  10"°sr“*  .somewhat  lower  than 
a  perfect  isotropic  reflector  at  0.3sr~*  could  be  accomplished  by  decreasing  the 
collection  angle  of  the  SS  with  an  aperture  stop,  and  increasing  the  collection  angle 
of  the  AZSCAT  with  a  collection  lens.  Other  such  maneuvers  are  possible. 


Instrumentation  for  Bulk  and  Surface  Measurements 

In  the  beginning  there  was  light,  and  specular  reflection  measurements  were 
made.  But  man  grew  sophisticated  and  wished  to  know  how  light  scatters  from 
materials.  And  now  he  would  like  to  measure  how  much  scatters  from  the  surface 
and  how  much  from  the  volume.  Not  an  easy  question!  (But  it  is  easier  than  how 
this  all  began.) 

"Traditional"  scatter  measurements  are  made  by  shining  light  on  a  transparent 
surface  and  viewing  the  light  that  is  either  back  scattered  or  forward  scattered  from 
that  "spot."  Of  course  every  optical  system  views  a  solid  angle,  so  the  scattered 
light  is  not  from  the  spot  but  from  the  entire  solid  angle  that  is  viewed.  This  means 
that  most  measurements  of  BTDF  include  surface  scatter,  subsurface  scatter  and 


volume  scatter.  Most  BRDF  measurements  of  reflecting,  opaque  surfaces  include 
surface  and  some  subsurface  scattering.  Most  BRDF  measurements  of  transparent 
samples  include  volume  scattering  as  well  as  the  scattering  from  both  surfaces. 
Separation  of  these  effects  is  possible  to  some  extent,  but  it  is  not  easy. 

Consider  the  geometry  of  Fig.  6.  The  beam  is  incident  at  some  carefully 
chosen,  arbitrary  angle.  The  light  that  gets  to  the  aperture  scatters  from  the  front 
surface  over  an  area  B1  with  a  certain  BRDF.  It  scatters  as  well  from  the  sample 
volume  between  B1  and  B2  in  the  beam.  It  also  scatters  from  the  back  surface  B2. 
These  are  the  primary  scatterers.  Next,  the  scatter  from  B1  is  specularly  reflected 
from  the  back  surface  and  scattered  from  the  back  surface.  The  incident  beam  is 
scattered  from  the  volume  to  the  back  surface  and  reflected  up.  These  are 
secondary  scatterers.  Tertiary  ones  include  »:atter  from  the  front  to  the  back  and 
back  and  then  up;  from  the  volume  to  the  back  and  up;  from  the  back  through  the 
volume  and  up,  back  to  front  and  up.  It  is  hoped  that  higher-order  scattering  can 
be  neglected.  The  message  here  is  that  scattering  in  a  non-absorbing  sample  is  a 
highly  complex  matter.  In  the  diagram  shown,  the  flux  on  the  aperture  comes  from 
the  first-surface  BRDF  at  one  angle,  second-surface  BRDF  at  another,  volume  scatter 
at  an  intermediate  angle,  and  then  there  are  second-order  effects.  Emphasis  can  be 
applied  to  front-surface,  back-surface,  or  volume  scattering  depending  upon  the 
geometry,  as  shown  in  Figs.  7,  8  and  6  respectively. 
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Figure  6.  A  possible  scattering  geometry. 


One  approach  to  modeling  all  this,  and  therefore  to  measuring  it,  is  to  calculate 
the  BRDF  from  the  first  surface  based  on  the  well-known  formulae  and  the 
geometry,  add  to  this  the  BRDF  from  the  second  surface— with  a  slightly  different 
angle  (based  on  Snell's  Law),  and  add  an  appropriate  amount  of  surface  scatter.  All 
these  must  be  weighted  by  the  solid  angles  of  collection.  Several  measurements  at 
different  incidence  and  scatter  angles  must  be  made  to  separate  the  effects. 

A  second  approach,  that  may  be  used  to  complement  the  first,  is  the  use  of 
frustrated  total  internal  reflection.  The  frustration  in  this  case  is  the  result  of 
scatter,  rather  than  the  close  juxtaposition  of  another  material.  The  geometry  of  the 
radiometry  is  critical.  The  collimated  beam  is  arranged  to  exit  the  sample  at  a 
critical  angle. 

The  specular  beam  is  then  totally  reflected  by  a  perfectly  smooth  surface.  A 
rough  surface  will  scatter  some  of  the  beam  out  of  the  surface— according  to  the 
well-known  equations  of  the  vector  EM  Theory.  Some  light  is  also  scattered  by  the 
volume  into  angles  that  are  less  than  critical.  By  careful  focusing  and  direction  of 
the  field  of  view  of  the  receiver,  the  scattering  from  the  near-surface  structure  can 
be  emphasized.  The  entire  secret  seems  to  be  to  consider  as  carefully  as  possible  the 
focusing  geometries  of  both  beams.  One  approach  is  shown  in  Fig.  7,  where  the 
collimated  beam  is  viewed  as  it  would  strike  and  where  it  doesn't.  If  both  beams 
are  focused  on  the  surface,  the  emphasis  is  on  the  surface  scatter.  It  may  not 
completely  separate  surface  from  volume  scatter,  but  it  should  emphasize  one  over  the 
other  to  an  extent  that  significant  measurements  can  be  made.  The  geometry  is 
shown  in  Fig.  9. 

A  third  approach  is  that  of  nephelometry.  The  input  beam  and  viewing  field 
are  arranged  such  that  the  beam  strikes  the  surfaces  only  outside  the  field  of  view. 
The  main  scatter  will  then  be  that  from  the  volume,  but  it  will  be  modified  by  the 
surface.  It  should  also  be  possible  to  move  the  field  to  measure  the  surface  to  obtain 
correlative  data.  This  is  shown  in  Fig.  6,  where  neither  surface  is  viewed. 
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Abstract 

We  demonstrate,  by  the  use  of  second-harmonic  generation,  that  a  significant  aniso¬ 
tropy  in  Langmuir-Blodgett  monolayers  may  be  imposed  by  the  deposition  process.  This 
represents  the  first  demonstration  of  a  systematic  influence  on  the  properties  of  a  pre¬ 
pared  monolayer  due  to  the  type  of  film  balance  employed. 
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Monolayer  studies  with  cyanine  dyes  at  the  water-air  interface  were  of  great  inter¬ 
est  in  the  past  because  of  their  technical  importance  as  sensitizers  in  photography. 
Recently,  they  have  been  used  more  generally  as  model  systems  for  the  study  of  the 
crystallization  and  self-organization  processes  in  quasi  two-dimensional  systems.*** 
Moreover,  these  monomolecutar  layers  are  the  basic  units  for  the  build-up  of  multilayer 
assemblies  by  the  Langmuir-Blodgett  dipping  technique*  aimed  at  tailoring,  at  the  molec¬ 
ular  level,  functional  superstructures*  with  special  properties.  Thus  techniques  are  nec¬ 
essary  that  allow  for  a  sensitive  characterization  of  the  geometric  as  well  as  electronic 
structure  of  mono-  and  multilayers.  Second-harmonic  generation  (SHG).  with  its  intrin¬ 
sic  surface  and  thin  film  sensitivity,  has  been  shown  to  fulfill  these  requirements.  Pre¬ 
vious  workers  have  derived  information  about  molecular  hyperpolarizabilities  or  orien¬ 
tational  properties  of  single  molecules  or  parts  of  the  molecules.**** 

We  will  show  in  this  letter  that  structural  information  is  also  available  at  the 
supramolecular  level  of  single  crystallites  and  demonstrate  this  by  elucidating  the  slight 
anisotropy  in  the  orientation  of  elliptically  shaped  monomolecular  crystallites,*  so-called 
J-aggregates,'**'*  on  the  water  surface.  We  show  that  this  anisotropy  is  maintained 
during  the  transfer  process  and  is  dependent  on  the  type  of  Langmuir  film  balance 
employed  for  the  preparation  of  the  monolayer. 

The  cyanine  dye  used  is  S120.  (l-methyl-r-octadecyl-2,2'-cyanine  iodide,  structural 
formula  given  in  Fig.  1.  kindly  provided  by  Dr.  K.  Wirthensohn.  University  of  Munich) 
which  can  be  spread,  compressed,  and  thereby  organized  at  the  water-air  interface.’* 
The  break  in  the  pressure-area  isotherm  (Fig.  I)  marks  the  onset  of  the  formation  of 
crystalline,  highly  fluorescing  aggregates  which  are  induced  by  increasing  the  lateral 
pressure  r  to  above  S  mN  m'*.  These  J-”  or  Scheibe‘*-aggregates  can  be  grown  to 
dimensions  of  10  to  100  /im.  A  typical  ensemble  of  elliptically  shaped,  randomly  ori¬ 
ented  crystallites  is  shown  in  the  inset  of  Fig.  I.  This  fluorescence-microscopic  picture 
was  taken  from  the  water  surface  at  ir  a  IS  mN  m'*  with  polarized  light.  Only  those 
aggregates  with  their  transition  dipole  moment  parallel  to  the  E-field  fluoresce  brightly. 
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For  the  preparation  of  the  J-aggregate  monolayers,  two  differently  designed  Lang¬ 
muir  troughs  were  used.  The  first  one  was  a  commercial  (Joyce-LoebI,  JL  for  short) 
tank  equipped  with  a  Wiihelmi  plate  for  the  pressure  reading.  A  special  feature  of  this 
system  is  the  half-immersed  teflon-coated  glass-fiber  ribbon  which  confines  the  spread 
molecules  to  a  well-defined  area  (Fig.  2a).  Compression  is  achieved  at  constant  perime¬ 
ter  by  two  counterpropagating  overarms  which  leave  the  center  line  of  the  working  area 
still  (see  Fig.  2a).  The  second  trough  is  home  built  (HB)  and  equipped  with  an  integ¬ 
rated  fluorescence  microscope."  Here,  compression  is  performed  by  a  single  moving 
barrier  which  slides  over  the  edges  of  the  teflon  frame  (Fig.  2b).  The  resulting  com¬ 
pression  speed  profile  across  the  monolayer  is  substantially  different  from  the  first  one 
(schematically  indicated  in  both  figures  by  the  length  of  the  arrows)  with  considerable 
shear  stress  along  the  long  edges  of  the  frame.  Fluorescence  microscopy  during  com¬ 
pression  shows  that  this  stress  relaxes  stochastically  in  the  densely-packed  J-aggregate 
layers  which  are  very  stiff.  The  result  is  a  slight  positionally  dependent  orientational 
redistribution  of  the  crystallites. 

The  monolayers  were  transferred  onto  hydrophobic  substrates  as  follows:  First, 
two  clean  BK-7  slides  were  mounted  together  with  edges  sealed  with  Teflon  ribbon  in 
order  to  prevent  subphase  liquid  from  penetrating  between  the  plates.  This  assembly 
was  then  coated  with  one  cadmium  arachidate  monolayer  by  the  usual  vertical  dipping 
technique.  Dye  monolayer  transfer  onto  the  front  side  was  subsequently  performed  on 
the  down  stroke  -  either  by  hand  or  by  a  mechanical  sample  lift  with  the  substrate 
being  mounted  nearly  horizontally.  The  angle  a  between  an  arbitrarily  chosen  axis  in 
the  plane  of  the  substrate  and  the  center  line  of  the  monolayer  could  be  chosen  freely 
(see  Figs.  2)  to  an  accuracy  of  ±5”.  While  submerged  the  substrate  was  rotated  to  the 
vertical  and  the  water  surface  aspirated.  A  protective  overlayer  of  arachidic  acid  was 
added  in  a  conventional  upstroke.  Finally  the  two  slides  were  separated  and  only  the 
front  one  was  used  for  the  experiments.  This  procedure  ensured  the  coating  of  only 
one  side  of  the  substrate,  with  the  back  side  being  absolutely  free  of  contamination. 
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The  laser  at  frequency  u  was  incident  on  the  center  of  the  substrate  at  an  angle  6 
relative  to  the  surface  normal,  with  the  dye  monolayer  on  the  second  surface  of  the 
glass  slide.  The  intersection  with  the  beam  is  adjusted  to  insure  that  a  fixed  spot  on 
the  sample  is  illuminated  when  the  sample  is  rotated  around  the  normal  to  the  surface 
by  an  angle  ^.**'**  The  second-harmonic  (at  2w)  generated  by  the  film  travels  colline- 
arly  with  the  laser  beam  and  is  detected  in  the  forward  direction.  The  laser  beam  con¬ 
sists  of  3S  ps  pulses  at  10  pulses  per  second  at  a  wavelength  of  1.064  nm.  Typical 
experimental  pulse  energies  were  1  mJ  spread  over  a  measured  beam  diameter  (1/e*)  of 
2.3  mm.  Polarizers  allowed  for  a  free  choice  of  $-  or  ^waves  for  both  frequencies,  u 
and  2u. 

We  describe  the  SHG  process  by  a  nonlinear  polarization.  P2«.  associated  with  the 
monolayer,  induced  by  the  electric  field  of  the  laser.  E„: 

Plu  -  .  (1) 

The  form  of  is  determined  by  the  crystal  symmetry  of  the  J-aggregates.  For  a 
single  oriented  J-aggregate.  the  second-harmonic  signal  may  then  be  calculated  in  the 
usual  way.*'  In  fact,  fluorescence  microscopy  shows  that  the  film  consists  of  a  macros¬ 
copic  distribution  of  J-aggregates  (see  inset  of  Fig.  1).  each  of  which  is  oriented  at  an 
angle  relative  to  some  9  of  the  film  as  a  whole.  The  laser  beam  samples  ailO* 
J-aggregates  whose  dimensions  are  large  compared  with  the  second-harmonic  wave¬ 
length.  We  therefore  average  the  second-harmonic  signal’  from  a  single  oriented  J- 
aggregate  over  the  macroscopic  distribution  function,  p(^i).  Here  is  the  number  of 
aggregates  oriented  at  ^i,  ner  unit  tj.  Thus  the  total  second-harmonic  signal  is  given 
by 

Si*)  -  p(*,)  S"««(*  -  ♦,)d*,  .  (2) 

Jo 

where  is  the  second-harmonic  signal  from  an  individual  J-aggregate  oriented  at  <>  - 
*i  relative  to  the  laboratory  frame  ^  -  0  axis.  We  now  expand  both  the  distribution 
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function  p(^,).  and  S'"'®  in  a  Fourier  series  and  obtain 


OO 

P(*i)  -  ^  PmCos(mdi  +  «m)  (3) 

ni"0 


and 


Smic(«  -  ^  S„cos[n(*  -  ♦,)  +  «J  .  (4) 

n<-0 

Generally  only  a  few  terms  in  this  last  series  are  required  to  specify  Sub¬ 

stituting  Eqs.  (3)  and  (4)  into  Eq.  (2).  performing  the  integration  and  setting  5o  -  eg  -  0 
gives 

OO 

S(t)  -  »  •  pgSg  +  »  XI  +  *m  +  <m)-  (5) 

m-0 

Gearly  S(4)  will  show  an  angular  dependence  containing  only  those  spatial  frequencies 
present  in  both  S'"^(4)  and  p(^,).  Since  we  expect  that  S*"^  from  a  single  J-aggregate 
has  lower  than  isotropic  symmetry,  we  can  deduce  from  S($)  information  about  the 
macroscopic  distribution  function  p(4|).  Alternatively,  given  p(^i).  it  is  possible  to  det¬ 
ermine  the  microscopic  symmetry  class  of  the  individual  J-aggregate.*" 

The  dependence  of  the  second-harmonic  signal  observed  under  pj„  Po„t  condi¬ 
tions  (i.e..  incident  laser  p  polarized,  analyzer  set  to  record  only  ^polarized  second-har¬ 
monic  output),  from  a  dye  layer  deposited  in  the  JL  trough  with  a  -  0".  is  shown  in 
Fig.  3a.  Gearly.  the  sample  is  not  macroscopically  isotropic,  but  shows  a  two-fold  rota¬ 
tional  symmetry.  This  result  is  somewhat  surprising,  as  it  has  generally  been  believed 
that  molecules  on  the  water  surface  of  a  Langmuir-trough  are  isotropically  distributed. 
We.  on  the  other  hand,  are  dealing  here  -  at  least  in  the  late  stage  of  the  isothermal 
compression  -  with  large,  elliptically  shaped  crystallites.  It  is  therefore  well  conceivable 
that  the  two  moving  belts  apply  an  uniaxial  force  to  the  aggregates  which  leads  to  sene 
orientation  upon  compression.  To  test  this,  and  to  rule  out  that  the  observed  anisotropy 
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in  the  macroscopic  distribution  function  was  induced  during  the  transfer  process,  we 
deposited  another  film  on  a  similar  substrate,  but  this  time  with  the  slide  rotated  to 
approximately  or  •  4S*  in  the  plane  of  the  water  surface  (horizontal  dipping  is  employed 
throughout)  so  that  the  compression  axis  was  rotated  4S‘  relative  to  a  marked  edge  of 
the  substrate.  The  SHG  measurements  of  this  film  are  shown  in  Fig.  3b.  As  expected 
for  our  interpretation,  the  origin  of  the  curve  is  shifted  by  approximately  45*  and  the 
direction  of  the  shift  is  correct.  The  slight  assymmetries  are  attributed  to  a  residual 
lateral  inhomogeneity,  which  can  be  directly  observed  by  scanning  the  laser  beam  across 
the  sample. 

When  the  monolayer  is  prepared  on  the  HB  trough  with  the  single  compressing  bar¬ 
rier.  Fig.  4  shows  that  for  a  >  0  the  twofold  symmetry  of  the  distribution  function  is 
again  found  under  pj^  ■*  conditions.  However  there  is  a  loss  of  the  correlation 
between  the  azimuth  dipping  angle  a  and  the  phase  angle  (5g,  +  e„)  of  the  d-scan.  S(d). 
In  order  to  ensure  that  the  observed  phase  angles  are  not  influencevi  i<\  the  hand-dip¬ 
ping  process,  a  mechanical  sample  lift  was  employed  which  could  be  i><>Miioned  prior  to 
the  monolayer  deposition  to  a  much  better  anpilar  accuracy.  Despite  the  obvious  loss 
in  phase  correlation,  these  samples  seem  to  be  more  isotropic  in  their  macroscopic  orien¬ 
tation  distribution  as  judged  from  the  considerably  reduced  modulation.  Also  the  lat¬ 
eral  homogeneity  was  significantly  better,  as  checked  by  a  separate  experiment.  As 
demonstrated  by  the  full  lines  in  Fig.  4,  the  angular  dependence  can  be  described  by  a 
single  cos  24  function,  with  negligible  improvement  by  adding  some  cos  44  contribu¬ 
tions  (not  shown). 

In  concluding  we  note  that  for  most  applications  of  LB  films,  (poly-)  crystalline 
monolayers  of  different  materials  are  grown  on  the  water  surface  and  then  transferred 
to  the  solid  support.  We  have,  we  believe  for  the  first  time,  demonstrated  a  systematic 
influence  of  the  type  of  film  balance  employed  on  the  properties  of  the  prepared  mono- 
layer.  And  we  have  shown  that  SHG  is  a  sensitive  tool  to  monitor  even  subtle  aniso¬ 
tropies  in  the  macroscopic  orientation  distribution  of  polycrystalline  monolayers. 
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Figure  Captions 


Fig.  1.  Pressure-area  curve  of  a  cyanine  dye  (SI 20,  structure  formula  given  in  the 
inset)  monolayer  on  pure  water  (Millipore  quality.  pH  5.S.  T  ■  22*’Q.  The 
onset  of  J-aggregate  formation  is  indicated  by  the  arrow.  The  inset  shows  a 
fluorescence  microscopic  picture  taken  from  the  water  surface  with  polar¬ 
ized  excitation  at  a  pressure  ir  -  15  mN  m~^  The  area  shown  corresponds 
to  75x100  /tm. 

Fig.  2a.  Schematic  top  view  of  the  Joyce-Loebl  Langmuir  trough  used  for  some  of 
the  J-aggregate  monolayer  preparations.  Special  features  are  the  two  coun- 
terpropagating  overarms  which  control  the  active  area  by  means  of  the 
Teflon-coated  glass-fiber  ribbon  at  constant  perimeter.  The  arrows  indicate 
the  resulting  compression  speed  profile  on  the  monolayer  with  the  charac¬ 
teristic  center  line  of  speed  zero.  There  the  (near)  horizontal  transfer  of  the 
condensed  monolayers  (at  v  -  40  mN  m*')  onto  the  hydrophobic  substrate 
was  performed  by  hand  at  various  angles  a. 

Tie.  2b.  Schematic  top  view  of  the  home  made  Langmuir  trough  which  was 
equipped  with  an  integrated  fluorescence  microscope.  Compression  is  per¬ 
formed  with  a  single  moving  Teflon  barrier  gliding  over  the  edges  of  the 
teflon  frame.  The  resulting  compression  speed  profile  with  its  shear  stress 
at  the  edges  is  again  indicated  by  the  arrows.  Deposition  of  the  monolayer 
onto  the  substrate  was  performed  by  a  mechanical  sample  lift. 

'r^.  3.  Angular  4  dependence  of  the  second-harmonic  intensity  of  a  SI  20  J-aggre¬ 
gate  monolayer  transferred  form  the  Joyce-Loebl  trough  at  a  -  0°  (a)  and  ct 
-  45®  (b).  Angle  of  incidence  35®,  laser  pulse  energy  1.5  mJ;  ^polarization 
for  u>  and  2tj)  radiation.  Note  the  corresponding  phase  shift  of  the  <^-scan. 
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Same  as  in  Fig.  3.  but  for  a  monolayer  prepared  in  the  home-built  Lang¬ 
muir  trough.  For  (a),  a  -  0"  and  the  modulation  can  be  fit  (solid  line)  to  the 
single  cosine  function  1.53  +  0.18  cos  (2<ft  -  32°).  For  (b).  a  m  45°  and  the 
full  line  corresponds  again  to  a  single  cosine  function  1.41  +0.15  cos  (2$  - 
157°)).  Note  the  uncorrelated  phase  of  the  <^-scan. 
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